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Introduction 


SCOLE stands for the "Spacecraft Control Laboratory 
Experiment". The objective of the SCOLE Program is to 
provide an example configuration and control objectives 
which enables direct comparison of different techniques in 
modeling, systems identification and control. The "SCOLE 
Design Challenge" was formulated in 1983 by L. W. Taylor and 
A. V. Balakrishnan. The details of this challenge are reprinted 
at the end of this document. 

Annual SCOLE Workshops have been held for specialists 
to share and compare their research results. This proceedings 
is a compilation of the material presented at the 5th 
Workshop held at Hilton Lodge at Lake Arrowhead, California 
on October 31, 1988. 
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DYNAMIC ANALYSIS OF THE JOINT 
DOMINATED BEAM (TRUSS BEAM) 

by 

Elias G. Abu-Saba 
Raymond C. Montgomery 
William M. McGinley 
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INTRODUCTION 


Construction in Outer Space of: 
Platforms for Space Stations 
Antenna Systems for Space Explorations 
Support Structures for Solar Panels 
Housing for Space Workers 


( 2 ) 



ABSTRACT 


Method Presents 


- Theoretical Analysis of the Vibrational 
Modes of the Joint Dominated Beam 


“ Cantilever Truss Beam is Used for the 
Analysis 


- Chord Members Contribute Most 
Deformations. Web Members are 
Ignored 


- Lumped Mass System is Used for 
Analysis 
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- Algorithms are Developed to determine 
Flexibility of the System 


- System is Analysed With and Without 
Joint Contribution 


- A Set of Joint Flexibility is Used 


- Computer Programs are Developed to 
Obtain Numerical Results: Frequencies 
and Mode Shapes 


- Conclusions and Recommendations are 
Provided 
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DYNAMICS 

EQUATION OF MOTION OF SYSTEM 
[B]{Y}= 0 

[B] = [A] - _i_[I] 

m w i 


3 3 

A i i = 1 ( 2i-l ) , i = 1 N 

24EI 


3 2 

Ri = _1 3( 2i -1 ) , i = 1 N 

24EI 


Aij = Ai,j-1 + Ri , j = i + 1 , . . . ,N 
Ai j = Aj i 


( 5 ) 




JOINT PARTICIPATION 


Assumptions 

Behavior of the Joint is Linear 
Damping of the Joint is Not Included 
Free Play is Not Considered 
Forces in Members are Axial 
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MODIFIED FLEXIBILITY MATRIX 

3 2 

A*ij = Aij +__J_ s(i) (2 j- 1 ) 

24EI 

i = 1 N 

j= 1 N 


s - 1 2 E I k 
2 

h 1 


k = Flexibility of the Joint 
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FREQUENCY ( RflD/SEC ) 


FREQUENCY CORRELATION 
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CONCLUSIONS 


A Practical and Simple Approach for 
Predicting Frequency Response of 
Structures 


Computational Cost, Time and Storage 
Requirements Provide a Clear 
Advantage over Finite Element Method 
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EXPERIMENTAL PROGRAM 


STAGE 1 Static Loading 


1 . Determine the flexibility of the 
joint assembly. 

2. Determine the combined flexibilty of 
the truss and joint assembly for 
increasing numbers of truss panels. 

3. Under the action of static loading 
measure the stress distribution in the 
truss panels. 
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EXPERIMENTAL PROGRAM 


Stage 2 Dynamic Response 


1 . Measure the deflection and the 

truss cord stress variations for a single 
truss panel and joint assembly. 

2. Measure the dynamic response 
(deflections and stresses) of an 
increasing number of truss panel 
and joint assemblies. Frequency 
response will be determined from 
these measurements. 


( 11 ) 
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DESIGN Oh' THE SCOi.E BOOM 

BASED ON THE DYNAMIC ANALYSIS Ob' THE JOINT DOMINATED 

BEAM 


Elias G. Abu-Saba 

North CaroJina A&T State University 

Raymond C. Montgomery 
NASA Langley Research Center 

William M. McGin.ley 
North Carolina A&T State University 
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ABSTRACT 


Mathematical models have been developed to predict the 
dynamic response of the joint dominated beam. Various 
assumptions of the force interaction between the beam 
elements and the joints have been studied using these models. 
However, the validity of these models have not been 
adequately determined. In any dynamic analysis of 
assemblies, the effects of joint imperfections must be also 
included. To accomplish these tasks, a combined analytical 
and experimental investigation will be conducted to evaluate 
the significance of the various joint behaviours with respect 
to the overall stiffness of the beam. Results of the 
experimental analysis will be incorporated in the dynamic 
analysis of the beam. 

Models of the truss-beam will be constructed at North 
Carolina A&T State University and tested to determine the 
validity of the analytical method. Test data will be 
obtained for a particular joint type with negligible free 
play. Using the measured material properties of the joints 
and truss model, the equation of motion will then be used to 
predict the dynamic behaviour of the truss-beam. A 
comparison of predicted and observed behaviour will be made 
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and the adequacy of the analytical methods evaluated. 
Modifications of the analytical procedures will be made as 
required by the comparison study. This approach of the 
dynamic study of the truss-beam system can be used in the 
dynamic analysis of the scole. 
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INTRODUCTION 


Manned space stations are on the drawing board, and 
structures to house astronauts and support personnel will be 
erected in outer space. These stations will be used tor 
communication, power generation, docking, and launching space 
vehicles in the future. Since all structural components have 
to be fabricated and preassembled on earth and packaged for 
shipment via the space shuttle, weight and volume become a 
determining factor in the design. Light weight materials as 
epoxy graphites are used in the members, and interlocking 
joints are provided to facilitate deployment. 

As the exploration of outer space is increased, the need 
for more reliable structures which perform various tasks in 
orbit becomes essential. Since these structures are 
deployable, large, and highly flexible, their dynamic 
analysis and control using established methods yields results 
that differ from those obtained on board testing. Hence 
further modifications of the mathematical model will be 
necessary to increase convergence between theoretical and 
experimental results. 

DYNAMIC EQUATION 

The general equation of motion for a linear structural 
system is expressed as 

2 

C Y } = w L'AJCMDCYJ (1) 
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where 


{ Y ) - Displacement vector 

C a 3 = Flexibility matrix ot the structure 

CM3 = Mass matrix of the structure 

w = Natural frequency of the system 

rad/ sec 

The dynamic study of the truss-beam is seen to depend on 
its f lexiblitity . 

FLEXIBILITY OF THE TRUSS BEAM 

Consider that the truss-beam behaves as a cantilever as 
shown in Figure (1). For a slender truss the transverse 
nodal displacements will be mainly attributed to the axial 



FIG. I - TRUSS IbE. AM 

strain in the chords. Thus the contribution from the web 
members can be neglected with very little error. If the joint 
response were to be excluded innitially, the flexibility 
matrix CAD for the structure can be readily obtained from 
the following algorithms: 
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3 

1 3 

A = ( 2i - 1 ) , i = 1, N 

ii 24EI 

3 

1 2 

R = 3 ( 2i - 1 ) , i = 1, N (2) 

i 24E1 


A = A + K , j = i + 1 , N 

i _i i j - 1 i 


A A 

i j i i 


When the bays are identical, the nodal masses will be 
equal throughout the structure. Letting m represent the mass 
of one bay, the mass matrix can be expressed by: 

CM3 = m C 1 3 < 3) 

Substituting Equaion 3 into Equation 1 and simplifying 
the result yields 

LBJCY) = 0 (4) 


where 


L B 3 


L A 3 


XC 1 3 


(*>) 
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24EI 


( 6 ) 


7 \ = 

2 3 

raw 1 

Given the material property E, the cross sectional area 
of the members of the truss-beam, the depth of the bay h, and 
the length 1, the number of bays N and the density of the 
material used in the structure, a computer program has been 
written to provide the natural frequencies of the system. 

FLEXIBILITY OF THE TRUSS BEAM WITH JOINT RESPONSE 

To simplify the response of the joint, it is assumed that 
the joint displacements are caused by the strain in the chord 
members. Another assumption which is used in the study is 
that the displacement in the joint is due to axial forces in 
the chords acting on the connecting pin as a spring. 'The 
joint is thus replaced by a flexibility k. No other factors 
that contribute to joint imperfections are included at this 
stage . 

Introducing the joint flexibility k, the modified 
flexibility matrix is given by 

3 

* s 1 2 

A = A + i ( 2 j - 1 ) (7) 

ij ij 24EI 

i = 1, N; j = 1, N 
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where A is given by Equation 2 and s is obtained from 
i D 

the following expression: 

12 El k 


APPLICATION TO THE METHOD 

An example is used with a typical panel of 20" in length 
and 20" in depth. The cross sectional area of the chord is 
0.25 sq. in. The values of E and the density are 3E+7 psi and 
0.282 pci, respectively. The finite element method is 
applied to the same example for comparison. Since the finite 
element method does not include joint imperfections, the 
value of k is taken as zero. The results of this comparison 
are shown in table 1 and Figure 2 for N varying from i to 80. 

EXPERIMENTAL PROGRAM 

To confirm the accuracy of the theoretical procedure, an 
experimental approach has been proposed to be performed at 
North Carolina Agricultural and Technical State University. 
The proposal has been submitted to NASA at Langley Research 
Center for funding. 

Since the accuracy of the analysis is highly dependent on 
the mass and flexibility of the truss panels, in addition to 
the flexibility of the joint, these quantities must be 
accurately measured before a reliable evaluation of the 
analytical procedure can be made. 


Table 1 - Natural Frequencies, Rad/sec 
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FREQUENCY (RADS /SEC) 
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Only one type of "joint will be evaluated during" the 
experimental phase. Furthermore, the joints will be 
fabricated so that free play in the joint is negligible 
and can be ignored. Once the accuracy of the analytical 
procedure have been continued for this type of joint, 
subsequent investigations can determine the effects of 
varying joint flexibility and free play. 

TESTING APPARATUS 

Truss joint specimen will be tested using a Forney 
Material Testing machine located in the Structures Laboratory 
at the North Carolina Agricultural and Technical State 
University in Greensboro, North Carolina. 

In an effort to reduce cost, the same apparatus will be 
used to test the truss panels and truss panel assemblies. A 
schematic of the proposed testing apparatus is shown in 
Figure 3. It consists of a truss support frame, air cushion 
table, actuator, and an adjustable actuator support and 
"hold-down". The apparatus will also include devices to 
measure the deformation and chord strains at various 
positions on the specimen, in conjunction with a high speed 
data acquisition system and a micro-computer. 
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REMARKS AND CONCLUSIONS 


The approximation of a large space truss by a beam model 
provides a practical approach for predicting the frequency 
response of the structure. The computational cost, time and 
storage requirements may render discrete finite element 
analysis impractical, particularly in the conceptual and 
design phases when parameter studies and alternate designs 
are being evaluated. 

The proposed beam method offers a structural model with 
significantly fewer elements (N) and degrees of freedom (N>, 
compared to a finite element model with (4N+1) degrees of 
freedom, where N is the number of panels in the actual 
structure. Frequencies for the first five modes obtained by 
the beam model are found in good agreement with those 
obtained by the finite element, especially as N gets large. 
The proposed algorithm saves significant programming effort 
as well as results in considerable economy of computat ional 
time, cost and storage. 
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Nonlinear Damping Model: Response to 
Random Excitation 

Weijian Zhang* 

Electrical Engineering Department 
6731 Boelter Hall 
UCLA 

Los Angeles, CA 90024 
November 10, 1988 


Abstract 

The objective of this study is to investigate stationary Fokker-Planck equa- 
tion corresponding to nonlinear random vibration problem. A method of en- 
ergy approximation (MEA) is proposed to obtain a modified model as an ap- 
proximation of the exact model. The closed form solution of the stationary 
Fokker-Planck equation corresponding to the modified model is obtained, and 
so are the various moments of the stationary response. Comparisons are made 
between MEA and MEL (method of equivalent linearization) to illustrate the 
advantage of MEA over MEL. The MEA also overcomes the shortcoming of 
non-Gaussian closure method in which the density might have negative value 
caused by truncation. 


* Research supported in part under AFOSR grant No. 83-0318. Presented at the 5th Annual 
NASA SCOLE Workshop, Lake Arrowhead, California. Nov., 1988. 
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1 Introduction 


Experimental evidence in SCOLE [l] seems to support the need for nonlinear 
models-the decrement is much smaller than predicted by linear models. Nonlin- 
ear damping model under random excitation becomes a random vibration problem 
which has received significant attention in recent years in both civil and mechanical 
engineering. When the excitation is a white noise process, the response is Marko- 
vian. The transition probability density function, which together with an initial 
distribution completely describes the response, satisfies the Fokker-Planck equa- 
tion which is a linear but degenerate parabolic type partial differential equation. 
A variety of approximate methods have been proposed to solve the Fokker-Planck 
equation corresponding to broad-band excitation of a nonlinear dynamic system. 
All of the methods investigated so far appear to involve prodigious amounts of 
labor. 

In this paper, a Method of Energy Approximation (ME A) is proposed to find 
the approximate solution of the stationary Fokker-Planck equation corresponding 
to one-dimensional nonlinear damping model under white noise excitation. 

2 Energy Approximation Model 

Let us first consider a concrete model due to Balakrishnan and Taylor [2]: 

x + 2£ijJqX + 7x 2m |x|"i 2n+1 |i| /? + w 2 x = crn(t) (l) 

where 0 < a, (3 < 1, £, 7 > 0, and m,n are nonnegative integers. n(t) is a Gaussian 
white noise process. 
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Notations: 


E = 


u>$x 2 + y 2 
2 


y = 

q = 


2 

m + n + 


a + (3 
2 


We observe the following: 


72 2m |2| a i 2n+1 |i;|* 5 


= 7|z 


2m+a U/l 2n +^, 




y| ""y 

) 

^ /r*77» ,,2\m+§ /n 77» ,.,2_2\n+4. 


2m+ 
^0 


-(2E-y 2 r + H2E-u,lx 2 ) n+ ?y 


We consider the following replacement of the above: 


7 

,2m+a 

W 0 


(2£)<y 


Here fi 0 is chosen such that 



\[2E 


sm ‘ip\ 2m+a \\/2E cos rp\ 2n+0 \/2E cos ip — fj, 0 (2E) q y/2E cos %p\ 2 dip 


is minimized. 

In the above minimization, we consider E as a constant because in situations 
where the response has reached stationarity and where the damping is small, the 
energy dissipated through damping and the energy input through the excitation 
during one cycle, will, on average, be small fractions of the total energy level in 
that cycle. 


29 



Computation gives: 


Ho = — f s'm 2m+a xjj cos 2n+2+ P ipdxjj 

7T Jo 

2 r(m + g f 1 )r(u + 1 + 

7r r(g + 2) 

Then, we obtained the following energy approximation model: 

x + 2£u 0 x H — "Im+a (^o X<1 + ^ 2 ) q x + u l x = on{t) (2) 

u> 0 

When there is no random noise excitation, computer simulation shows the 
Krylov-Bogoliubov approximation of (l) is a very good approximation when 7 
is small. It happens that (l) and (2) have exactly the same Krylov-Bogoliubov 
approximation, which is an evidence to support the validity of the energy approxi- 
mation model. Of course, this fact only referrs to no noise models. 


Proposition 1 Both the exact model (1) and the modified model (2) have the same 
Krylov-Bogoliubov approximation given by 


da(t) 

dt 

d<fi(t) 

dt 


-£w 0 a(t) - ^Wo” +/ V 9+1 (0 
0 


( 3 ) 


Proof: Straightforward computation based upon Krylov-Bogoliubov approximation 

[ 5 ]: 


da(t) 

dt 

d<f>{t) 

dt 


1 [ 2 * 

/ D(a sin ip,aojo cos ip) cos ipdtp 

27tl> 0 Jo 

1 [ 2 * 

— / D( a sin ip, aoj 0 cos tp) sin ipdip 

2ttu>o a Jo 
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3 Stationary Fokker-Planck Equation 


In this section, we solve the stationary Fokker-Planck equation corresponding to 
the energy approximation model. The technique is used in [3]. 

The stationary Fokker-Planck equation is 


0 = -T x [yp) ^Yy Hxp) 


+ ^:K 2 ^o V + (2 E) q y)p} + 


dy 


Wn 


d ,dE , d ,dE , 
~ dx^ dy P + dy^dx^ 


9 [(2£w 0 y + ^j^{2 E) q y)p + 


dy 


Wn 


2 3y J 


0 = y[(2£w 0 + w 2^°i'(2^) 3 )p+ y^] 


In the above, we used the following fact 
p{x,y) = p[E) => 


-£(!!*) + £(ffri = ° 

— yi£- 


dy 


dE 


Therefore, the stationary density is given by 

2£u> 0 


p{x,y) = C Q exp[- 
= C 0 exp[- 


( 2 E) - 7 ! ( 2 E) 9+1 


2£u>o 


(wqX 2 + y 2 ) - 'ni^x 2 + y 2 ) ,+l1 


where 


7Mo 

71 ~ a 2 Wo 2m+ “(y + l) 


C, 


L = — f c e -^r'-™"‘dp 

7n CJn «/0 


7T 


— ^(«o) 


Wo7i 


<7 + 1 
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F{s) = f°° e~ st e~ tq+l dt 

Jo 

2£u> 0 

s 0 — ~T~ 

ct 2 7 i ,+1 


4 Moments of Stationary Response 


~Ex k y l = 0 if either k or l is odd. If both k and l are even, then 

1 


k n l 


Ex y 


;E(w 0 x) V 


Wn 


c„r(^)r(^) 


x p^ L exp[-‘—^-p - hp' 1+l }dp 

r(¥)r(¥) 


M S(¥)i 


"Ti 


k+l 
2 « + 2 


*■(*>) 


5 Comparison with MEL (Mehtod of Equivalent 
Linearization) 

Consider £ < 0 case of (l). The Krylov-Bogoliubov approximation of (1) is given 
by (3), which can be solved as 

a(0)e^ u ' ot 


a(t) = 


2n+0-l X 

1 + a 2 «(0)2^b (e 2 !^ 4 - 1)]*? 


(4) 


We have 


lim a(f) 


m 

2 |£| 


*- °° ' ' ‘7MoWo n+ ^ 1 


J- def _ 

2« = a 


for any a(0) ^ 0. This fact implies that the exact model (l) has a stable limit cycle 
with radius between a and u 0 a (approximately). 
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(5) 


Actually, the modified model (2) does have a limit cycle 

f x 0 (i) = asinw 0 £ 
y°(t) = uoacoswot 

or, on phase plane (§) 2 + (^) 2 = 1. 

Proposition 2 The limit cycle (5) of the energy approximation model (2) is asymp- 
totically stable. 

Proof: To prove the asympotic stability of the limit cycle (5), we consider the 
corresponding linear variational equation 

( dx - y 

{ I = £(x°(t),!/ 0 (<»z + f(x°M (6) 

where 

Y{x,y) = - 2 £ Wo y - + y 2 ) q y - "0* 

u> Q 

The linear system (6) with periodic coefficients always has a non-trivial ~ 
periodic solution which is given by (i 0 (i),y°(i)). Therefore/?! = 1 is a characteristic 
multiplier. By Theorem 2.1 of [4, page 217], it is sufficient to show the other 
characteristic multiplier p 2 < 1. 

Since 

r 2z. ay 

P2 = Pi • P‘1 = ex P [jJ° -^-(^ O (0»2/ O (*)) rff ] 
it is sufficient to show 

r 2* ay 

f U ° ^-{x°{t),y°(t))dt < 0 
Jo dy 

In fact, by noticing the definition of a we have 

^(x°(t),y°(*)) = -2^o-^oWo” + V 9 (l + 2 9 cos 2 u;oO 

dy 

= -2£u> 0 - 2|£|w 0 (l + 2gcos 2 w o 0 

= — 4y|^|w 0 cos 2 u 0 t 
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and hence 


[ — 4 g | fiascos 2 u 0 tdt = — 47 rq\£\ < 0 

Jo 

which establishes the asymptotic stability of ( 5 ). □ 

What is more interesting is that the stationary solution of the corresponding 
Fokker-Planck equation given by 

p(x, y) = C„ e xp[^(u>’x J + y 2 ) - 7,kV + y ! ) ,+1 ] (7) 

achieves maximum on the ellipse (limit cycle). 

(~) 2 + (~~) 2 = 1 
fl W 0 G 

and achieves local minimum at the origin. 

The following Van-der-Pol self-excited oscillator 

x + (7 x 2 — r])x + x — on{t) (8) 

is a particular case of the above discussion with £ = *7/2 (*7 > 0), m = 1, n 

a = p — 0 , dq — 1 . 

One can easily find n 0 = 1 / 4 , a = 2^77/7 311(1 the ener gY approximation model 
is given by 

x + [7/4 (x 2 + x 2 ) - rj]x -t- x = an(t) 
with the corresponding stationary density 

, ..2\ _J_f^ _i_ ,,2^,21 


where 


p(x,y) = C 0 exp[^(x 2 + y 2 ) - + y)] 


— = + *{-%)] 
Co V X/T 


$(x) = f \J2jn j* e~^^dt 


J 4 


6 Generalization 


The MEA described above can be easily generalized to nonlinear random vibration 
problem in the following general form 

x + D(x, x) + WqX = on(t) (9) 


Our conclusion is that the exact model (9) can be approximated by the following 
modified model 

x + //( — --- - -- )x + UqX = an(t) (10) 

where fi(E) minimizes 

r 2w \/2 E 

/ [Z?( sin ip, \/2E cos ip) — p(E)y2E cos ip] 2 dip 

Jo UJq 

and is given by 

4 f — \/2E 

p(E) = — f D{ sin ip, \/22£cos ip) cos ipdip (11) 

7 t\ 2E Jo c j 0 

Proposition 3 If D(x,y) is even with respect x and odd with respect to y, then both 
the exact model (9) and the modified model (10) have the same Krylov-Bogoliubov 
approximation. 


Proof: For the exact model (9), 

da(t ) _ 1 r 2 * 


dt 27TW, 

By the assumption on D(x,y), 

d<p{t) 


1 [ 2ir 

— / D(a sin ip , au>o cos ip) cos ipdip (12) 

ru)Q Jo 


dt 


1 f 2 * 

/ D( a sin ip, au>o cos ip) smipdip 

2n<jJoa Jo 

1 rf 

/ D( a cos ip, auo sin ip) cos ip dip 

TTU! 0 a J-i 



0 



For the modified model (10), the Krylov-Bogoliubov approximation is given by 

MU = cos 1 ^ 

dt 2nu>o Jo 2 

= 

2^' 2 ' 


1 /' 2?r 

— / D (a sin ip, Uoa cos ip) cos ip dip 

rwn Jo 


27TCJ0 


d(f>(t ) _ 1 

dt 2Tvu>o a 

= 0 


- — [ p(JJ^-)u> 0 a cos ip sin ip dip 

Una Jo 2 


Therefore, (9) and (10) have the same Krylov-Bogoliubov approximation. □ 


The solution of the stationary Fokker-Planck equation corresponding to (10) is 

given by 

2 r" oW 



p(x,y)=Cexp[ 2 p{z)dz\ 

(J J 0 

(13) 

where 

l^rex P [--J\{z)dz)d P 

C uq Jo a 1 Jo 

(14) 

Example: 

Consider the following saturation type active damping model 



x + 2£u> 0 x + A tan -1 (6i) + ujqX = on(t) 

(15) 


By the identity 
/•*Y 2 


r * ’/2 7f 

/ tan -1 ( b cos x ) cos xdx = — - 

Jo 2 


/ITT*- l 


beR 1 


one can compute 


4 

TTy/2E 


2£u>o + A 


- JT 

/ 2 [2fw 0 \/2£ cos ip + A tan _1 (6\/2i? cos 0)] cos ip dip 
Jo 

. v/I + We - i 
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And consequently, one has 


2 f E 

2 / n{z)dz 

o L Jo 

= ~^{2E) - ^Vi + WE 

a 2 o L b 

+-^- ln[l + Vl + 2b*E] + const. 

0*0 

The solution of the stationary Fokker-Planck equation: 
p{x,y ) = C\ 1 + y/l + t 2 (wqX 2 + y 2 )]^ 
x exp[— ( w o x2 + y 2 ) _ ^\A + b^(u lx 2 + y 2 )] 

where 

1 7T , 2fu>o, 

c = Z^ exp{ ^W } 

f°° . r 2£o> 0 4A 

x J 1 (1 + \ffi^exp\-—x-—y/Z\dx 

It is easy to realize that p(x, y ) achieves maximum at the origin. 


7 Concluding Remarks 

A Method of Energy Approximation(MEA) is proposed for the investigation of 
nonlinear damping problem under random excitation. Closed form steady state 
density is obtained for the energy approximation model. It is shown through an 
example that MEA gives better accuracy than MEL and MEA reflects the nonlin- 
ear nature in the damping. This problem is pending further research. Problems 
such as spectral density of the nonlinear damping model, the absolute difference 
between the exact steady state density and the density obtained by MEA are under 
investigation by the author. 
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Abstract 


A new, explicit representation of damping operators for strictly proportional damping 
for the torsion mode of a finite beam with end mass is presented. The damping operator is 
the square root of the stiffness operator (enhanced to include the boundaiy) and is calcu- 
lated using the Balakrishnan formula. It is nonlocal, and turns out to be a finite-limit 
version of the Hilbert transform for the clamped-clamped case. If strict proportionality is 
required, the operator is more general and involves boundary terms which, however, tend to 
zero as mode frequency increases. 
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1. Introduction 


In the design of active controllers for stability agumentation of flexible structures it is 
naturally important to have a model for the inherent passive damping already present. Often 
it is assumed that the damping in each mode is strictly proportional to the mode frequency. 
If, as in [1], we formulate the problem in the abstract wave equation form as: 

Mx{t) + Ax(r) + 2^Dx(t) + Bu c {t) = 0 , 

for strictly proportional damping we will need to have 

D = 

where 

t = Va/" 1 a Vat 1 . 

Usually M commutes with A so that 

D = <M<A . 

In any event we are left with the problem of calculating 'lA . In this, the first part of a 
two-part paper we calculate the square root for the case of a beam torsion model — or the 
familiar “string” equation. When the controls are on the boundary — as in current large 
space structure control design [2] — the square root introduces terms on the boundary. 
This should certainly be considered much more unnatural than the fact that the damping 
operators are nonlocal ( not differential operators). Asymptotically however the modes do 
approach those of the clamped-clamped (or fixed-ends for the string case) beam. In 
particular for the clamped case the square root has the form of a “finite-limit” Hilbert 
transform. Our point of departure is the fractional-power formulas due to the author [3]. 
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2. Beam Torsion Model 


We begin with the torsional mode of vibration of a uniform Bernoulli beam of finite 
length, with one end fixed and the other end with control force. Thus we have: 


pa - G/ v = 0, 0 < r; -Z < s < Z 


u{t,-Z) = 0, 0 < r 


(2.1) 


GfyuXt.l) + / 4 ii(r, £) + u c (t) = 0 


where for our purposes we need only note that pa, Gly and I 4 are given positive constants. 
The primes denote derivative with respect to s, and the superdots derivative with respect to r. 
The abstract formulation is obtained by taking 


X = U[-t, Z] x R l . 


Let us use the following notation for elements in X: 


w € H ; w 



«(•) g LA.-Z, £], b <= R l . 


We shall shorten LA-Z, Z] to L 2 in the sequel. Define the operator A with the domain in 
X given by: 


2)(A) = f w- U( ' } , «"(•) ls 2 , u(-Z ) = 0 1 

u(Z) J 


(where primes denote derivatives with respect to the space variable), and 


Aw = 


-«"(•) 

u'(Z) 


mapping 2)(4) into X. 

Then it is readily verified that A is closed, self-adjoint and, [ , ] denoting inner-product 
in X: 

1 t 

[Aw, w] = J -u"(s) u(s) ds + u(Z)u'(Z) = J u\s ) 2 ds . 

-e -t 
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Hence A is nonnegative definite. Also A has a compact (Hilbert-Schnudt) resolvent and zero is 
in the resolvent set. As in [1] we may reformulate (2.1) as an abstract wave-equation in H: 


Mx(t) + Ax(t) + Bu c (t) = 0 


( 2 . 2 ) 


where B maps /?' into H by 


Bu = w = 


0 

u/GI 


v 


Mw = 


(pa/G/ y )u() 

(/ 4 /G/ v )b 


(2.3) 


Our first objective is to calculate the (positive) square root of A denoted 'Ia . For this 
purpose it is helpful to use the notation 


R(X, -A) = (XI + A) 

so that we can write the Balakrishnan formula [3] for the square root: 

VAw = 1 f°°X~ I/2 (XJ + A)~ 1 Aw dx, for w in 3) (A). (2.4) 

*8 

It is convenient now to introduce the operator a D with domain in La with 
SDfao ) = [«(•)! u"()eL2, «(-£) = 0 = «(£)] 


and 


c^u - -u". 

Then a 0 is closed, self-adjoint and nonnegative definite and hence we may proceed to define 

by 


w J ^ V2 (M + °o) 1 q o U ^ 
n 0 


(2.5) 


Our first result is: 


Theorem 2.1. 


5D(VO = [«(•) I «'(•) e W, «(-«) = 0 and “(*) = °1 


( 2 . 6 ) 



and 


u = v 


has the representation 


^ = h 


cos (rc.y/2£) 


sin (7W/2£) - sin(7tc/2£) 


-u'(a) da , a.e. -I < s < l 


(2.7) 


Proof. We use (2.5). Let u e 3)(a 0 ) and let: 

u 0 {l) = QJ + 0 ^,) ck) u • 
Then « 0 (^) is the unique solution of 


Xuo (X, s) - Uq(K, 5) = -u"(s) 1 


U 0 (K, -£) = uo(k, «) = 0 


and is given by: 


m 0 (X, s) 


f sinh V/((£- s ) sinh V\ (£+c) ^ u"(q)) Jo 

J VX sinh 2VX £ 

-r 


+ 


f sinh VX(£ + .s) sinh VX(£-q) , ^a 

J VX sinh 2VX £ 


Integration by parts yields 


) _ r’ sinh VE(t- si cosh ^t±gl ( „ (g) ) * 

J sinh 2VA, £ 

_t (2.8) 

r* sinh VX(£ + s) cosh VI(£-q) (u^q)) do. 

J sinh 2VX £ 


Using the known formula: [6, p. 344] 
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J 


sinh ax cosh bx , 

r— r ax 

sinh cx 


jz_ sin (a7t/c) 

2c cos ( aide ) + cos (bn/c) ’ 


c > a + b , 


we can see that for c < s 


If ,- 1/2 sinh s} cosh V\(£+ a) ,, 

*J sinh 2'JX Z 


2 f sinh x(t-.r) cosh *(£+0) , 

71 J sinh 2 x 1 

0 

(2.9) 

1 COS (7t$/2£) 

21 sin (ks/2Z) - sin (KG/2Z) 


and for s < a 


X f ,-1/2 sinh VX(£ + 5) cosh VX(£-q) « 

*J sinh 2VX Z 


2 f sinh *(£+$) cosh x(£-o) , 

7t J sinh 2x1 

0 


1 cos (tls/2£) 

21 sin (ns/2Z) - sin (kg/21) 


Hence for a.e. in (-£, £) 


— f x 1/2 u(X, j) dx 

K i 


J_ f cos (7U/2£) , . 

2Z J sin (ks/2Z) - sin(7to/2£) “ w 


( 2 . 10 ) 


providing of course that we can justify the change in order of integration. Let 


H f (s, a) 


_1_ 

21 


r cos ks/2Z 

.sin ks/2Z - sin no/2l 


—Z < s,o < Z . (2.11) 


and let 


t 

v F (s) = J H F (s, a)u'(a) da 

-t 

where «(•) e X(a 0 ). We can decompose H F (s, a) as 
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H f (s,o) = ^cot C) 


— tan K(5 + q ) 
41 1 4 i 


For any «(•) in Lj , 


i 

v i( J ) = ilj ( cot TC( 4£~ °^] “(°) da ’ -t<s<l 


the integral taken in the Cauchy sense at s = a, is defined wherever the Hilbert transform 


if _«£) 

71 J (s - 


is defined, and hence we note that vi(-) is actually also in Z^(-£, £). Hence 


- k\ 0) u ' (0) da 


which is defined in the open interval -t < s < t, must be in Lz, if Vi(-) + v 2 (-) is. 
Next, for a finite L, we can clearly change the order of integration, and hence let: 


. L t 

v L (s) = — J «o(X, s)X~ l/2 dk = J H(L, s, a)u'(o) da . 


2 f si 

s, S, O) = I 


sinh x(£-s) cosh x(fc+a) 


sinh 2x1 


2 f si 

= 


sinh *(£+$) cosh x(l-o) 


sinh 2 xt 


so that 


v F (s) - v L (s) = J R(L, s, a)u’(s) ds . 


Now for x > L and L sufficiently large, 


sinh x(£-,r) cosh x(£+o) _ -*(j-o) 

sinh 2x1 


4 $ 



sinh x(£+s) cosh x(g-0) 
sinh 2 xt 


a > s . 


-x(a-s) 

€ 9 


Hence 


R(L, s, a) 


2 

JC 


-L(j -a) 

€ (s- a) 


2 

71 


-L(o-j) 


j > a 


a > s 


and hence 

s, O)u'(o) do = J* R(L,s,C)u'(o)da + J R(L, s, G)u'(O) do . 

-I |s — O |> £ |j-Oj<£ 

The first term goes to zero as L -> for any e > 0. The second term can be written 

2 f E wTr+x) - k'(s-x) -Li dx 

K J X 

0 

and since «(•) e £(a 0 ), 

X 

u'(s + x) - «'(5 - X) = J u"(5 + o) da 

-X 

goes to zero with Vt . Hence v L (s) converges to v^-Cr), a.e. in (-£, £). But v L (-) converges 
to in Li and hence 

V/r(‘) = ^ u , 

or we have justified changing the order of integration. Of course, we have also proved that 
\ F (•) defined by (2.12) is actually in L for u in 50(0,,). 

Let us next establish (2.6). The orthonormalized eigenfunctions of a Q are 

1 k%(s + t) , > „ _ o 

<M*) = ^ sm 2 1 ’ -t<s<l 

and 

^ <t>* = §f <t>* • 


Hence for X > 0, 
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(Va^" + Xl)g = f 


has the solution 



[/. 4>J 

\ kK 

K+ 21 




where the series is absolutely convergent. Hence 

8(i) = 0. 

Also term-wise differentiation is valid, showing rb.ii g'(‘) G ^ 2 - 

Since the domain of is precisely the range of the resolvent (A, + Va^)” 1 , it follows 
that for g in the domain of we must have that 

g(D = 0 ; g ' (•) g . 

Conversely for any g in with these properties we can find a sequence u n in X(a 0 ) such 
that converges to g'(-) in Lj , and u„ converges to g. We have only to take, for 
example 

u n = n(n + Oo) - g = I rkn.'\^ • 

1 " + I2IJ 


Then 


Ih^oK - W » +P )ll 2 = Ihi -“•♦pll 2 • 

Hence u„ converges. Hence it follows that g g jD(Vc^~ ), or (2.6) is proved. Also 

= v » 


where 

1 

V n (s) = J H F (s,o) u ;(o)do 

-t 

= iej 001 n( , 4t 0) “" (c) *** “ 130 ^4£ + “" (G) 

-t -« 

where the first term converges in and the second term for each \s\ < Z. Hence it follows 



that if 


v = limit v n , 
i 

v(s) = J H F (s y a)g'( a ) do y a.e., -Z < s < Z . 

- £ 

Hence for every g in SD(Va^) we have the representation above. 

Let us note in conclusion that 

i 

v(s) = f H F (Sy o)u(o) do 

-t 

is defined a.e. in (-£, £) for every u(-) in Z^, but v F (*) is in if (and i n fact only if) 

e 

J u(o) do - 0 . 

-e 

Finally, for £=+<», it is known (see [4, 5] for example) that is Hu where H is the 

Hilbert transform. Hence (2.12) may be viewed as the Finite-limit Hilbert transform, new with 
this paper. 

We can now get back to Va . 


Theorem 2.2. 


r «(•) 

1 

w: 

, u (•) e L2, uc-e) - oj 

eW 


and we have the representation: 

r~ r- w ( ‘ ) 

<Aw = ^Ja = D 0 w + 7w (2.12; 

<Z) 

where 7 is a compact linear bounded operator on X into X, and Dq has the same domain 
as VZ and is closed thereon, and has the form: 


Luv) + u m , 

D 0 w = 

I [hy u] I 

where L is linear-bounded on Z^ into L 2 , and /i(-) e L 2 , <}>(*) e Z^ 
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Proof. We use (2.5). Let we D(A); and 


w 


«(•) 

u(t) 


Let 


w(X) = (X/ + A) 1 Aw , X > 0 


Then 


is the unique solution of 


w(X) = 


m(X, •) 
M(X, i) 


X«(X, s) - u"(k, s) = -«"(*) , 
Xm(X, t) + u'(k,l) = «'(*) 


-t < s < i 


u(X, -£) = 0 . 


Let, as in the proof of Theorem 2.1, 


-l 


and let 


«o(X, •) = (X + c^) a 0 u (-) , 


z(X, s) = u(X, j) - «o(X, s) 


Then z(X, •) is the unique solution of 

Xz(X, s) - z"(K,s) = 0 
z(X, -£) = 0 

Xz(X, t) + Z'(X, l) = U'(£) - u 0 '(X, *) 


Hence 


z(X, 5) = O(X) sinh VX (£ + s) 

<j(X)(X sinh 2^X£ + VX cosh 2VX £) = «'(£) - “o(X, £) 


(2.13) 


so 



But 


•*) = 


_ f' ^ cosh sft(t- s) cosh 'l\(t+o) Ja 

J sinh 2 VX i 

-( 

cosh VX(e + s) cosh U ' (CT) 

J sinh ifk Z 


Hence 



cosh VX (£ + o) 

(sinh 2 VXE)(VX sinh 2 VX£ + cosh 2 VX£) 


h'(O) do . 


( 2 . 14 ) 


Hence 


-f r W 2 z(X, 5) dCk 

*o 

= IfY 1/2 <ft I"' Sinh C ° Sh 5 L - + g) 7= — «'(C) do ( 2 . 15 ) 

71 q J (sinh 2 'TkZ)('JX sinh 2 VX£ + cosh 2 VXt) 

* ;I (sinh 2rf)Csinh ft Z cosh 5t) ( j'^ *1”) l«*» »°)W°) *■ 

0 ( 2 . 16 ) 

and the second integral in parenthesis can be expressed 

. t , t 

2 J (sinh x( 2 £+s + o))u'(o) do + (sinh x(s-o))u\o) do . ( 2 . 17 ) 


The second integral in ( 2 . 17 ) can be integrated by parts to yield 


1 lr 

= 2 (sinh x(s- £))«(£) + 2J x ( cos ^ *( 5_0 )) u (°) do . 


Hence ( 2 . 15 ) can be expressed 

oo l 

— f — ■ — r — ~ — . - ~~r . . .. f (sinh x( 2 £+s + a))u'(°) do 

tt J (sinh 2 jc£)(x sinh 2 x£ + cosh 2 x£) ■ J v 

0 ( 2 . 18 ) 

t 

+ J M(s,o)u{o)do + <>(s)u(fc) 

-t 


FJ 



where 


m - ij 


sinh x(.y-£) 1 

sinh 2x1 (xsinh 2xZ + cosh 2xt) 


dx, -Z<s<Z (2.19) 


and is an absolutely continuous function which goes to zero as Z and the kernel M( -, •) 
is given by 

,,, . If * cosh xj s - a) ^ (2.20) 

M(s,o) - jjJ ( sinh 2 xZ)(x sinh 2xZ + cosh 2 xZ) 
o 

and is absolutely continuous in -Z < s,o < Z, (and also goes to zero as Z ->°°). Hence 


— f X 1 12 u(X, s ) dk 
*o 

t 

= <t>(j)«(e) + f M(s, a)u(o) do 

-i 

oo £ 

.if — — — — r- f sinh x(2Z + s + o)u'(o) do 

+ n J (sinh 2x£)(xsinh 2 xZ + cosh 2 xZ) 
o 

+ -f dx sinh ^ f cosh x(Z + o)u\o) do 

7t J sinh 2xZ J 

0 

_ 2 f”^ sinh x(Z + s) r * cosh x(Z - o)u'(o) do . 

u J sinh 2xZ J 

o 

We now combine the last three terms and then justify interchanging the order of integra- 
tion. The sum of the last three terms can be expressed as: 

~ t 

-If . r 0 pr f sinh x(2l + 5 + o)u'(c) do 

“ kJ (sinh 2x£)(xsinh 2x1 + cosh 2x1) ^ 

i 

J sinh x(s + C)u\o) do 

-t 

s 

J sinh x(2 Z + o- s)u\0 ) do 

-t 

e 

J sinh x(2 Z + s - o)u’(o) do . 

S 

S2 


l 

dx 

K , 
i 

sinh 2xZ 

) 

1 

oo 

dx 

K * 

sinh 2 xi 

0 

1 1 

oo 

dx 

K J 

sinh 2 xZ 


The sum of the first two terms 


1 P dx * 

v 3 (*) " n) (* sinh 2x£ + cosh 2x£) J ( cosh “ *sinh **+o))ii'(o) da (2.21) 

0 “ c 

For w(-) in 5D(/4), w"(-) is in Z^, and hence we can change the order of integration in the 
last two terms and express their sum as: 

V4 (-s) = J u'(cr) cot n ( s 4 ~ do , -l < s < t . (2.22) 

-t 

Defined wherever the Hilbert transform is, 

v(i) = J U(C) cot do , -l < s < i (2.23) 

-t 

yields a bounded linear transformation mapping L 2 into Li . It follows in particular that 
v 3 (-) is in L 2 . Let 


K(L, s, o) 



cosh xfo + q) - 
xsinh 2xt 


0 


xsinh jc( 5 + q) 
+ cosh 2x1 


dx 


and 


m 


if f cosh t_x - xsinh tx ) , 

nj [x sinh 2xt + cosh 2xt) ’ 
0 


-2i<t<U. (2.24) 


Then 


t 

J K(L, s, o)u'(o) do , 
-* 


converges in Li to v 3 (-) as L->°°. Now 


-£ < s < t 


\K(s + o)| = 0 ( 2e _ j _ ~ ) , (s+c) < 21 


and hence the integral 


S3 



J* K(s + a)u'(a) da 


IS 


defined for each s , s * £. Also for some M < 


t K 

J \K(s+a) - K{L, s, o)| |u'(o)| da < m\ 

-i -e 


* -L(2t-s-cs) 


21 - s - a 


l«'(o)l da 


-U 2e-j) t 


21- s 


M J |m'(o)| da 


-i 


and hence 

t 

v 2 (s) = J K(s + a)u'(a ' la , 
-t 

Hence finally, we have, for w in 3)(A): 

y[A w = 'Ja 

where 


«(•) 


v() 

«(£) 


v(£) 


where 


v(s) = Vi(s) + V 2 (s) + v 3 (s) 

l 

vi(s) = ^(5)H(t) + J M(s, a)u(a) da 

-t 

V 2 W = J u '(o) cot 

-e 

e 

v 3 (j) = J + o) «'(<*) ■ 


-* 


Next let us calculate 


-t < ^ £ 


J A, 1/2 u(A., £) = J X 1/2 z(X, £) <&, (since uo(M) = 0) 


-\{*l 

0 -i 


“ s .n do 


X sinh 2xt + cosh 2x1 


(2.25) 


(2.26) 




Now 


t 

J cosh x(l + o)u'(a) do 

-i 


< 11“ '(Oil 'sj J t cosh 2 ar(^+o) do 


t 

< J cosh x(l + o)|u'(o)| do 

-i 


and 


t 

J cosh 2 x(£ + o) do 

-t 


sinh 4x1 (. 2x1 1 

2x [ + sinh 4xij 


Hence (2.26) is smaller in absolute value than 


Vsinh 4xZ -j 1+ . . 

V2x (x sinh 2x1 + cosh 2x1) 

o 


and because of the presence now of >/x in the denominator, the integral is readily verified 
to be finite. Hence we can change the order of integration in (2.26) and express it as 

t 

J h(s)u\s) ds (2.27) 

-t 

where h(-) is defined for -i < s < t by 



cosh x(£ + s} 

x sinh 2xt + cosh 2xt 


dx 


(2.28) 


and is positive for every s. Moreover since 



£ ||u'C)ll • const. 


it is clear that 


t 

J h(s)u(s) ds 
-t 




defines a continuous linear functional on Li and in particular 


J h(s ) 2 ds < °° . 

-t 

Define the operators on Li into Li : 

v - Tu ; v(j) = Ms)u(l) + J M(s 9 c)u(o) do. 

-t 

and 

e 

v — Hu j v(j) = J M(o) cot da , a.e. -£<*<£, 

-t 

where the integral is defined in the Cauchy sense at s - a. Then for w in J)(A): 


ylA 


w = 


v(-) 

I. AC0«'C0 ds 




where 

v(s) = V](j) + v 2 (s) + v 3 (s) , a.e.. 


where 

Vi = Tu 

v 2 = Hu’ 

l 

v 3 (s) = J K(s + a)u\a) da 

-t 

Next we let w € 3 )(V^). Then we can find w n € 3)(A) such that 

w n w 


and 


"4X w n -* w . 



(2.29) 


(2.30) 


(2.31) 


Let 


\\<Aw n - <Aw m \ | 2 = \\u‘ -u’\\ 2 


and hence u„'(-) converges in L 2 to u\ where 


and u' e Li . 


Tu„ + Hu ; -> Tu + Hu'. 


Hence 


v 3.n( J ) = J AT(s + a)u„'(a) da , -i<s<t 


is such that 


V 3.n() 

converges in L 2 . But for each s, s < i, we do have that 
t t 

f K(s + CT)u n '(o) da -> J K{s + a)u\a) da . 

-t -t 


Hence 


* 

J K(s + a)u'(a) da , 


is in L 2 . Now given any z(-) in L 2 , we can define 


S 

u(s) = J z(o) da 


J K(s + o)z(a) da = f K(s + a)u'(a) da 


Hence it follows that 
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Ku = v ; 


(2.32) 


t 

v(s) = J K(s + o) u(a) da 

-t 

actually defines a linear bounded operator on La into La- 
Next let us consider the sequence 

t 

J h(s)u;(s) ds . 

-t 

This sequence converges to 

[h, u '] 


since u ' n converges to 

Let us next show that any w in X of the form 


w = 


«(*) 

u&) 


U(-t) = 


where u'(-) g La belongs to 2 )(Va). Then we can find u n in La 


«,(-«) = o 

«„'(•) € La 


and such that «„'(•) converges to u‘(- ) in La . Let 




«»() 

M*> 


Then w n g I (.4) and 


-> w . 


Also 

t 

||VAw n - VAwJ| 2 = f \u^(s) - u^(s)\ 2 ds 

-t 

and hence 

Va w n converges . 


0, 

such that 
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Hence w e SD(>/i4 ). 

Thus finally for any w 


in 3)('lA) 


w 


we have the representation 

«(•) 

«(*) 


where 


'lA w 


v() 
[h, «'] 


= Hu' + tfu' + 7u. 


Lu‘ + «(£)$ 
[h, u'] 


We can decompose VA w in various ways. For instance let L = H + K, and 

Dqw = 

on 3 )(VA ) . Then D 0 is closed thereon. Let 

Mu 

Tw = 

0 


where 

t 

Mu = v ; v(s) = J M(j, o)u(o) da . 

-t 

and M(-, •) is given by (2.20). Then T is compact and we have representation (2.12). 




3. Strictly Proportional Damping 

It is known [1] that to obtain a “strictly proportional” damping model where the 
damping is strictly proportional to the mode frequency, the abstract wave-equation (2.2) 
must have an additional damping operator and be reformulated as: 

Mx(t) + Ax(r) + Dx(t) + Bu c (t) = 0 (3.1) 

where, as in our case, M commutes with A, we must have 

D = 2^<M<A (3.2) 

where £ is the damping constant, 0 < £ < 1 . The eigenvectors of A, defined by 

M >* = 


and 

4> t (s) = a k sin co*Vpa/G/ v (£ + s) , -t < 5 ^ £ 


where co t satisfies 


and as a result 


Rewriting (3.1) as: 



Y(t) = ^r(0 + Su c (t) 


(3.3) 


(3.4) 


where 


A 


0 I 

-M~ l A -AT'd 




C 


0 

-M~ l Bu c 


where 


BO 


x\ e 3)(-Ja) 

x 2 e X 

and we introduce the energy inner product on 3)(^JA) x X by 

|| y || 2 = [^Ax,,y/Ax 2 ] + [Mx 2 ,x 2 ] 
the eigenvectors of A are given by 



‘ 4> 

. YM * J ’ 



and the eigenvalues are 


K - -£©* ± - C 2 


where 


y = -£ + Wl - £ 2 . 


We have thus strictly proportional damping in the sense that 


R e- K 

Im. X k 

The concrete version of (3.1) is 


= constant 




pa £ML-H _ + [♦(,)!„(, ,t) 

' l 2 1 A 

R(s, a) da + J M(s, a) ^ u(t , a) <io 


(3.5) 


-t 


= 0 


GIyu\t, l) + I 4 U(t, Z) + 2^V/ 4 G/ v 


t -o 

J *»co afaj ^ + u cW = ° 

-t 


(3.6) 


where 


R(s.o) - ^ cot ’t' ^ - £- g + AT(j + o). 


6/ 



The main unnatural feature of this operator is the presence of the new terms in the 
boundary equation (3.6). Unfortunately, this is essential for strictly proportional damping. 
It would be interesting to search for an operator without this feature that yet retains 


asymptotically proportional damping: 


Re. X k 


Im. X k 

vi - 


In particular this would also retain the main feature of (3.4) in that A generates an 
analytic semigroup [1]. Of course the eigenfunctions <>*(•) approach those of a,, as k -> «, 
because of (3.3) — <!>*(*) goes to zero as * 
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Abstract 


A special class of nonlinear damping models is investigated in which the damp- 
ing force is proportional to the product of positive integer or the fractional power 
of the absolute values of displacement and velocity. For a one degree of freedom 
system, the classical Krylov-Bogoliubov “averaging” method is used, whereas 
for a distributed system, both an ad hoc perturbation technique and the finite 
difference method are employed to study the effects of nonlinear damping. The 
results are compared with linear viscous damping models. The amplitude decre- 
ment of free vibration for a single mode system with nonlinear models depends 
not only on the damping ratio, but also on the initial amplitude, the time to 
measure the response, the frequency of the system, and the powers of displace- 
ment and velocity. For the distributed system, the action of nonlinear damping 
is found to reduce the energy of the system and to pass energy to lower modes. 

1.0 Introduction 

One of the major challenges remaining in the development of large space struc- 
tures is to determine a damping mechanism in order to stabilize flexible flight 
structures such as solar arrays, antennas and platforms. As the size and flexi- 
bility of space structures increase, the need to characterize energy dissipation in 
a more appropriate and accurate manner also increases. Under the assumption 
of linear viscous damping, the amplitude decrement of free vibration depends 
only on the damping ratio, regardless of what the frequency or initial condi- 
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tions might be. Numerous experimental results, such as those in the Spacecraft 
Control Laboratory Experiment (SCOLE)^ 1,2 ^, indicate that this is far from suf- 
ficient and that there is a great need for understanding the damping mechanism 
which may be inherently nonlinear. 

Various nonlinear models, such as linear dampers with clearance, Coulomb 
friction dampers, velocity-nth power damping, etc., have been investigated in the 
past In many cases, these models can be represented by a damping force 

that is proportional to the product of integer or fractional powers of the absolute 
values of displacement and velocity, Balakrishnan introduced this nonlinear 
model in i 6 l and obtained approximate solutions using the Krylov-Bogoliubov 
“averaging” method i 7 h He also showed that these results can be quite useful to 
study the response of flexible structures to nonlinear boundary feedback control. 
In this paper we further study this special class of nonlinear damping models. 
We use the Krylov-Bogoliubov “averaging” technique for a one degree of freedom 
system and employ both an ad hoc perturbation method and a finite difference 
technique for a distributed system. 

This paper is organized as follows. In Section 2, the approximate equations 
of amplitude are derived for a single degree of freedom system with nonlinear 
damping. In Section 3, the transient response of free vibration of a single mode 
nonlinear system is compared to that of a system with linear viscous damp- 
ing. In Section 4, the perturbation solution is derived for the vibration of a 
pinned-pinned beam with nonlinear damping. In Section 5, the vibration of the 



pinned-pinned beam with nonlinear damping is simulated via finite-difference 
methods and the results are compared with those obtained using the perturba- 
tion solution discussed in Section 4. 

2.0 Single Degree of Freedom System 

The classical Krylov-Bogoliubov “averaging” method, introduced in 1947, is 
basically a method of variation of parameters. Over the decades, this averaging 
technique has been employed to study nonlinear mechanics and solutions can 
be found in the literature [8] for special cases of nonlinear differential equations. 
Balakrishnan applied this averaging method to a particular class of nonlinear 
damping models [6] which will be discussed in detail in this section. Since this 
damping model is representative of a variety of nonlinear damping mechanisms, 
we further study the effects of the model on vibrating structures of the special 
class of nonlinear damping represented as 


m'x + c\x\ a \x\ b x + kx = 0 

(la) 

x + 7 |x| a |i| 6 i + w 2 x = 0 

(16) 

x(0) ~ Aq } and i(0) = 0 

(lc) 


where -y = — and w 2 = — , c is damping constant, and “a, b”> 0 

• m m 

Note that the term "f\x\ a | i J 6 x represents the dissipating effect of a nonlinear 
damper with “a” and “b” both being positive integers or fractions. 



When 7 << 1 (c small relative to m) we may apply the averaging method 


of Krylov-Bogoliubov to obtain an approximate equation 


x(t) = A(t)sin(wt + 


where the amplitude / 1 ( 2 ) and phase angle satisfy the following equations 


dA{t) 

dt 



(2a) 


and 




(26) 


The functions Ko(t) and Pq( 0 in equations (2a) and ( 2 b) are defined as 



r 2?T 

/ D(A sin <j>, Acj cos <fi) cos<j)d<j> 
Jo 

(3«) 

P °M = T*. 

r 

/ D(Asm<f>,Awcos<f>)s\n<j>d<f) 
/ 0 

(36) 


where D(x,i) is equal to \x\ a |i| fc for the choice of nonlinear damper in 
equations (la-lb), that is, D(x,i) = D(Asin<j), Aojcos<f>). Substituting D(x,x) 
into equations (3a) and (3b) we obtain 


K 0 (A) = u; 6 + M a+6+1 p 
Po(A) = 0 

where 


(4a) 

( 46 ) 





| sin <f>\ a \ cos<f)\ b ^ 2 d(f) 



(4c) 



Equation (4b) implies that for the choice of a nonlinear damper represented 
by 7|x| a |z| 6 x, the phase angle <f>(t) does not, on the average, change over time. 

The positive number // in equation (4c) is called the nonlinear damping 
factor. By changing variables in equation (4c), it can be shown that 


o r 7r / 2 

— / | sin <£| G | cos<^| &+2 c/0 

7T J 0 


m 


The above expression for // is very similar to that for the so-called "damping 
force amplitude ratio y n ” proposed by Jacobsen [9] when studying equivalent 
viscous damping. By employing the properties of the Gamma function, it is 
found that 


r[°±i] r[*±3] 

^ r js±i + 2] 


(4c) 


where T(-) represents the Gamma function. Furthermore, when “a” is an 
odd number , it can be shown that 


/i = 


2 ( n T 1 ) T2 ! 


(6 + 3)(4 + 5)-*-(H-3 + 2n) 


for any b > 0 


With the above information we are in a position to derive equations for the 
amplitude A(t) and displacement x(t) for a system with nonlinear damping. 

For the nonlinear damping system, (“a 4- b” > 0), we substitute equation 
(4a) into equation (2a) to obtain a differential equation for the amplitude ,4(/) 


dA(t) _ 


dt 


= -jfiu b A(t) a+b+l 


7 0 


It can be shown that 


/ m \ 

\cfi(a 4* b)u> b (t + t c )) 


where 


tc = 


cfuv b (a + 6)y4o a+& 


and 


c (0 = 


T^W> 


cos ut 


v cfi(a -f b)uj b (t -f tc), 

Recall that for linear damping viscous damping (“a 4- b” — 0) 


(5a) 


(56) 


(5c) 


x(t) = Aoe ( ‘^ t cosuit 


In equation (5c), the quantity t c is associated with the initial conditions but 
is not the initial time. The constant t c has units of time and is never equal to 
zero. 

Defining t + t c = and replacing A(n^~ — t c ) by ^(n) in equation (5a), 
we obtain an expression for the amplitude in terms of the number of cycles of 
nonlinear vibration 


i(n) 


(b+1) 


K cfi(a + b)k LL s 11 (2nn ) , 




where n is the number of cycles of oscillation. 


(5 d) 


3.0 Transient Response of Nonlinear Damping 

Using the solution for the nonlinear damping system in equations (la) -( 1 c) , 
it can be shown that the logarithmic decrement, 6 , which is the ratio of two 
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successive amplitudes, is given by 


8 = — ^ln f 1 4- 7 - 7 - 7 -') ( 6 °) 

a + b \ t\ ~h tc ) 

In equation (6a), ti is the time when the response is first measured to com- 
pute the amplitude ratio, 6, and T is the period of oscillation. Assuming 1 1 = 0 
and using equation (5b) for t c , we obtain 

6 = —In ( 1 + 2x/ic(ct + 
a + b l m 

For linear viscous damping it is well known that the logarithmic decrement is 

[10] given by 

4 = GtiTTrj) = 2,< (<k) 

Comparing equation (6b) with (6c) one may conclude that while the rate of 
amplitude decay depends only upon the damping ratio, for the linear damping 
model, the amplitude associated with the nonlinear damping system decreases 
more rapidly as: 

a) the initial amplitude Aq increases 

b) the frequency increases 

c) “a + b”, especially “b”, increases 

d) damping ratio £ increases 
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Numerical Examples 


Some numerical results are summarized in Figures 1 through 4, which contain 
the time histories for the nonlinear and linear damping models (a — b ~ 0) and 
for the single degree of freedom system with mass m = 1, stiffness k = 4, and 
damping coefficient c = 0.01. 

For the nonlinear damping cases, “a + b” = 1 in both Figures 1 and 2, 
but “a + b” = 2 in Figure 3, The initial amplitude Aq is 25 in Figures 1 
and 3, and 50 in Figure 2. It can be seen that while the amplitude decrement 
remains constant for ail of the linear damping cases in Figures 1 through 3, it 
decreases more rapidly for the nonlinear damping model. This condition is more 
evident: (a) at the initial time than at a later time, (b) as the initial amplitude 
increases and (c) as “a + b” , especially “b” increases. This is in agreement with 
observations made in practical engineering problems. 

Figure 4 compares the logarithmic amplitude decay in terms of the number 
of cycles of vibration represented by equation (5d) for a system with the same 
parameters as in Figures 1 through 3. The sum of “a 4* b” remains constant 
(“a _{_ b” = 3) as “a ” and “b” vary individually. It is found that the rate of 
amplitude decrease is greater as the value of “b” increases, even though the sum 
of “a” and “b” remains constant. 
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4.0 Perturbation Solutions of Nonlinear 
Damping for Distributed Systems 

Consider the vibration of a pinned-pinned beam with non-linear damping 


pu + C ab \u"\ a \u"\ i u" + EIu iv =0 

( 7a) 

(0, t) = u(L, t) = u"{ 0, t) = u"(L , <) = 0 

(76) 

. mnx 

u(x, 0 ) = A sin — j — 

Li 

(7c) 


where p is the mass per unit length, L is the length of the beam, and C a b is 
the nonlinear damping constant. 

In [6], Balakrishnan applied the Krylov-Bogoliubov method to a multi- 
dimensional system and obtained expressions similar to those in the single mode 
case. A common approach to treat an elastic system, such as the pinned-pinned 
beam discussed here, is to use a modal expansion method to convert a partial 
differential equation into a series of ordinary differential equations; however, 
great difficulties were encountered in the modal expansion due to the presence 
of the absolute value function in equation (7a). For this reason, both an ad hoc 
perturbation technique and the finite difference method presented in Sections 4 
and 5 were utilized to study the effects of nonlinear damping for a pinned-pinned 
beam. 

For small values of nonlinear damping coefficient, C a b , the system will oscil- 
late at the frequency 


74 


and period 


/m7r\ 2 (El) 


V L J 


(Id) 


Pm =2n/w r] 


Let the perturbation solution after one period be 


(7e) 


u(x,p m ) = u°(x,p m ) -f At*(x,p m ) 

where u°(x,t) is the unperturbed solution of equations (7a-7c) with C Q b = 0, 
and A u(x,t) the perturbed solution. In other words 


u°(x, pm) — A sin • 


and 


kirx 


(8a) 


( 86 ) 


A u(x,p m ) = ^2 Aksin - 

k = l 

The values of the coefficients Ak in equation (8b) indicate the degree to 
which the initial mode is damped out and the other modes are excited. The 
perturbed solution in equation (8b) is given approximately by 


Ati= — \u"\ a \u"\ b u" 

P 


(9a) 


with u(x,<) on the right hand side of equation (9a) being further approxi- 
mated by tx°(x,2), that is, 




u(x,t) = u°(x } t) = A sin— —cos u m t (96) 

Li 

hence 

4 ( 771 7T \ 2 . mitX . / 0 n 

/(x, f) = -A (*77; 5m— Y^ C0SUmt 
Substitution of equations (9b-9c) into equation (9a) gives 

Ati = 0 m \sin^^\ a+b sin^^\cos u m t\ a • \sinuj m t\ b sinw m t (9d) 

L L 

where f3 m is a constant that depends on many of the system parameters, such as 
the nonlinear damping coefficient, C a t, the amplitude A, and the mode number 
m, etc., that is, 

Q A Q ( }\ A a + h + 1, ,fe+l (™Y f 

Pm — Pm\ a yV) — A y £ ) 

where / = a + 6 - 1-1 

The coefficient p m is proportional to 2(a + 6 + 1) powers of the mode number 
m , hence, Ati in equation (9d) may become quite large for higher modes. This 
suggests that the perturbation solution may apply only to low frequency modes 
and not to higher frequency modes. 

Because we are mainly interested in the maximum perturbation of the ampli- 
tude decay due to nonlinear damping after half a period, we integrate equation 
(9d) twice from 0 to t(t < p m / 2) to obtain 
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* / Prn \ a , • m * x 

Au(x, — ) - p m CT\sin — 


a + 6 « 


(9c) 


where 


Ct 


g r'\tf 

J O Jo 


cosu/ m J| a | sin u> m /| 6 sina; m rc/r 


The quantity ct is a constant that depends on the half period p m as well 
as “a” and “b”, but does not depend on x. In order to examine the interaction 
between different modes in Au(x,^) in equation (9e), we expand 
into a sine Fourier series in terms of the multiples of the initial mode sin 
that is 


oo 

\ / P m \ A • nm7ra? /An 

Au(x, — ) = 2_s A nm sin — - — (9 / ) 

n-1 

It can be shown that if the pinned-pinned beam in equation (7) is excited 
only by a single mode A sin , then for any positive numbers “a” and u b”, 
all of the even multiples of the initial mode will not be excited; only the odd 
multiples of the initial mode get excited. In other words 


Anm = 0, for n- 2 j 


Anr, 


2/? m Cr 


/'I- 


m 7T3C 

L 


I a+^sin mpsin ™s*dx 


| for n = 2 j - 1 

Now expand Au(x, Ej 5 -) into a sine Fourier series to obtain the Fourier coef- 
ficients A nm , as represented in equation (9f). 
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It can be shown that for a = 0, 6=1 


( 10 ) 


Anm = o, for n = 2 j 

. ^ = [ (2j-3K2^lK2;>l)]p^ ' f° r U = 2j ~ 1 
Similar expressions can be derived for any positive integers “a” and u b” . 
Equation (10) relates the initial mode A sin to the distributed nonlinear 
damping coefficient A nm . Because the perturbation solution contains compo- 
nents of other modes, they too will be excited. The degree of excitation is 
illustrated by the values of A nm . It is obvious that the even-multiples of ini- 
tial modes are not excited while all of the odd-multiples are. It will be shown 
later that for a beam excited by fundamental modes, equation (10) provides an 
estimate for the amplitude decay after a half period due to nonlinear damping. 


5.0 Finite Difference Simulation for Nonlinear Damping 

In this section, the transversal vibration of a pinned-pinned beam excited by 
a single mode is simulated via the finite-difference method. We will present 
a finite difference scheme for solving the partial differential equation, derive 
the stability conditions and then discuss numerical results. Specifically, we 
will compute amplitude decay after half a period and compare the results with 
both linear and nonlinear damping. We treat only the case of a = 0, b = 1 
for nonlinear damping. Finally, we compare the amplitude decrement using 
the finite difference method with the amplitude decrement obtained using the 
perturbation method presented in Section 4. 



Finite Difference Formulations 


When modeling the nonlinear differential equation (7a) with the finite dif- 
ference method, we represent beam displacement as u(xi,tj) = where Xi 
is defined as Xi = iAx , and tj is defined as tj == jAt. Both grids Ax and At 
are defined as Ax ^ —.and At = — , where m and n are the numbers of grids 
used for L and T. L is the beam length and T represents the half period for the 
initial mode. With the above definition, the nonlinear finite difference equation 
(corresponding to equation (7a)) for the displacement at tj + \ is given by 


Ui,j + 1 = 7( u ;+2,j - 4« i+ ij + 6 Uij - 4u,_i i j + + (2 M.-J - i) 

l) 2(u, J + ( w t— l p j u i— lj — l)l 

■ [( + 1 ,j — l) 2(t/<j t/ij' — i) + w t— i,j — i)] 

(11a) 

If linear damping is assumed (“a” = “b” = 0), then equation (7a) becomes 

pii + cu" + EIu iv =0 (116) 

and the corresonding finite difference equation is 

1 = 7(u< +2 j - 4u i+ ij 4- 6ui i; - - 4u,_ij + «.-2,j) + (2m., ; - u»,;-i) 
+Q-[(Mi + l i j — i) 2 (m,j M.,j_x) 4" (m.— 1 , j M,_i ( j_i)] 

(lie) 
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Finally, if no damping is present, equation (11c) is further simplified as 


t/jj+i = 7(u i+2 ,j - 4u,-+i j + 6 Ui t j — 4uj_ij + + (2«»j 


where 


(lid) 


El At 2 


C ab 


cAi 


7 = - 7 -, /? = ~z~~a > ar, 4 « = -- T— r 

? pAa; 4 pAx 4 Ax 2 


(lie) 


Formulas for central differences have been used to approximate both u 
and ii. Forward differences for t and central differences for x are used to obtain 
the mixed derivative ii" in equations (11a), (11c), and (lid). Two fictitious 
boundary conditions have been created such that the zero moment conditions 
are satisfied at both ends for the pinned-pinned beam. They are 


ttm+lj — an< l U -1 ,j — U 1 J 

j = 0, 1,2, ...n 

Stability Conditions 

It should be noted that equations (11a), (11c), and (lid) belong to the 
explicit forms of finite difference formulation, for which there always exists a 
stability problem for the specific finite difference scheme [11]. In other words, 
when the finite difference scheme is stable, there exists an upper limit to the 


80 


extent to which any error arising during the simulation can be amplified. This 
implies that the numerical solution will not diverge. Clearly, stability alone 
does not necessarily guarantee that the deviation between the true solution to 
a certain partial differential equation and its finite difference approximation 
will be small in any sense. Stability only implies the boundedness of the finite 
difference solution, at a given time, as At approaches zero. 

In the case of a beam without damping, the equation of motion is 

pii + Elu iv = 0 ( 1 2a) 

Assuming harmonic motion, that is, 

u(x t t) = e* Xx T(t) (12 b) 

Using the finite difference method, equation (12a) can be replaced by 
T(t + At) - 2 T(t) + T(t - At) + i[(e 2jAx + c~ 2 ^ T ) 

-4 + = 0 (12c) 

Since [e 2j — 4(e** x + e” J Aar ) -f 6 ] = 16cos 4 (4p), equation ( 12 c) 

is simplified to obtain 


T(t + At) - [2 - 16 7 cos 4 (— )]T(t) + T(t - At) = 0 (12rf) 

Substituting 7 in equation (lie) into equation ( 12 d) and utilizing the stabil- 
ity criteria for the eigenvalues of a difference equation, we find the stability 
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condition for the finite difference equation (lid) to be 


EIAtr 
p Ax 4 


cos4 (4r) < 


1 

4 


(12c) 


or 


El At 2 1 


( 12 /) 


pAx 4 — 4 

For a beam vibrating with linear damping, represented by equation (lib), 
the difference equation for the stability condition is 


A y, /\ Jp 

T(t + At) - [2 - 167cos 4 ( — ) - 4asin 2 ( — )]T(<) 


+[1 — 4a sin 2 (-^j- )]T(t — At) = 0 


(13a) 


This can be further simplified to 


AEIAt 2 

pAx 4 


4 ,& x \ , 

cos (— ) + 


2(^y/ E Ip At . 2 /Ax . 


or 

El At- CVEI pAt ^ 1^ (136) 

pAx 4 2pAx 2 “ 4 

which is similar to equation (12f). 

Notice that the stability condition similar to equations (12f) and (13b) is 
not available for the nonlinear difference equation. This is mainly because of 
the difficuties encountered due to the presence of the absolute value function in 
equation (11a). 
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Numerical Examples 


The numerical example used for the simulation is a pinned-pinned beam which 
resembles the SCOLE project mast [1] with L — 130 ft, El — 4E07 lb. ft 2 , 
p — 0.09556 slug/ft. The beam is initially at rest and is excited by a single 
mode /fsin 1 ^ with A— 1.3 ft and i = 1,2, 3, 4. The first frequency is 11.95 
rad/sec and its corresponding period is 0.5258 second. We proceed to evaluate 
the displacement of the beam for the first half period T. 

Because the stability condition ( 1 2f ) implies that Ax can not be arbitrarily 
small, we first chose m = 10 (Ax* = L/m = 13/2) when verifying the stability 
condition. It was found that if n > 64, At = T/n < 0.0041, |y| = 0.2473 < 0.25, 
then the numerical scheme was always stable. On the other hand, if n ~ 63, \y\ — 
0.2552 > 0.25, the numerical scheme was found to be unstable, which confirms 
our stability criteria represented by equation (12f). 

Before we compare the amplitude decrement using the perturbation and the 
finite difference methods, it is necessary to verify the finite difference scheme 
used in the simulation. We compare the finite difference displacement results 
with those obtained using an analytical solution after a half period for a beam 
vibrating without damping (corresponding to equation (12a)). We also com- 
pare the finite difference results with that for the logarithmic decrements of the 
amplitude for a beam vibrating with linear viscous damping (corresponding to 
equation (lib)). We believe that as long as the finite-difference scheme con- 
verges reasonably well for both equations (12a) and (lib), the scheme should 
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work for equation (11a) with nonlinear damping. 

The relative errors were computed for the displacement of a beam vibrating 
without damping after half a period when 20 intervals in L and 1100 intervals in 
T were used. The relative errors are defined as the difference between the dis- 
placement using finite difference method and that using an analytical method. 
The difference is then normalized by the analytical values. The relative errors 
are quite uniform for each mode over all locations of the beam with their maxi- 
mum values being 0.0025%, 0.037%, 0.23%, and 0.54% for modes 1, 2, 3, and 4 
respectively. This implies that the finite difference method provides quite accu- 
rate results for beam vibrations without damping. Notice that the errors can be 
greatly reduced if more intervals are used for L and T, subject to the stability 
condition of equation (12f). 

For free vibration of a pinned-pinned beam with linear viscous damping, 
represented by equation (lib), it can be shown that by assuming harmonic 
motion (as in equation (12b)), the damping ratio, for each mode depends 
only on damping constant c, beam properties p and El , but not on the mode 
number i, that is, 


2 


(14a) 


According to classical vibration theory [10], it is well known that the ampli- 
tude decay after half a period for single degree of freedom system is 
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( 146 ) 


A Ayl 

< = — = * 


which also does not depend on the mode number i. Specifically, if ( = 0.003, 
then theoretically S should equal to 0.0094 after half a period for any single 
DOF system. Using the same numbers of grid points for L and T as before, 
the amplitude decay, 6, after half a period for modes 1 and 2 at a variety of 
locations on a beam with linear damping (£ = 0.003) were computed. For both 
mode 1 and 2, the amplitude decay, 6, remains almost identical at all locations 
of the beam, as if each point on the beam vibrated as a single DOF system. It 
is interesting to note that the amplitude decay, 5, is about 0.0094 for mode 1, 
and 0.0096 - 0.0097 for mode 2, both of which are close to 0.0094 (— 7 t£). These 
results are in good agreement with equations (14a- 14b) and further confirm that 
the finite-difference algorithm used is quite reliable. 

Figures 5 and 6 compare the values of the amplitude decay for modes 1- 
3 after half a period T for the same beam with: (1) linear damping using 
finite-difference method, (2) nonlinear damping using finite-difference method, 
and (3) the perturbation method discussed in Section 4. The abscissas are a 
dimensionless quantity damping ratio £. Linear damping coefficient c is equal 
to 2QyJpEI y and the nonlinear damping coefficient C a b si equal to 4c. These 
values are used for both finite-difference and perturbation methods. For the 
case of linear damping, the amplitude decay takes almost the same value for 
modes 1-3 (e.g., 6 = 0.0064 for modes 1-3 when ( — 0.02), which is consistent 
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with equation (14a). 

From Figures 5 and 6 we observed that: (1) the effect of linear damping 
exceeds that of nonlinear damping for mode 1, (2) the effect of linear damping 
is about the same as that of nonlinear damping for mode 2, and (3) the effect of 
nonlinear damping outweighs that of linear damping for mode 3. These results 
imply that in the case of nonlinear damping the amplitude of higher modes 
will be damped out faster than the lower modes, whereas in the case of linear 
damping, the modes are all equally damped. This is consistent with the results 
of Section 3 for a single degree of freedom system, as well as with common 
engineering judgement. 

Finally, when comparing the results in Figures 5 and 6 for nonlinear damping 
using either finite-difference or perturbation methods, it is observed that for 
fundamental modes the perturbation method provides a very conservative upper 
bound of amplitude decrement after half a period. This result might be useful in 
preliminary assessment of the impact of nonlinear damping effect, since in many 
occasions, we are mainly concerned with the system’s fundamental frequencies. 

Conclusions 

A space systems dynamics and controls analyst is often confronted with the 
problem of gaining a better understanding of the damping mechanism which is 
inherently nonlinear. Fortunately, some of the difficulty in handling nonlinear- 
ities is offset by the fact that damping is still small. This makes it possible to 
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obtain approximate solutions using the classical Krylov-Bogoliubov “averaging” 
technique to study a class of nonlinear damping models. 

In this paper, the damping force is assumed to be proportional to the product 
of positive integer or fractional powers of the absolute values of displacement 
and velocity. As is expected that for a typical nonlinear system, the amplitude 
decrement of free vibration with nonlinear models depends not only on damping 
ratio, but also on the initial amplitude, the time to measure the response, the 
frequency of the system, and the powers of displacement and velocity. For a 
pinned-pinned beam, both an ad hoc perturbation method and a finite difference 
technique are used to study the vibration of a beam with nonlinear damping. 
The action of nonlinear damping is found to reduce the energy of the system 
as well as to pass energy to lower modes. As a result, the amplitude of higher 
modes will be damped out faster than the lower modes. All of these results 
are very useful to study the response of a flexible structure to the action of a 
nonlinear boundary feedback control. 
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Figure 1 Time History for Linear and Nonlinear Damping: a+b — 
A n = 25, Mass = 1, Stiffness = 4, Damping Corf. = 0.01 

Figure 2 Time History for Linear and Nonlinear Damping: a+b = 
A 0 - 50, Mass = 1, Stiffness = 4, Damping Coef. = 0.01 

Figure 3 Time History for Linear and Nonlinear Damping: a+l>=2 
A 0 = 25, Mass = 1, Stiffness = 4, Damping Coef. = 0.01 

Figure 4 Amplitude vs. Number of Cycles for Linear and 
Nonlinear Damping 

Figure 5 Amplitude Decay vs. Damping Ratio after Half Period 
FDS: Finite Difference Simulation 

Figure 6 Amplitude Decay vs. Damping Ratio after Half Period 
FDS: Finite DifTcrence Simulation 
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ABSTRACT 

Large spacecraft such as Space Station Freedom employ large trusses in 
their construction. The structural dynamics of such trusses often exhibit 
nonlinear behavior and little damping which can impact significantly the 
performance of control systems. The mini-MAST truss was constructed to 
research such structural dynamics and control systems. The mini-MAST 
truss is an object of study for the Guest Investigator Program as part of 
NASA's Controls-Structures Interaction Program. The Mini-MAST truss is 
deployable and about 65 feet long. Although the bending characteristics of 
the Mini-MAST truss are essentially linear, the angular deflection under 
torsional loading has exhibited significant hysteresis and nonlinear stiffness. 
It is the purpose of this study to develop nonlinear and distributed 
parameter models of the truss and to compare the model dynamics with 
actual measurements. Distributed parameter models have the advantage of 
requiring fewer model parameters. A tangent function is used to describe 
the nonlinear stiffness in torsion, partly because of the convenience of its 
easily expressed inverse. Hysteretic slip elements are introduced and 
extended to a continuim to account for the observed hysteresis in torsion. 
The contribution of slipping to the structural damping is analyzed and found 
to be strongly dependent on the applied loads. Because of the many factors 
which affect the damping and stiffness in a truss, it is risky to assume 
linearity. 


INTRODUCTION 

Future missions in space require spacecraft which are considerably 
larger and more flexible than current spacecraft. Large spacecraft such as 
Space Station Freedom employ large, complex trusses in their construction. 
The structural dynamics of such trusses often exhibit nonlinear behavior 
and low structural damping which can impact significantly the performance 
of control systems. For example, in reference 1, Lallman studies the effect 
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of damping on the performance of the attitude control system of the Space 
Station Freedom. The mini-MAST truss was constructed to research the 
interaction of such structural dynamics and control systems and is an object 
of study for the Guest Investigator Program as part of NASA’s Controls- 
Structures Interaction Program. 

The Mini-MAST truss was designed to be deployable to a length of 
66.14 feet when fully extended. The bending characteristics of the Mini- 
MAST truss are essentially linear. The angular deflection under torsional 
loading, however, has exhibited significant hysteresis and nonlinear stiffness 
during laboratory tests. 

The complexity of such structures create a burden to optimal design 
and to systems identification for upgrading dynamic model parameters by 
analyzing experimental test data. The large number of model parameters 
which results if each structural mode is assumed to be independent can be 
greatly reduced if distributed parameter models are used. 

It is the purpose of this study to develop distributed parameter 
models of the Mini-MAST truss and to compare the model dynamics with 
the actual dyilamic characteristics. A second purpose is to model the 
nonlinear stiffness and damping properties of this joint-dominated truss. It 
is hoped that the study results will be useful in designing control systems 
for large spacecraft such as Space Station Freedom which employ similar 
trusses. 


DISCUSSION 

Because the Mini-MAST truss is representative of structures that will 
be used for large spacecraft such as the Space Station Freedom, the study of 
its structural dynamics is valuable in assuring the dependability and high 
performance of spacecraft control systems. Figure la. pictures the Mini- 
MAST truss being deployed. The reduction in volume is striking when 
compared to the deployed truss shown in figure lb. Reference 2 describes 
in detail the design of the Mini-MAST. Because of the complexity of the 
truss it is important to study simplifying models of its dynamics. Figure 2 
shows how many modes are required to depict accurately the static 
deflection of a cantilevered beam. The problem is compounded if the modal 
parameters are considered to be independent. Because of the resulting 
complexity there is considerable advantage in using distributed parameter 
models. Due to the greatly reduced numbers of parameters required for 
such models as shown in figure 3, the ability to employ systems 
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identification (Reference 3) and optimal design techniques is greatly 
facilitated. Because of these advantages it is valuable to determine the 
accuracy with which distributed parameter models can represent the Mini- 
MAST truss. For example, can such simple models predict accurately the 
peaks of the frequency response shown in figure 4? If distributed 
parameter models represent accurately the dynamics of the Mini-MAST 
truss, then the model equations can be used to upgrade the model 
parameters using systems identification. Also the models will be useful in 
integrated control-structures design because their form provides easy access 
to global varibles such as the modulus of elasiticity. 

The Mini-MAST truss, being deployable, requires a large number of 
joints. The compliance and possible slippage of the joints may affect the 
overall stiffness of the truss when viewed as an equivalent beam. The 
action of the joints may also affect the damping of the truss as well. It is 
important to know accurately the damping of a spacecraft in order to assure 
reliable and high performance control. It is also important to understand 
and to model any nonlinear behavior caused by the numerous joints. 


Distributed Parameter Bending Model 

The Mini-MAST truss is modeled as a cantilevered beam with an 
added tip mass as depicted in the schmetic in figure 5. The partial 
differential equations (Euler beam equation) and boundary condition 
equations (Cantilevered and tip mass) are solved thereby determining the 
modal characteristics. First, the calculated static deformation resulting from 
a constant 15 pound for applied to the tip is compared to actual test results 
in figure 6. The value of the stiffness parameter, El, for an equivalent Euler 
beam derived from this test is 27.6 x 10 6 pound feet squared. The 
comparison suggests that the model deformation matches the actual 
deformation within the measurement error. The resulting modal 
frequencies in bending are then compared with experimental results and 
those for a finite element model (Reference 4) in figure 7. The frequencies 
for the first few bending modes of the distributed parameter model 
accurately match the actual bending frequencies of the truss. At higher 
mode numbers, however, the actual modal frequencies are lower than the 
theoretical values for the Euler beam model. Belvin, in reference 5, showed 
that the shear deformation of a similar truss cannot be ignored as is done in 
the Euler beam model. Belvin used the techniques of reference 6 in his 
study. The Timoshenko beam, in contrast, accounts for the shear 
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deformation and more accurately models the frequencies in bending as 
shown in figure 7. 

Figure 7 also shows the accuracy with which the frequencies of a finite 
element model match the actual frequencies of the truss. The finite element 
model is reasonably accurate even at high mode numbers. The parameter, 
El, used in the finite element model equals 29.8 x 10 6 pound feet squared. 

In figure 8 the bending mode shapes generated by the same finite element 
model exhibit shapes similar to Euler beams with one exception. Examina- 
tion of the third mode reveals that the shear deformation is significant 
enough to give a change in slope almost at the bottom of the truss. The 
general contour of the mode shapes in figure 8 compare well with those of 
the Timoshenko beam (not shown) but the irregularities which show 
significant local deformation will be missing from the distributed parameter 
models. It is possible that overlookong such local deformations could cause 
control system instability. 

The effect on the first bending mode frequency of changing the mass 
at the tip of the equivalent beam is shown in figure 9. The frequency 
response measurements of figure 4 had a tip mass which weighs 70.125 
pounds (mass ratio = .31). The Mini-MAST truss excluding its tip mass 
weighs 229 pounds. The Euler beam model depicts accurately the change in 
frequency when the tip mass is removed. The assembly for the active 
control of the Mini-MAST is expected to weigh in excess of 300 pounds. The 
frequencies for higher mode numbers will not change as much as that for 
the first mode because as mode number increases the motion of the tip mass 
diminishes, thereby approaching a pinned end condition. 


Distributed Parameter Torsion Model 

Similar to the bending case, the truss is modeled in torsion as an 
uniform shaft which is fixed at one end and has a tip body attached to the 
other end. Based on the angular deformation due to an applied moment the 
torsional parameter GI po iar equals 2.16 x 10 6 pound feet squared per 
radian. The partial differential equations and end conditions are solved and 
in figure 10 the model's torsional frequencies are compared with 
experimental results and the finite element model of reference 4. The close 
comparison indicates that the modal frequencies for both the distributed 
parameter model and the finite element model compare closely with the 
actual frequencies. 
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Nonlinear Torsional Stiffness 


Because the Mini-MAST truss exhibits significant nonlinear stiffness 
and hysteretic behavior in torsion, it is necessary to model these 
characteristics. The nonlinear stiffness model will be discussed first. 1 he 
hysteretic model will be treated in the next section. 

Although the form of the nonlinear stiffness is approximately cubic, 
a tangent function is used because (1) its form gives the nearly linear plus 
cubic relationship that is needed, and (2) the tangent has a conveniently 
express inverse. Figure 11 depicts the tangent model of the nonlinear 
stiffness in torsion and introduces the parameters, K, and, B, which govern 
the linear and the cubic contribution, respectively. The parameter, K, then is 
the usual torsional stiffness. 

In figure 12 it is evident that the tangent relationship compares well 
with the experimental results. The data shown is believed to not involve 
any slipping as it represents the relaxation from a load having been applied. 
As the load is increased slipping does take place and will next be considered. 


Torsional Slip Model 

The torsional hysteretic model is comprised of an infinite number of 
slip elements. An individual slip element is assumed to slip instantaneously 
upon reaching a particular moment threshold. A reverse slip is assumed to 
take place at a moment of equal level but opposite sign as depicted in figure 
13. A slip distribution function is introduced which describes the 
probability density function of the values of moment threshold. The second 
order exponential form of the function, shown in figure 14, was chosen to fit 
the experimental data. Effort is underway to link this distribution function 
to the vertical loading of the joints. The total deflection amplitude consists 
of (1) the deflection due to compliance without slipping plus (2) the 
deflection due to an accumulation of slips due to the applied moment. The 
expected value of the accumulation of slips is given by the integral of the 
slip distribution function between the last moment reversal or zero and the 
current applied moment. The deflection equation is depicted in figure 15. 

The total hysteretic model which contains both the nonlinear stiffness 
and the hysteretic slipping is compared with actual test results in figure 16. 
The close comparison of the model results and the actual hysteretic behavior 
gives validity to the model for torsional deflection due to applied moment. 
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The hysteretic behavior is expected to be dependent on the vertical loading. 
When the 300 pound plus active control assembly is attached to the top of 
the Mini-MAST the total angular deflections are not expected to change 
significantly, but an increase in the moment threshold is expected. Because 
of the effect of gravity it is difficult to determine the hysteretic behavior in 
an unloaded condition as in space. 


Structural Damping 

The damping for the first bending mode is affected by the mass of the 
tip body as shown in figure 17. The damping ratio which was measured for 
the truss without tip mass was about 3.3%. This value was about three time 
the value expected based on the assumption that the dimensional damping 
of the truss would not change. The damping ratio would be expected to 
double from the value of about .45% for the 70 pound tip mass. This 
discrepency is probably due to slipping being affected by vertical loading, as 
is the case for torsion. 

In torsion it is possible to link slipping to damping by accounting for 
the loss of energy due to slipping. Figure 18 shows that the expected 
contribution to damping from slipping for oscillations about the unloaded 
condition reflect the shape of the slip distribution function. The damping 
contribution for oscillations in torsion about a loaded condition may be as 
low as zero because of the complete lack of slipping. 

The statically determinant truss to be used on the Space Station 
Freedom can be expected to involve internal loading. As a result the 
damping of the truss for small amplitudes is not expected to involve slipping 
and will consequently exhibit very low damping. 

Laboratory tests have revealed a damping ratio for bending modes for the 
cantilevered truss to be about .0045. The damping ratio will decrease when 
large bodies are added to the truss. In the absence of air, the damping can 
be expected to be even smaller, perhaps approaching .002. 

Past practices of using a constant damping ratio of .005 for space station 
studies does not represent the worst case. Lower values of damping should 
be used which reflect mass loading and internal loading effects. 



CONCLUDING REMARKS 


The Mini-MAST truss has been tested and analyzed for the purpose of 
understanding the dynamic characteristics, nonlinear stiffness and 
hysteretic damping of large spacecraft. 

It was necessary to use a Timoshenko beam model for bending to 
account for the shear deformation of the Mini-MAST truss. The modal 
frequencies of the Euler beam model were higher than the actual values. 

A tangent function model of the nonlinear torsional stiffness was 
developed and its parameters estimated to match experimental results. 

A hysteretic slip model for torsion was developed using the experimental 
test data. The slip distribution function used has a second order, 
exponential form. The hysteretic behavior is expected to be affected by 
changes in the vertical loading due to gravity. 

The damping contribution in torsion of the hysteretic behavior was 
deduced by analyzing the torsional slip model. The damping due to slipping 
was determined to be quite dependent on loading conditions. A steady load, 
for example, might eliminate slipping and consequently any damping 
contribution due to slipping. 

Future studies of control system performance should use lower values of 
structural damping than the .005 used in the past, and should consider the 
nonlinear effects. 
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Figure 3. Comparison of the Number of Model Parameters for Modal Models 
and Distributed Parameter Models. 
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Figure 4b. Frequency Response of the Mini-MAST Truss in Term of Inches 
of Response per Pound of Input. Frequency from 20 to 40 Hz. 


112 


Log Amplitude 


t 1-1 ! L A l I 1 : 1 * 1 * 1 1 1 k L 

“t • j 

-i I 1 1 



40 50 60 

Frequency - Hz 


Figure 4c. Frequency Response of the Mini-MAST Truss in Term of Inches of 
Response per Pound of Input. Frequency from 40 to 60 Hz. 
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Figure 4d. Frequency Response of the Mini-MAST Truss in Term of Inches of 
Response per Pound of Input. Frequency from 60 to 80 Hz. 
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Figure 5. Schematics of Distributed Parameter Models for Bending and 
Torsion. 
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Figure 6. Comparison of the Model and Actual Static Deflection in Bending of 
the Mini-MAST Truss Subjected to a 15 Pound Force at the Top. 
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Figure 7. Comparison of Distributed Parameter Model, Finite Element Model 
and Actual Normalized Bending Frequencies. 
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Figure 8. 


Finite Element Model Mode Shapes for Bending. 
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Figure 9. Comparison of the Model and Actual First -Mode Bending 

Frequencies as a Function of Tip Mass to Truss Mass Ratio. 
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Figure 10. Comparison of Model and Actual Torsion Mode Frequencies. 
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Figure 11. Nonlinear Stiffness Formulation Used for the Torsion Model. 
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Figure 12. Comparison of Model and Actual Static Deflection in Torsion. 
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Figure 13. An Individual Nonlinear Slip Element for the Torsion Model. 
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Figure 14. Nonlinear Hysteretic Model for Torsion. 
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Figure 15. The Slip Distribution Function. 
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Figure 16. Comparison of Hysteretic Model and Actual Static Defelection in 
Torsion. 
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Figure 17. Comparison of Model and Actual Damping Ratios in Bending as 
Function of Tip Mass to Truss Mass Ratio. 


127 




\ 

Damping 

Ratio 

J 



Figure 18. Effect of the Hysteretic Behavior on Damping in Torsion. 
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EFFICIENCY AND CAPABILITIES OF MULTI -BODY SIMULATIONS 


R.J. VanderVoort 
DYNACS Engineering Co., Inc. 
Clearwater, FI 


ABSTRACT 


Simulation efficiency and capability go hand in hand. The more 
capability you have the lover the efficiency will be. Section 1 of this 
paper discusses efficiency and section 2 deals with capabilities. The lesson 
we have learned about generic simulation is: Don’t rule out any capabilities 
at the beginning but keep each one on a switch so it can be bypassed when 
warranted by a specific application. 


1. EFFICIENCY 


Efficiency means different things to different people. For the person 
running simulations interactively on a terminal quick turn around time is 
efficiency. For the person making 10,000 Monte-Carlo runs low cost is ef- 
ficiency. For the person running real time simulations minimum CPU time is 
efficiency. 

Three aspects of a simulation should be considered when dealing with 
efficiency; hardware, software and modeling. 

Hardware A fast processor will reduce CPU time for a given simulation but 
this doesn’t necessarily equate to improved efficiency. For example, the 
Monte-Carlo simulation may take 10 minutes on a super computer and 2 weeks 
on a PC but if time is free on the PC then that may be an efficient solu- 
tion. We will not discuss hardware related issues except for two points. 1.) 
Fast hardware is of primary importance to the real time simulation because 
it means higher fidelity models can be incorporated 2.) Vector processors 
and parallel processors should use custom algorithms that take full ad- 
vantage of the special machine architecture. 

Software A fast algorithm will also reduce CPU time but again this doesn't 
necessarily equate to improved efficiency. For example, it is generally 
accepted that an ad-hoc simulation is much faster than a generic simulation. 
The cost of developing and testing the ad-hoc simulation may exceed the run 
time saving thereby reducing overall efficiency. 
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Recent work in the area of symbolic programming has shown that sig- 
nificant savings can be achieved by symbolicaly forming t e eqiia on 
motion and numerically solving them. Other aigor i thmshavebeenproposed 
that promise similar savings. There is one point that software 
should keep in mind. With generic simulations the user must have complete 
flexibility in retaining or deleting different parts of his mod . 
because generic simulations are often used for model development and valida- 
tion In^Sat environment an analyst will add or delete certain features to 
determine the effect on performance and whether or not the feature should be 

retained in the model. 

More on this subject in section 2. 

Modeling This is the domain of the simulation user and the area In which 
many improvements in efficiency can be made. For exampie, deieting a high 
order mode in a flexible body model has a compound effect. It reduces the 
model complexity and at the same time allows a bigger l»t'grrtlon step size 
both of which reduce run time. Often times the reduced fidelity is justified 
by the savings in run time. 

The point to be made is that the analyst is the end authority on the 
"correct" P model for a given application. The more flexibility he has in 
S hi" model the easier It Is for him to select the best model for the 

job. 


2. CAPABILITIES 


Capability in our context is synonymous with flexibility and not with 

complexity. A simulation may be very detailed and complex s and has 

changed then it’s only useful in a narrow range of applications and has 

limited capability. 

In our experience with TRBETOPS and DCAP we have found that it is much 
easier to generate a model and obtain a response than it is to predict the 
correct response. In other words, when we don’t get the expected response 
the simulation is usually correct and our expectation is wrong. This is not 
entirely unexpected because it is very difficult, even for an expert, to 
solve the equations of anything but the simplest dynamical systems. The 
solution to this dilemma is flexibility. Start with simple models that have 
known analytic solutions. Then add complexity one step at a time while 
gaining confidence in your model and insight into the behavior of your 

system. 
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For multibody systems vith flexible bodies the same arguments apply but 
the complexity of the model increases more rapidly than for rigid bodies. 
The person doing software development makes assumptions that simplify the 
resulting equations of motion. If this is done carelessly then terms are 
dropped that may prove essential in specific applications. On the other 
hand, if simplifications are not made then the computation burden becomes 
too great. 

The lesson ve learned is that you must retain as many terms as possible 
in the kinematics but they must have associated switches so you can easily 
add or delete them from a specific application. This is done for two 
reasons. 1.) to give you Insight into the effect of various model elements 
on system response and 2.) to allow the selection of the most efficient 
model for a given application. 
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SIMULATION EFFICIENCY 



• BYPASS TERMS • BYPASS TERMS 

• MULTI-RATE ALGORITHMS • INTEGRATION TYPE S STEP SIZE 

• SYMBOLIC PROGRAMMING • REDUCED ORDER 

• ADHOC SIMULATION 



SPEED-UP OPTIONS 
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SOLUTION ALGORITHM OVERVIEW 
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1 . INTRODUCTION 


Supported by the Spacecraft Control Branch of NASA Langley Research Center 
under the Spacecraft Control Laboratory Experiment (SCOLE) program, the Con- 
trol Research Corporation continued the investigation into the control design 
challenges of large space systems and Spacecraft Control Laboratory Experi- 
ment* This study concentrated on the second stage of a two-stage approach to 
active control of the flexible orbital configuration of SCOLE. The principal 
objective was to investigate if the structural vibrations excited by time- 
optimal line-of-sight pointing slew maneuvers of the bang-bang type could be 
quickly suppressed via "modal-dashpot" design of velocity output feedback 
control. 


Structural vibrations in future large space systems such as space anten- 
nas, space platform, space station, or of deployed flexible payloads attached 
to the space Shuttle orbiter, and their interaction with on board controllers 
have become a major concern in the design and operation of such control sys- 
tems as, say, for pointing and stabilization. The natural vibration frequencies 
of such systems are unconventionally low (tenths of 1 Hz in many cases) and 

closely spaced, many of which lie inside or nearby the bandwith of various 

traditional (rigid-body) control systems. In the past few years, many 
approaches were proposed for designing advanced control systems that would 
suppress vibrations In large flexible space structures, and various In— house 
laboratory experiments were also conducted, each being specifically set up for 
demonstrating some particular design techniques of interest. In 1983, the 

Spacecraft Control Branch at NASA Langley Research Center initiated the Space- 
craft Control Laboratory Experiment (SCOLE) program and the NASA/IEEE 

Design Challenge [1] to promote direct comparison and realistic test of dif- 
ferent approaches to control design against a common open-to-public laboratory 
article. As shown in Fig. 1-1, the article was intended to resemble a large 
space antenna attached to the Space Shuttle Orbiter by a long flexible mast, 
similar to the proposed space flight experiments and various space-based 
antenna systems, and to have a truly three-dimensional complex dynamics. 

As stated in Ref. [1], the primary control task of the Experiment is to 
rapidly slew or change the line— of— sight (LOS) of an antenna attached to 
the space Shuttle orbiter, and to settle or damp the structural vibrations 
to the degree required for precision pointing of the antenna. The objective is 
to minimize the time required to slew and settle, until the antenna line-of- 
sight remains within a specified angle. 

Research on a practical two-stage approach and some associated control 
design challenges in the context of SCOLE had been conducted earlier by Lin 
[ 2 ] — [ 5 ] . His initial efforts, also supported by the SCOLE program, were con- 
centrated on "Stage 1" while the flexible-body dynamics of the configuration 
with a flexible mast beam was being actively developed at the Langley Research 
Center. 


Among the most commonly held ideas for pointing/retargeting of large flex- 
ible space systems is the following intuitively appealing and rather practical 
two-stage approach: (Stage 1) slew the whole structure like a rigid body in a 
minimum time under the limited control moments and forces first, and then 
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Fig. 1-1 Spacecraft Control Laboratory Experiment (SCOLE) 
the orbital Shuttle-Maat-Antenna configuration. 




(Stage 2) damp out the excited structural vibrations afterwards. Such an 
approach undoubtedly will be a very relevant, and realistic as well, to study 
with SCOLE.* 

To slew a spacecraft for a given angle in a prespecified ti«e, there are 
many ways to command the slew actuators on board. The one that is easy to 
implement is a bang-bang control. That is, a constant force at its allowable 
maximum is applied in one direction half of the time and then in the opposite 
direction the other half. Such is the most convenient and common with reac- 
tion-jet thrusters, and most spacecraft including the space Shuttle use thrus- 
ters. As the structure considered for future space antennas and optics was 
becoming larger and more flexible, structural dynamicists suggested modifying 
the constant profile by a sine or versine function so as to smooth the switch- 
ing. To explore further in their theorectical and experimental investigations, 
control engineers also started to apply Pontryagin's Maximum Principle of the 
optimal control theory [6] to develop open-loop profiles that would minimize 
excitation of the first few vibration modes.** Including more than a few 
modes generally will make it almost impossible, even with the aid of powerful 
digital computers, to carry out the complicated computations necessary for 
applying the optimal control theory. To implement any such slew profile other 
than the bang— bang type will also require that the thrusters be at least 
"throttleable” in fine steps, which is still beyond the current state of the 
art. 


To slew SCOLE for the desired 20° under the specified limits on control 
moments and forces in a minimum time , instead of some arbitrarily fixed time, 
application of the well-known time-optimal bang-bang control theory [6]-[7] 
was considered the most appropriate for the Stage— 1 design. The theory, how- 
ever, is not directly applicable to SCOLE: due to the asymmetrical configura- 

tion and the moving coordinate frame that is fixed on the Shuttle body axes, 
all axes are tightly coupled through nonzero products of inertia as well as 
through different moments of inertia. After examining the major assump- 
tions in the theory, Lin [2]- [3] was able to develop a useful practical design 
technique for time-minimized single-axis bang-bang slew maneauvers. This in- 
cludes the possible "bang-pause-bang control” when some judicious slew rate 
limits are imposed on the slew design. 


Analytical and numerical studies were then conducted on the implicit tran- 
scendental nonlinear expression initially provided by NASA Langley Research 
Center for SCOLE's line-of-sight error. A designer's choice of allowable in- 
itial alignment to take advantage of the low moment of inertia in the roll, as 
suggested by Taylor [1], was determined directly analytically. The slew 
angles to achieve the desired LOS pointing were thus determined. [4]-[5] 


A computer program for SCOLE 's complete 3-axis rigid-body dynamics was 
developed and used to simulate numerically the application of various time- 


* The Space Shuttle, while in orbit, is under the single-axis "phase-plane 
rigid-body attitude control of "Digital Auto-Pilot” (DAP). If the two-stage 
approach is applicable, then the current DAP can be used conveniently with 
various proposed flexible-structure flight experiments in space without having 
to make a major specific change in operation or design to suit each different 
experiment setup. 

** Usually all but 2 to 3 modes of the structure were ignored. 
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minimized bang— bang type attitude slew maneuvers. The numerical simulation 
test results indicated that the single-axis bang-bang or bang-pause-bang slew 
maneuvers work fairly well for pointing the LOS of SCOLE under the specified 
conditions. In particular, applying a maximum allowable control moment (i.e., 
10,000 lb-ft) on the Shuttle and a maximum allowable control force (i.e., 800 
lb) on the Reflector, plus imposing 5 deg/sec slew rate limit on the design, 
yields the beat pointing accuracy (0.097°) with minimized slew time (3.733 
sec) and least sensitivity to nonzero products of inertia. Such is a best 
design for LOS pointing slew maneuver for the SCOLE configuration so far as 
the Stage 1 is concerned. [ 5 ] 

For designing vibration control systems (the Stage 2), a standard choice 
would be to apply the modern control and estimation theory, namely, the 
linear-quadratic-Gaussian (LQG) state-feedback control technique. The LQG 
technique has been well accepted because of its success in various other 
applications. Control spillover and observation spillover, however, have 
surfaced as major roadblocks to successful application of such a state-of-the- 
art state-feedback design technique to control vibrations in large flexible 
space structures. Current spacecraft and many other engineering systems on 
which the LQG technique has been very sucessful are of the rigid-body type 
that do not have as many closely spaced low-frequency vibration modes as 
there are in a future large flexible space system. Earlier, Balas [8] showed 
by an example and Herrick [9] followed by a hardware experiment that, because 
of control and observation spillover, even a simple flexible beam, which was 
initially stable in the open loop, became unstable when the "modern modal 
control” loops were closed. 

On the other hand, dynamic properties of large flexible space structures 
can be enhanced by active augmentation of modal damping and stiffness through 
proper output— feedback control [10]— [24]. Lin [20]— [23] showed analytically 
that an appropriate output feedback control system, particularly when it is of 
the type of "modal dashpots" and/or modal springs [21], can even ensure 
full-order closed-loop asymptotic stability of a very general class of 
lightly damped large flexible space structures while improving their dynamic 
characteristics. 

For Stage-2 design to damp the excited vibrations in the SCOLE configura- 
tion, one can consider using a modal-dashpot type of output feedback control 
system first. One may then consider using a modal-dashpot augmented LQG 
optimal state feedback control system, if the LOS stabilization performance is 
not enough to satisfy the specified stringent accuracy requirements. 

Before proceeding to designing a vibration control system for SCOLE, many 
technical issues need to be addressed. For example, one needs to characterize 
SCOLE’ s vibration modes with respect to (i) the excitation by the rapid slew 
maneuvers, (ii) their contribution to the vibration (jittering) of SCOLE's line 
of sight, and (iii) the control authority of the control actuators. Which 
modes need to be controlled directly? Which modes only need some additional 
damping? Which are more likely to cause serious control spillover if not in- 
cluded as "modeled modes"? Those are among many technical questions one gen- 
erally should look into before starting out a meaningful design for SCOLE's 
vibration control system. 


163 



2. MATHEMATICAL MODEL OF THE ORBITAL SCOLE CONFIGURATION 


To assist our quantitative assessments of the vibratory impact of rapid 
bang- bang slew maneuvers on the flexible SCOLE configuration and the perfor- 
mance of proposed vibration control designs, we have developed a computer pro- 
gram to simulate various vibratory responses of the configuration. The com- 
puter simulation was based mainly on the modal data set D3D585 provided by 
Dr. Suresh M. Joshi of NASA Langley Research Center as the flexible-body dyn- 
amics, and the nonlinear LOS error expression formulated by Mr. Larry Taylor 
[1]. We extended a portion of the expression by including a few more terms to 
take a better account of the effect of bending in the mast. 


2.0 Outline of the Orbital Shuttle-Mast-Antenna Configuration 

As shown in Fig. 1-1, the configuration of the SCOLE represents a large 
antenna attached to the Space Shuttle Orbiter by a flexible beam as the Mast. 
The configuration was chosen for its similarity to proposed space flight exper 
iments and various space antenna systems. 

The dynamics of the SCOLE configuration are described [1] by a distributed- 
parameter beam equations with rigid bodies in the three-dimensional space, each 
having mass and inertia at either end. One body represents the space Shuttle 
Orbiter, having the mass, inertia, and dimensions typical of the real one. The 
other body is a large antenna reflector. The equations of motion for the com- 
plete configuration are formed by incorporating the three-dimensional rigid- 
body equations into the partial differential equations of beam bending and tor- 
sion. The flexible mast is treated as a standard slender beam. The boundary 
conditions at the ends of the beam contain the forces and moments applied to 
the rigid Shuttle and reflector bodies. The mast is not attached to the 
mass center of the reflector, but rather significantly away in both x and y 
directions. The nonlinear kinematics of the two sizable bodies and the offset 
attachment of the reflector couple the three otherwise uncoupled beam equa 
tions. The reader is referred to Taylor and Balakrishnan’s paper [1] for the 
details. The rigid-body part of the mathematical model was used by Lin ear- 
lier in his studies on the LOS pointing (i.e., the Stage 1) of the configuration. 
The studies on vibration control reported here were based on a most recent 
version of the flexible-body part available; see Section 2.1 below. 

The line-of-sight (LOS) error of the SCOLE configuration is a highly nonli- 
near implicit expression. The line of sight is defined by a ray emited from 
the feed on the Shuttle which is reflected at the center of the Reflector. It 
is affected by the pointing error of the Shuttle, the offset attachment of the 
Reflector, and the misalignment due to the deflection and torsion of the Mast. 
The reader is again referred to Ref. [1] for the original formulation of the 
LOS error. An equivalent expression having a simple modification, which is 
more convenient than the original for both efficient numerical computations 
and in-depth analytical investigations, was used by Lin in his earlier rigid 
body studies [ 3 ]- [ 5 ] . For the current flexible-body studies, the nonlinear LOS 
error expression also needs some more terms in order to have a better 
accounting for the bending of the mast beam; see Section 2.2 below. 
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2.1 Flexible-Body Dynamics 


The bending and torsion characteristics of the SCOLE configuration were 
originally formulated in partial differential equations by Taylor and Balak- 
rishnan (IJ. Robertson [ 2 5 J derived the corresponding equations of free motion 
taking into account the kinematic coupling resulted (i) from the offset attach- 
ment of the Reflector to the Mast and (ii) from the nonzero products of 
inertia of both the Shuttle and the Reflector. He then solved the equations in 
terms of trigonometric and hyperbolic functions and computed a set of natural 
frequencies and mode shapes. Such results are not readlily useful for control 
studies. 

To facilitate the control analysis and design for SCOLE, Joshi [26] first 
derived a state-space model from Robertson's results; the data set was named 
"BMDT3D". Later, he improved Robertson ? s results, and also derived another 
state-space model, named M D3D585 M . Our computer simulation program was in- 
itially based on the data set BMDT3D, which contained only the first five flex- 
ible-body vibration modes. It was then updated when Dr. Joshi furnished us 
with the set D3D585 later. 

The set D3D585 provides modal data in the state-space (A,B,C) form. It 
contains only the first 10 flexible-body modes but no rigid-body modes nor any 
nonlinear rigid-body dynamics. This set was quite appropriate for our purpose 
of assessing the vibratory impact on SCOLE. We found it more effecient and 
convenient, however, to compute the time transition of the states using the 
second-order modal equations directly, because of the decoupled nature of the 
former, than to do so using the first-order state equations. We thus con- 
verted the furnished data back to the following standard modal form: 


m 

n i + ^ = 1>***»N (2—1) 

k-1 

n 

y_j - I< c vM + VjVP j = (2-2) 

i=l 

where = 2 £.&>., = Uj. 2 (2-3) 

are, respectively, the damping and stiffness coefficients of the unit-mass 
linear oscillator representing the ith vibration mode; and denote the 

natural frequency and mode shape, respectively, of mode i; qi denotes the in- 
herent damping ratio of mode i, which had been assumed to be 0.3% for all 
flexible-body modes of SCOLE [1]. m and n i denote the coordinate and velo- 
city, respetively of the ith mode. 

The kth force (torque) Input is denoted by u^, with column vector bp^ rep- 
resenting the corresponding actuator influences on SCOLE. The jth measurement 
output is denoted by yj, with row vectors Cyj and Cpj representing, respec- 
tively, the velocity and displacement sensor influences. 

Putting (2-1) and (2-2) into a matrix form, we get 
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(2-4) 


n + An+Eri = ^Bp,u 

y = C v *fi + C D *n (2-5) 


where 


A = diag [ 2 C i w i ]» 1 



( 2 - 6 ) 


■» *n ] ’ 


[b 


FI* 


b FmJ 



(2-7) 


( 2 - 8 ) 


In accordance with Robertson's formulation, we also assume that the bend- 
ing and torsion in SCOLE are referred to the coordinate system defined on its 
initial undeformed configuration. Thus, before any deformation, the center of 
mass of the Shuttle is at the origin of the coordinates; the roll, pitch, and 
yaw axes (i.e., body x, y, z axes) of the Shuttle, align with the x, y, z coor- 
dinate axes* respectively; and, in particular, the straight mast beam coincides 
with the z coordinate axis. Note that, since the flexible mast was not tre- 
ated as a cantilevered beam in Robertson's derivation, not only the mast may 
not be tangential to the z coordinate axis, but the center of mass of the 
Shuttle also may not remain at the origin, nor may the Shuttle body axes 
remain parallel to the coordinate axes, when a significant deformation of the 
mast occurs. The line of sight of the SCOLE configuration will thereby be sig- 
nificantly affected. 


2.2 Line-of-Si ght Error Expression with More Bending Terms 

In order that the jittering of the line of sight (LOS) due to excited vibra- 
tions can be more accurately evaluated, we used almost the same nonlinear 
expression for the line-of-sight (LOS) error of the SCOLE configuration as ori- 
ginally given in Ref. [1]. Unlike the orginal, however, our improved version 
also takes into account the z— axis dislocation of the Reflector due to 
bending of the mast, and like the one Lin used earlier [4]-[5] the LOS vector 
R l03 is not normalized. Note that the LOS error expression could be expanded 
in a Taylor series and a linearized version could be obtained by taking the 
first-order terms. A linearized version, though useful in linear-quadratic 


* Robertson's y and z axes are opposite in sign to those defined by Taylor for 
SCOLE. We continue to adopt Taylor's definition for consistency with Lin's 
earlier studies [ 2 ] — [ 5 J on the Stage 1. 
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optimal control designs, is not appropriate for our current use since the 
excited vibrations are sufficiently large in magnitude and the second- and 
higher-order terms of the series expansion may not be negligible at all* 

The location of the center of the Reflector, represented by Rr ,is defined 
by the location of the joint where the Reflector is attached to the Mast; see 
Figure 2-2. Denote by Rj the location of the joint relative to the center of 
the Reflector, and by R^ the location of the same point (also the tip of the 
mast) with respect to the center of the Shuttle. Then the vector R^ is given 
by 


R R ■ *T - T 1 T 4 R J 


(2-9) 


where Rj 


-18.75 

32.5 

0 


( 2 - 10 ) 


The vector is constant in magnitude because of the rigid reflector, but its 
orientation with respect to the Shuttle is affected by the deflection at 

the tip J. The product Ti^T 4 of coordinate transformations T^ and T 4 is to 
take care of the angular change. As in Ref. 1, T^ denotes a direction-cosine 
transformation from the Shuttle to the Earth (inertial) coordinates, and T 4 one 
from the Reflector to the Earth coordinates. 


A reasonable approximation for the tip location is given by 

Bend x 

Bend v 

R T a 

-/ 130 2 - Bend x 2 -Bendy 2 


( 2 - 11 ) 


where Bend x - u xg - u xR Bend y = u yg - u yR 

u x s and Uyg denote the deflections of the mast at the Shuttle end in the 
xz and yz planes, respectively; u x r and UyR are the corresponding deflections 
at the Reflector end. 

Eqs. (2-9)-(2-ll) constitute our additional modification to the LOS error 
expression. Note that the vector R R originally given as (18.75, -32.5, -130) in 
Ref. [1] corresponds to the undeforned case. To see it, assume that ther 

are no deflections at all. Then Bend x and Bendy are both zero, and T 4 is 
equal to T^. Therefore, R^ » (0, 0, 130). Consequently, Rr * Kj -Rj * (18.75, 
-32.5, 130). 

In the analytical studies on the LOS error of SCOLE [4]-[5], we found it 
more convenient not to normalize the LOS vector R^qs first, although the res- 
ulting error expression is the same since division by its norm is still made 
later at the end. The LOS error with such a trivial modification is given by 

e L0S ‘ * sI “' 1 [ I I“t * T 1 R I-0S| I / I I*L0SI I ] (2-12) 
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with the unnormalized LOS vector being defined as 


r los ■ RF '- r r - r f - 2 1< r r - r f>- r a1 r a «- 13 > 

Where as defined in [1], Rp is the vector representing the feed location (3-75, 
0, 0). is a unit vector in the direction of the Reflection axis in the Shut- 
tle body coordinates, i.e.. 


R 


A 


tT j 


0 

0 • 

1 


For the target direction specified in [1] as Dj * (0, 0, 1), Expression 

(2-12) reduces to 


'LOS 


- ± sin 


-1 


/ (T lrl + ^ T lr2 ^LOS^ ^ ^ 


(2-14) 


where and Tir2 denote respectively the first and the second rows of 

matrix T^. 
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3. VIBRATORY RESPONSES TO BANG-BANG TYPE RAPID SLEW MANEUVERS 


Several LOS pointing slew maneuvers of the bang-bang type were applied to 
our computer simulation of the SCOLE flexible-body dynamics. The resulting 
responses range from excessive to minimal, depending on the magnitude of the 
applied force at the Reflector. Note, however, that all these slew maneuvers 
were designed to provide minimized slew time under the increasingly tight limit 
imposed on the respective applied force. 

The slew maneuver that excited the most violent vibrations in SCOLE was 
chosen for studying the control design and for generating in-depth insights 
into the vibration control challenges. On the other hand, the least violent 
one deserves further exploration in the future, since it may potentially 
require a smaller total time for both slew and stabilization . 

In assessing the impact of structural vibrations on SCOLE, we view the 
slew maneuvers as time-dependent disturbances instead, and only the vibra- 
tory portion of the time-domain responses are of real interest. Therefore, t 
is reasonable that we concentrate only on the flexible-body and temporarily 
ignore any rigid-body dynamics in this study. This assumption is equivalent to 
the absence of rigid-body dynamics. It is also reasonable to assume that, 
before being subject to such disturbances, SCOLE was initially at rest and a 
no deformation nor LOS error . The former assumption is equivalent to setting 
to zero the initial conditions on the normal coordinates and velocities of all 
modes, and the latter equivalent to aligning the undeformed SCOLE configura- 
tion with the attitude (♦i,0i,*i) that corresponds to zero LOS etror. Such 
roll-pitch-yaw Euler angles, calculated and used by Lin earlier [4] [5], are 

listed below for reference: 

4^ - -14.03624347°; = -6.38707294°; <1^ = 0°. 


3.1 Excitation by the Rapid Time-Minimized Bang-Pause-Bang Sl ew Maneuver 

We first examined, through numerical simulation, the SCOLE flexible-body 
dynamics under the excitation of the rapid time-minimized roll-axis bang- 
pause-bang (BPB) slew maneuver that was considered a best candidate for 
pointing the line of sight of the SCOLE as a rigid body [4]-[5]. Among many 
other single-axis LOS pointing slew maneuvers of the bang-bang type previous y 
studied, this BPB maneuver was judged to be the best compromise in terms LO 
pointing accuracy achievable, slew time required, and performance robust 
ness to nonzero products of inertia. It was designed to slew the SCOLE confi 
guration about the negative roll (i.e., -x) axis for about 20 to correct the 
initial 20° LOS error specified in [1]. This slew maneuver requires that the 
maximum allowable moment (10,000 lb-ft) be applied to the Shuttle about the 
negative roll axis and simultaneously the maximum allowable force (800 lb) at 
the Reflector center along the negative y axis, both for only 0.86 sec *> 
then, after a long pause of 3.158 sec., these maximum moment and force be 
applied again for only 0.867 sec. but in the opposite directions (i.e., positive 
roll and y axes, respectively). 


Such a BPB slew maneuver was applied to our computer simulation of t e 
SCOLE flexible-body dynamics. The simulation results are summarized by the 
plots in Fig. 3-la, which show that such a maneuever would cause excessive 
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Fie. 3-1 Vibratory responses to Rapid Time-minimized Bang-Pause-Bang Slew; 

b. Attitude deviations at the Shuttle (S) and the Reflector (R) ends 
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vibrations in SCOLE! Observe that: the line of sight vibrated with error 

between 89.8° (or 133.3° if not taking on the principal value of the arcsine in 
Expression (2-14)) and 14.7°; the tip of the mast vibrated in the yz-plane 
between +114 ft and -113 ft. 

Fig. 3-lb show the deviations in Euler attitude angles of the Shuttle (S) 
and the Reflector (R) from their "nominal** alignment of zero LOS error. These 
deviations correspond to the bending slopes and the torsion at the respective 
end of the mast. Observe that the Shuttle rolled to the right and the left 
between +17.16° and -17. 06°, while the Reflector rolled to the left and right 
between -86.96° and +88.35°. There were virtually no pitch and yaw motions of 
the Shuttle, but the Reflector pitched between -10.63° and 8.75° and yawed 
betweeen -32.27° and +27.97°. 

In general, some significant excitation of the vibration modes of a flexi- 
ble space system, such as the orbital SCOLE configuration, should be expected 
when large moments and forces were used to their limits in a bang-bang manner 
to minimize the slew time. The appalling magnitude of the vibratory impact, 
however, was indeed a surprise. 

Such excessive vibrations certainly post serious challenges to the Stage-2 
design, i.e., the control design for suppressing such vibrations after the exci- 
tation. Can such large-magnitude vibrations be brought down to some tolerable 
level in about the same length of time (say, 5 sec.) as the slew maneuver? 
How to design such a fast effective vibration controller? We shall continue 
to address such design challenges in Section 4. 


3.2 Excitation by Other Rapid Time-Minimized Bang-Bang Slew Maneuvers 


Are all slew maneuvers of bang-bang type so terrible to flexible space 
systems? Why are the excited vibrations in SCOLE so large in magnitude? Even 
when one can design a powerful fast vibration controller capable of damping 
out such vibrations, one still cannot stop thinking of these and other puzzling 
questions* To investigate further, we conducted the following numerical 
experiments on our computer simulation of SCOLE flexible-body dynamics. All 
were the same as before, except that a different bang-bang slew maneuver was 
applied. 


3.2.1 Experiment F10 — No force on Reflector. First we tried to use 
only the 10,000 lb-ft moment on the Shuttle. The same roll-axis bang-bang 
slew maneuver using only such a moment for accomplishing the same 20° pointing 
task in the minimum time as was previously designed and evaluated on the rigi- 
dized configuration in [ 4 ] — [ 5 J was tried. This maneuver requires that the max- 
imum moment be applied first about the negative roll axis for 6.307 sec, and 
then switched to the opposite directions (i.e., positive roll axis) for another 
6.307 sec. It was truly a bang-bang (BB) control. 

The simulation results, as shown by plots in Fig. 3-2, clearly show that 
the vibratory impact was greatly reduced. The LOS error was only 6.25° at 
most, and the mast tip vibrated only between +5.06 ft and -5.18 ft. 

Of course, the (minimized) slew time is much longer; it is a main reason 
why this maneuver has been rejected earlier [4]-[5] as a Stage 1 design for 
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SCOLE. This simulation is useful only when its results are compared with the 
foregoing case of using additional 800 lb force on the Reflector: it serves as 
an opposite extreme, since no force was applied to the Reflector at all. 

By a careful inspection of the time histories of the tip deflection in both 
(see Figs. 3— lb and 3— 2b), we can make the following interesting observe 
tions. While the moment was being applied to the Shuttle about the negative 
roll axis without any force on the Reflector, the beam bent backwards and the 
Reflector lagged behind*. On the contrary, the addition of the maximum force 
on the Reflector reversed the situation, even though the additional force had 
exactly the same purpose of rolling the configuration to the same side as the 
moment on the Shuttle! The Reflector then became leading instead of lagging. 


3.2.2 Experiment F180 — 80 lb Force on Reflector. The leading of 
the Reflector might be responsible for the huge increase in LOS error, as imp- 
lied by the above observations. It is therefore reasonable that reducing the 
applied force might reduce the lead and hence reduce the LOS error. A second 
experiment was thus conducted with an 80 lb maximum force, which is only one 
tenth of the original allowable maximum. 

A new roll— axis slew maneuver was designed, in the same way as the first 
BPB slew maneuver; but only 80 lb, instead of 800 lb, force was to be used in 
conjunction with the same 10,000 lb-ft moment to accomplish the same 20* LOS 
pointing in a minimized time. It turned out to be a bang-bang maneuver in- 
stead, since the slew rate would not reach the imposed 5 deg/sec limit. In 
almost the same way as in the case of 800 lb, the slew maneuver requires that 
both the moment and the (tighter-limited) force be applied with respect to the 
corresponding negative axes for 4.416 sec, and then reversed to the corres- 
ponding positive axes for another 4.416 sec, but with no pause in between. 

The simulation results, as summaried by plots in Fig. 3-3, confirmed what 
we thought. The lead by the Reflector is now greatly reduced, and so are 
the LOS error and the mast bending, compared to the case of 800 lb (Fig. 3-1). 
The LOS error was only 24.7* at the highest peak of its time history; the tip 
deflected only between +20.59 ft and -10.83 ft; and the Reflector rolled only 
between +15.98* and -8.31*. 

These results have clearly shown that the 800 lb force was direct ly res— 
ponsible for the excessive vibrations and the unreasonable LOS error. 


Next, compare these results with those of Experiment F10 (Fig. 3-2). A 
peak LOS error of 24.7* is fairly large compared to only 6.25° of Experiment 
F10; so is a maximum deflection of 20.59 ft compared to only 5.18 ft of F10. 
Does this mean that no force should be applied to the Reflector at all? No, 
we did not think so! Instead, we reasoned that if one could reduce the lead 
slightly further, one could further reduce both the LOS error and the tip def- 
lection. So a third experiment with a slightly smaller force was performed. 


* Note that when a negative moment is applied to the Shuttle, a positive def 
lection indicates the lagging of both the mast tip and the Reflector. 
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pig. 3-3 Vibratory responses to Rapid Time-minimized Bang-Bang Slew: 80 lb; 

b. Attitude deviations at the Shuttle (S) and the Reflector (R) ends. 
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3.2.3 Experiment F125 — 25 lb Force on Reflector. Since the Shuttle- 
attached SCOLE configuration was chosen because of its similarity to proposed 
space flight experiments [27]-[29], we thought it would also be more realistic 
to consider a force of about the same level as the vernier RCS thrusters on- 
board the Shuttle Orbiter. Since the existing vernier thrusters generate 24 to 
24.5 lb thrust each [30], we simply selected 25 lb for the third experiment. 

Again a new roll-axis slew maneuver was designed in the same way as 
before for accomplishing the same 20° LOS pointing in the minimum time. It 
certainly is a bang-bang maneuver, like the case with 80 lb force. This BB 
slew maneuver requires that both the 10,000 lb-ft moment and the 25 lb force 
be applied with respect to the negative roll and y axes, and then switched to 
the positive axes, as before, for 5.479 sec each time. 

The results, as shown in Fig. 3-4 by plots, are very pleasing, indeed. The 
largest LOS error was less than 0.51° ; the tip deflected only between 40.25 
and -0.3 ft; and the Reflector rolled only between 40.16° and -0.3° "i All are 
one order of magnitude smaller than those from applying no force on the 
Reflector! Of course, the time required for completing the 20° slew of the 
line-of-sight is also shorter. In summary, for a BB slew maneuver of the 
flexible SCOLE configuration, using a force of 25 lb on the Reflector in 
addition to a 10,000 lb-ft moment on the Shuttle is In all aspects supe- 
rior to using no additional force there. 

The force of 25 lb is simply a rather arbitrary trial value. One could 
continue to search for an optimal value that would result in still smaller tip 
deflection, but we did not do so because we felt that our original purpose had 
already been served very well. 

If LOS error were the only concern and time were not so important, then 
one should immediately stop studying the use of 800 lb force on the Reflector. 
On the other hand, since time is at least equally important for SCOLE, it is 
not clear at all that 25 lb might be preferred outright to 800 lb: the mini- 

mum time required for the same 20° slew is 10.959 sec for the case of 25 lb 
but only 4.892 sec for the case of 800 lb, that is, more than twice longer. 
Moreover, in both cases, some active vibration controllers are still needed to 
damp out the excited vibrations; and hence some additional time is required in 
order that the required LOS accuracy of 0.02 6 can be met. 

To damp out excessive vibrations, such as excited by the BPB roll-axis 
maneuver using both an 800 lb force and a 10,000 lb-ft moment, can be serious 
challenges to the Stage-2 control design. Insight and techniques generated 
from dealing with such challeanges certainly will be useful in designing 
effective vibration controllers for the case of using a smaller force, such as 
25 or 80 lb. 


199 





ATTITUDE DEV. DEG ATTITUDE DEV. DEG 

.80 - 0.40 0.00 0.40 „- 0.04 0.00 0.04 


00 

CO 


f > 


ROLL, S 



PITCH, S 


TRW, S 

X K )()()( 


. 00 


4 . 00 


8. 00 

TIME 


12.00 


SEC 


16.00 


F 1 25 

20. 00 



TIME SEC 


Fig. 3-4 Vibratory responses to Rapid Time-minimized Bang-Bang Slew: 25 lb; 

b. Attitude deviations at the Shuttle (S) and the Reflector (R) ends. 
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4. ACTIVE VIBRATION CONTROL FOR SCOLE 


4.1 Direct V elocity-Output Feedback Control 

Let yy denote the velocity-sensor outputs. Then from the system Eqs. 
(2-4)-(2-5), we have 

y y * CyOn (4-la) 

The general form of direct velocity-output feedback control is 

u - - Gy y (4-lb) 

where G denotes a matrix of constant feedback gains. Substituting (4-1) in 
(2-4) results in the following closed-loop system 

rf + (A + * t B f GC v $) n + Zn » 0 (4-2) 

The modal stiffness matrix Z of SCOLE flexible-body dynamics is positive 
definite, since no zero-frequency rigid modes were included in dataset D3D585. 
By applying the classical Kevin-Tait-Chetaev theorem, its extensions, or Liapu- 
nov's second method, one can show (see, e.g. [17]— [23]) that the closed-loop 
system (4-2) is: 

(I) stable (in the sense of Liapunov) if the augmented damping matrix 
(A +$^BpGCy$) is symmetric and nonnegative definite, and 
(ii) asymptotically stable if the augmented damping matrix is positive definite. 

When the velocity sensors are, as generally assumed, co-located with the 
sensors, i.e., Cy * Bp T , the additional damping matrix ($ T B F GCy$) is always 
nonnegative definite*, whether the gain matrix G is positive or merely non-neg- 
ative definite. In other words, direct velocity-output feedback control at 
least will never destabilize the system, even when no inherent damping exists 
(i.e., A * 0). 

For most practical cases where there are less actuators than vibration 
modes and there are virtually no inherent damping (i.e., A is small and some of 
its diagonal elements are virtually zero), the existing theory cannot help det- 
ermine whether a closed— loop system is asymptotically stable or not, though 
numerical results can [17]. 


The theory is not enough to help design the feedback gains, either. Usu- 
ally designers simply restrict the gain G to be diagonal matrix, and therefore, 
make each co-located pair of actuator and sensor act like no more than a (pas- 
sive) dashpot. Having no systematic method to help calculate the required or 
desirable values for the feedback gains, some designer even set the diagonal 
elements rather arbitrarily to some trial positive numbers. A practical ques- 


* Only in a rare special case, which is rather unrealistic to truly flexible 
large space systems, where there are as many independent actuators (and co- 
located independent sensors) as there are vibration modes and the influence 
matrix Bp is nonsingular, will a positive definite gain matarix G quarantee 
that the product S^BpCy^ is also positive definite. 
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tion is: how to design the gains so as to add more damping selectively to some 
specific modes than others? How not to restrict the resulting design to be 
strictly local feedback? How to design the direct velocity-output feedback as 
a really multi-variable control system? A systematic design method is needed. 


4.2 Concept of Modal Dashpots 

The diagonal form of feedback gain matrices spreads the control effort 
thin over all the vibration modes. One cannot design the diagonal form for 
adding desirable amounts of damping respectively to certain selected vibration 
modes. On the other hand, when one wishes to add a certain amount of damping 
to each mode, one might consider computing the gain matrix G as the general 
solution of the following NxN matrix equation 

$ t B f G C = A* < 4-3 > 

without restricting it to be diagonal, where A* denotes the matrix of desired 
additional modal damping. Expressed in terms of additional damping ratio ^ 
desired of each mode, the matrix A* may take on the same simple form as Eqs. 
(2-6), (2-3), i.e., 

A* = diag£ 6* J (4-4a) 

and 6* * 2c^u^ i = 1, •••, N (4- 4b) 

Note that for a realistic flexible space structure there are much more vibra- 
tion modes than there are locations for placing actuators or sensors (i.e., 
N»* and N»m). Thus, if one wishes to augment some indeterminate amount of 
active damping to all the modes, then one may try to obtain an approximate 
solution of Eq.(4-3), such as of the least squared error like the following 

G P - (.*»/ A* (C v ») P <«-5) 

where A P denotes the Moore-Penrose pseudo- inverse. No conditions on the 
matrix or Cy^ need to be satisfied, and the pseudo— inverses can be ealeu 

lated numerically using the singular value decomposition [31] -[32]. 

Solutions of the form (4-5) have three major practical drawbacks. First, 
the number N of vibration modes in a realistic flexible structures is enor- 
mously large, making it impractical, if not impossible, to calculate the 
pseudo-inverses of the extremely large matrices * T B F and C v *- Secondly, one 
still cannot really focus a specific subset of the modes, since the solution G 
is merely a least— square approximation, with errors spread all over the 
modes. Thirdly, also because of approximation errors, the resulting product 

* T B F g p C v * 

might not be symmetric, and hence stability might not be guaranteed. 

In practice, one needs to concentrate on a relatively small number of 
important modes. In many cases, one cannot care less for those modes which 
are less important when one cannot even get what is required for suppressing 
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the more important ones* Thus* assume that some n modes (n<<N) are the most 
important, and a reduced-order dynamic model is formed by selecting only those 
n modes* Call those modes modeled modes and the rest unmodeled modes * Par- 
tition the matrices n, ♦, A, and Z accordingly into the modeled (M) and the 
unmodeled (U) parts, i.e*, 



* - A * black “ dia 8t A M ,A U^ 1 " black-diaglEn.Zjj] (4-6) 


Then, the closed-loop equation for the reduced-order model is 


"" M + + ♦Svsv + Vm ■ 0 <w> 

Now, let a reduced matrix be given that corresponds to the desired 

additional damping for the n modeled modes. Then the design is reduced to 
solving the following much smaller nxn matrix equation, instead of the NxN Eq. 
(4-3), for the gain matrix G: 


♦S b f g Vm 


M 


(4-8) 


As before, a solution in the same general form as (4-5) can be obtained numer- 
ically by computing the pseudo-inverses of influence matrices ($M T Bp) an< * 
(C v *m)* It is still an approximate solution unless some rank conditions are 
satisfied by the influence matrices. 


Of the particular interest is when the control influence matrix (* M T Bp) 
has the full row rank and the observation influence matrix (Cv$m) has che 
full column rank. In other words, 


rank('t^Bp) ■ row(#^Bp) ■ n (4-9a) 

rank(Cy$^) = colunn^Cy*^) = n (4-9b) 


Such a special case requires that n<i and n<m, i.e., the number of modeled 
modes do not exceed both the number of actuators and the number of sensors to 
be used in the feedback control. Under the full-rank conditions (4-9), the 
pseudo-inverses are also generalized Inverses. That is. 


<*MV f 




) R 




T 

right generalized inverse of 

<*5v T [“Sv< 4 5 b f> t ] 1 

left generalized inverse of 


1-1 




( C V *M^ 


(4-10) 


(4-11) 
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The gain matrix G computed therewith solves Eq. (4-8) exactlyt. The closed- 
form expression is given by 


(; * = a m ( Vm )L 


(4-12) 


Consequently, the reduced-order closed-loop system equation (4-7) thereby 
simplifies to 

+ + *M + = ® (4-13) 

The desired damping is thus added ^to the reduced-order model exactly as 

specified. For stability, the matrix A^ of additional damping only needs to be 
nonnegative definite; it need not be diagonal. 


When Am is chosen to be a diagonal matrix, as it is often convenient and 
reasonable to do in practice, the resulting velocity-output feedback control 
will perform like a separate "dashpot" attached to each mode of the 
reduced-order model. Specifically, let 


A* = diag 



(4-14) 


Then (4-13) can be rewritten in the component form like (2-1) as follows: 


n Mi + 


( 6 Mi + 


* x 

6 m J 

Mi' 


'Mi 


“Mi n Mi = ° 


(4-15) 


where ^i denotes the normal coordinate of the ith modeled mode. Eq. (4-15) 
obviously means that the ith modeled mode, like an independent linear oscil- 
lator, is augmented with an additional dashpot whose damping coefficient is 
5 Mi* This is why Canavin called such s design a "decoupled controller”, or 
"modal dashpots" [11]. 


The diagonal elements 6^ should be nonnegative to jiake a practical sense. 
Like (4-4), it can be given in terms of damping ratios and natural frequen- 

(4-16) 


cies as 


Mi 


* 

p Mi 


2 CwJ w 


Mi 


t Given whatever values to the matrices, it is mathematically an exact solu- 
tion so far as the equation (4-8) is concerned. Of course, it may not be an 
exact solution so far as the system (4-2) or even (4-7) is concerned, when any 
matrix, for example the modal matrix $ as usual, contains some modeling or 
computational errors. Small errors in 4> may invalidate stability results of 
general feedback gains but not the modal-dashpot type [20]-[23]. 
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4.3 Improvement on the Design Method 


One can design very effective vibration controllers by the method of modal 
dashpots, as demonstrated by our applications to SCOLE. The interesting simple 
formula (4-10)-(4-12), however, does not by itself complete the design 
method for an effective control of structural vibrations. In fact, when 
the concept of modal dashpots was initially formulated by Canavin [10] — [12] as 
"decoupled controller", it was accompanied by two major technical drawbacks 
that almost rendered itself practically useless. Later, through various 
numerical evaluations and theoretical analyses, Lin and his associates 
[ 20 ] — [ 24 ] identified the underlying causes of these problems, and greatly 
enhanced the utility of this concept. In the course of applying it to the 
challenging SCOLE vibration control design problem, we also made some addi- 
tional improvement on this design method. 


A first initial technical drawback was the high-gain low-damping problem. 
After he applied it to a representative large space structure (of which 37 
vibration modes were considered), Canavin concluded that "the decoupled con- 
troller may be of limited utility due to the high gains produced by this 
approach "[11]. The feedback gains were mostly in the orders of 10 li to 10 ia , 
while only additional 10% of critical damping was designed for each of the 12 
modes he had selected to be "controlled" (i.e., modeled) modes. 

Aubrun [13] proposed the approach of low-authority control (LAC) by limit- 
ing to 10% modal damping and by using sufficiently small gains so that the 
amount of active damping achievable is predictable. Since then, direct velo- 
city-output feed back control has been commonly thought to be of only low 
authority, low performance, and secondary importance. However, the vibration 
controllers of Aubrun* s design should be of low authority, not because of 
direct velocity-output feedback, but rather because of the applicability of 
Jocobi's root perturbation formula on which he based his theory. For his use 
of the perturbation formula to remain valid, the control authority (and 
specifically the feedback gains) must be sufficiently low so that the closed- 
loop eigenvalues and eigenvectors would be resulted from only infinitesimal 
perturbations, i.e., only very small increase in damping ratios. 


A second initial technical drawback of the basic design method was severe 
interactions between modeled and unmodeled modes. When the method was 
applied to another representative large space structure (i.e., ACOSS model 2 
[33]), the interactions were so severe that the desired damping performance on 
the modeled modes was degraded very badly, although the closed-loop system 
remained stable [24]. 

The following common causes were discovered. 

(1) Some modeled modes had too small control influences or too small 

• observation influences (Cv<t>Mi) # This made the generalized inverses ($M T Bp) R 

or (Cy$ M ) L , and hence the resulting gain matrix, unnecessarily large. 

Theoe low-influence modes should he deleted from the reduced-order model, 
or else some actuators or sensors should be relocated to improve their in- 
fluences on these modes. 

(2) Some of the rows in matrix <fy T Bp had too small degree of independence from 

the others, or some columns of Cy% had the similar situation. This also 
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(3) 


made the generalized inverses, and the gain matrix, unnecessarily large in 
magnitude. Like (1) above, these modes should be excluded, or the location 
of some actuators or sensors be improved. 


Some Low-frequency umnodeLed modes had too large control influences 
(d> n t B p ) or observation influences (C v <|>uj) compared to those of the mode e 
modes. This made excessive spillover. These modes should be added to the 
reduced-order model; otherwise, some actuatore or sensors should be relo 
cated, or their influences be properly synthesized [35}-[36] - 

Some of the desired additional damping coefficeints (« M i*) were to ° large 
for some modeled modes, even all were set equal to the same small design 
value (say Cmi* = 0.1). This made some part of the gain matrix unneces 
sarily large, and hence increased interactions with some unmodeled modes. 


Open-loop responses of individual modeled modes should be analyzed and the 
need for additional damping realistically guesstimated with respect to the 
control/observation influences on each modeled mode. For properly desdgtl ® d 
modal dashpots, e.g., our design for SCOLE, the additional damping c °" ld * “ 
high as 67% for some modes or as low as 3 % for some others, depending 
ability of the actuators as well as on the individual open-loop responses. 


We have begun to develop the concept of modal dashpots into a useful sys- 
tematic design method for direct output feedback vibration control. Although 
the closed-form formula has reduced the design of modal dashpots to simple 
cranking of numbers, yet to make it really work for effective control of iarge 
excited structural vibrations in flexible space systems, such as the SCOLE con 
figuration, many careful pre-design steps have to be taken. 

The design method was initially formulated by Canavin without explicit 

consideration of limitations on the requirement for c ° ntro1 fl orce * !°J2 olE* 

Now the explicit limits must be considered when applying the method to SCOLE. 

Al”o some saturation "circuitry- must be Imposed on the feedback control so 
that’ the magnitude of the forces or moments generated by the modal dashpots 
tould automatically be limited to 800 lb and 10,000 Ib-ft, ^pectively. 
Saturation may not destroy stability when actuators are co-located *th sen 
sors [37], but would somehow limit the performance of the feedback con 

troller. 
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5. DESIGN OF MODAL DASHPOTS FOR SCOLE 


The vibrations in the SCOLE configuration excited by the rapid time-minim 
ized BPB LOS pointing slew maneuver, as reported in Section 3.1, posted three 
serious vibration control design challenges: 

(1) The excited vibrations were excessively and unrealistically large in magni- 
tude: the line of sight once had an error of 89.9° (or beyond) and the 

130-ft mast once had a tip deflection of 114 ft. 

(2) The allowable time was extremely short: it should be minimized, so only 

an equally short time (specifically, only 5 sec, which was approximately 
equal to the maneuver time) was allowed. 

(3) The available control forces and moments were limited: the 800 lb and 

10,000 lb-ft limits were imposed the same way as on slew maneuvers. 


In order to design effective modal dashpots for suppressing such excessive y 
large vibrations in SCOLE in a very short time, we conducted careful pre design 
analyses on the vibration modes and their influences by the actuators and sen- 
sors. The candidates for modeled modes were selected, and then divided into 
two groups according to the actuator influences. The design of the modal 
dashpots was therefore divided into two parts accordingly. 


5.1 Analysis on Vibration Modes 

Initially, two different numerical analyses of SCOLE vibration modes were 
made each with a different standard measure of importance. The results were 
inconsistent. Then a third measure was developed and used; the results were 
finally fair and satisfactory. 


5.1.1 Measure 1: LOS Error due to Initial Modal Displacement 

"LOS error contribution" is a common measure used by many structural dynami- 
cists for determining if a vibration mode is "critical" or not, i.e., if it nee s 
active control or not. It was used by Draper Laboratory [33]-[34], and 
accepted by other ACOSS* and VCOSS** contractors [38]-[45] as the standard 
approach, in the modal analysis of both Model No. 1 (namely, the Tetrahedron) 
[34] and Model No. 2 [33] of representative large flexible precision space 
structures. The standard approach is to express the LOS error as a linear 
function of physical coordinates under the assumption that all the displace 
ments ate sufficiently small. When the physical coordinates are transforme 
into the normal coordinates of the structure, the LOS error become a linear 
function of the normal coordinates. The "critical modes are then determined 
by comparing the modal coefficients of the LOS error. 

Such a measure is not directly applicable to rapid pointing of the SCOLE 
configuration nor, in general, to large space structures that are su ^® ct 
rapid slew or retargeting maneuvers. First, the displacements (deflections and 
..rsions, for example) generally are large, hence the linearization of the LO S 
error is not valid. Thus, for SCOLE, we used the original nonlinear expression 


* Active Control of Space Structures, a DARPA technology program. 

** Vibration Control of Space Structures, sponsored by Air Force Wright 

Aeronautical Laboratories. 
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without linearization. Secondly, the LOS error is a dynamic vibratory res- 
ponse^ instead of being simply a static displacement of the line-of-sight . 
Thus, Instead of comparing only the LOS error coefficients, we compared the 
time his Lories of the LOS error the individual modes would separately cause if 
they were intially excited alone. 

For this analysis, the SCOLE configuration was assumed to be initially at 
rest with no LOS error*, and only one mode was excited each time because of a 
unit initial displacement in its normal coordinate. Specifically, for the ith 
time history, the initial condition was assumed to be: 

Hi = 1, and n j = 0 for all j ± i ; n j - 0 for all j. 

For each such initial modal displacement, the time history of the resulting LOS 
error was calculated separately using our computer simulation program. 

The results of 10 separate cases (one for each mode) are shown together by 
the overlapped plots in Fig. 5-1, where each curve represents a completely 
separate time history of LOS error. Listed below are the highest peak value 
of each time-history curve. 

Mode: 1234567 8 9 10 

Peak: .37 .53 .54 .93 1.3 .14 ,51 .002 .18 .03 

The relative importance of the 10 modes is thus given in the descending order 
as follows. 

Mode: 5, 4, 3, 2, 7, 1, 9, 6, 10, 8. 

5.1.2 Measure 2: Modal Response to the Rapid Pointing Maneuver 

By intuition, a vibration mode is more in need of active control than others 
when its magnitude of excited vibration is larger. Thus, a second measure of 
importance for the SCOLE configuration naturally is the vibratory response of 
each mode to the rapid pointing maneuver. For this analysis, the configuration 
was assumed, as before (in Section 3.1), to be initially at rest without any 
LOS error or any nonzero initial conditions, and the same BPB slew maneuver 
was the source of excitation. The time history of the resulting modal res- 
ponse ni(t) was calculated for each mode separately . 

The results are shown by the plots in Fig. 5-2, with each curve represent- 
ing an individual mode. Listed below are the highest peak value of the curves. 

Mode: 123456789 10 

Peak: 21.6 603 41.2 13.7 0.49 0.48 0.28 .058 .041 .001 

Accordingly, the relative importance of the 10 modes is thus given by the fol- 
lowing descending order: 

Mode: 2, 3, 1, 4, 5, 6, 7, ... 


* As stated in the beginning of Section 3, we assumed that, before any of its 
vibration mode was subject to excitation, SCOLE was initially at rest and had 
no deformation nor LOS error. Specifically, the undeformed configuration was 
assumed to have been aligned with the attitude angles of zero LOS error. 
Therefore, if all the normal coordinates and velocities were zero, the LOS 
error would remain zero. 
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Fig. 5-la Line-of-sight error due to unit initial modal displacement in Modes 
1, 2, 3, 4, and 5, respectively. The separate time histories are plotted 
together for easy comparison and use with Measure 1. 
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Yig. 5-lb Line-of-sight error due to an initial modal displacement in Modes 6, 
7, 8, 9, and 10, respectively. The separate time histories are plotted 
together for easy comparison and use with Measure 1. 
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Pig. 5-2a Normal coordinates of Modes 1, 2, 3, 4, and 5, respectively. The 
cause of excitation was the rapid time-minimized bang-pause-bang pointing slew 
maneuver. These time histories are plotted together for easy comparison and 
use with Measure 2. 
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Fig. 5-2b Normal coordinates of Modes 6, 7, 8, 9, and 10, respectively. The 
cause of excitation was the rapid time-minimized bang-pause- bang pointing slew 
maneuver. These time histories are plotted together for easy comparison and 
use with Measure 2. 
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5.1.3 Measure 3: LOS Error Solely due to Each Mode Excited by the Man- 

euver 

To measure by the LOS error a vibration mode could cause, or to measure by the 
extent to which a vibration mode could be excited, seems to be a rather rea- 
sonable technique by itself, but the resulting rankings were inconsistent and 
rather confusing. For example, Mode 5 is the most important one by Measure 1 
but only the fifth by Measure 2. Moreover, Mode 5 could even be ignored 
because of its insignificant Measure-2 value (about two orders of magnitude 
smaller than the fourth). Similarly, Mode 1 ranks number 3 by Measure 2 but 
only number 6 by Measure 1. It was hard to determine rationally which modes 
would really need active control. A third measure was then developed. 

For this analysis, the LOS error caused by a single mode alone was calcu- 
lated separately, like for Measure 1, but the mode causing the error was 
excited by the very maneuver of concern, instead of initial conditions. All 
the initial conditions were assumed to be zero. On the other hand, the excita- 
tion of the vibration modes was exactly the same as for Measure 2, but the 
resulting LOS error, instead of the modal response, was taken as the measure. 

This measure is a sound rational combination of the cause (slew exci- 
tation) and the effect (LOS error) with respect to each vibration mode. 

It can appropriately indicate for each mode individually the extent to which a 
single mode could be excited, and the degree of LOS error this mode alone 
could cause if it alone were so excited and, hypothetically, no other modes 
were present at all. 

The 10 separate numerical results are shown together by the plots in Fig. 
5-3. Each curve represents the LOS error caused solely by a single mode 
while the mode was being excited by the rapid slew maneuver. The table below 
lists the highest peak value of the each curve. 

Mode: 123456789 10 

Peak: 3.26 88.6 9.57 6.53 0.33 .036 .077 .002 .004 .0002 

The relative importance of the 10 modes is thus given by the following des- 
cending order: 

Mode: 2, 3, 4, 1, 5, 7, 6,... 

An inspection of this ranking and the peak values will show that a signifi- 
cant break between the fourth- and fifth-ranked modes (i.e., modes 1 and 5, 
respectively). We thus selected the four top-ranked modes, i.e., modes 2, 3, 
1, and 4, as the primary candidates for modeled modes. 

Mode 5 is marginally important compared to other modes, but is the fifth in 
the rank and has a much higher value than the remainder. We therefore consi- 
dered it to be a secondary candidate for modeled modes. 

Mode 1 could have been ranked higher than Mode 4 if the time average were 
used instead. This would make no significant difference, however, since both 
were among the top four modes anyway, and these four had all been selected to 
be primary candidates for modeled modes. 
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Fig. 5-3 Line-of-sight error caused by a single mode alone while the mode was 
being excited by the pointing maneuver. The separate time histories are plot- 
ted together for easy comparison and use with Measure 3. 
b. Modes 6, 7, 8, 9, and 10. 









5.2 Analysis on Modal Control Influences of Actuators 

We recently discovered that the high-gain problem Canavin encountered in 
his first modal dashpot design [11] would not have existed if he had paid 
attention to the ill conditioning of the coefficient matrix ^^ T Bp he used in 
his numerical example. Among the 12 "controlled modes" he selected to form 
his reduced-order design model, several have very little control influences 
from the 32 actuators he used on the structure*. It is intuitively apparent 
that actuators having smaller influences on a given mode are less effective in 
controlling the mode, and thus require to be compensated with larger gains . 
Mathematically speaking, when the smallest "singular value" of the coefficient 
matrix one orc ^ er of magnitude smaller, the largest singular value 

of the resulting gain matrix G* as a solution of Eq. (4-7) generally is two 
orders of magnitude larger. This means that not all his "controlled* modes 
should be Included in the modal-dashpot design without any discrimination 
against excessively small control influences by the actuators. In other 
words, all the available actuators need not be lumped together to control 
all his "controlled* modes through one large feedback gain matrix. 

To make an effective design for the SCOLE configuration, we analyzed the 
control influences of the actuators first and match those modes in need of 
active control with the right actuators. 

For evaluating and comparing their modal influences properly, we grouped 
the actuators according to their location on the SCOLE configuration as well 
as their type. As a result, the actuators** were divided into the following 
four different groups: 

Group 1: Actuators 1 to 3, for applying moments on the Shuttle about its body 
x, y, and z axes, respectively; 

Group 2 : Actuators 4 to 6, for applying moments on the Reflector about its 
body x, y, and z axes, respectively; 

Group 3 : Actuators 7 and 8, for applying forces at the Reflector mass center 
in the x and y directions, respectively; 

Group 4 : Actuators 9 to 12, for applying forces at two specific points on the 
— Mast beam in the x and y directions, respectively. 

The control influences on each mode, say mode i, from all the actuators in 
a specific group can be summarized by calculating their RMS (Root-Mean-Square) 
value 

^(♦Jb pi )* +(*?b F2 ) 2 ... +(*Ib pk )* J/k 

over the group. Listed in Table 3-1 are these RMS values in the descending 
order. 


* The antenna-like structure consisted of a large dish in the forward section 
and a gimbaled equipment section to the aft. It had 32 member dampers (as the 
co-located actuators and rate sensors). Its finite-element model has 35 deg- 
rees of freedom. 
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Table 5-1 RMS Actuator Influences on first 10 Modes 


Group L 

Mode Act. 1 - 3 


2 0.30019961E-02 

4 0.41220308E-03 

1 0.40146321E-03 

3 0.19184369E-03 

5 0.11186474E-03 

6 0.69881789E-04 

7 0.36829457E-04 

8 0.26261532E-04 

9 0.15072107E-04 

10 0.13497747E-04 


Group 2 

Mode Act. 4 - 6 


5 0.36487188E-01 
4 0.25172627E-01 
3 0.15999462E-01 
2 0.15595800E-01 

7 0.14711439E-01 
1 0.13037169E-01 
9 0.57048416E-02 

6 0.34139471E-02 

8 0.12352261E-02 
10 0.63637015E-03 


Group 3 

Mode Act. 7-8 


2 0.14311218E+01 

1 0.14061384E+01 

3 0.81986851E+00 

4 0.39743480E+00 

7 0.30395976E+00 

9 0.25503686E+00 

6 0.21852742E+00 

8 0.14623879E+00 

10 0. 10801 539E+00 

5 0.74399590E-01 


Group 4 

Mode Act. 9-12 


2 0.10711402E+00 

1 0.10338868E+00 

3 0.10171293E+00 

4 0.69439910E-01 

9 0.68373762E-01 

8 0.67025743E-01 

5 0.63191518E-01 

10 0.46103600E-01 

6 0.39935779E-01 

7 0.32263912E-01 


Table 5-2 RMS Sensor Influences on first 10 Modes 


Mode 

Group 1 

Sen. 1-3 Mode 

Group 2 
Sen. 4-6 

Mode 

Group 3 
Sen. 7-8 

Group 4 

Mode Sen. 9-12 

2 

0.28690067E-03 

5 

0.34890966E-02 

2 

0.13466856E+00 

2 

0.10711402E+00 

4 

0.39390128E-04 

3 

0.32387748E-02 

1 

0.12736945E+00 

1 

0.10338868E+00 

1 

0.39113598E-04 

4 

0.24055073E-02 

3 

0.12407852E+00 

3 

0.10171293E+00 

3 

0.18940789E-04 

2 

0.15346858E-02 

4 

0.38158901E-01 

4 

0.69439910E-01 

5 

0.10689848E-04 

7 

0.14879148E-02 

7 

0.36793593E-01 

9 

0.68373762E-01 

6 

0.66779044E-05 

1 

0.13531352E-02 

9 

0.30879460E-01 

8 

0.67025743E-01 

7 

0.35194691E-05 

9 

0.67911280E-03 

6 

0. 2083207 5E-01 

5 

0.63191518E-01 

8 

0.25095521E-05 

6 

0.32624099E-03 

8 

0.14005536E-01 

10 

0.46103600E-01 

9 

0.14402938E-05 

8 

0.11804001E-03 

10 

0.10299906E-01 

6 

0.39935779E-01 

10 

0.12898448E-05 

10 

0.60812166E-04 

5 

0.90737212E-02 

7 

0.32263912E-01 


Observe that Actuators 1 to 3 (Group 1) have an RMS value for Mode 2 that 
is one order of magnitude higher than all other modes, and hence are most 
effective in controlling Mode 2 than controlling other modes. Observe also 
that Mode 2 ranked the highest in RMS value with respect to Group-3 actuators 

7 and 8. In addition, this RMS value is two orders of magnitude higher than 
that with Actuators 1 to 3. Consequently, Actuators 7 and 8 should be more 
effective for controlling mode 2 and require much smaller feedback gains. 
Note that Actuators 9 to 12 (Group 4) are less effective than Actuators 7 and 

8 in controlling Mode 2. 

With a similar argument, Actuators 7 and 8 are also most effective in con 
trolling Mode 1. Therefore, Modes 1 and 2 and no more others should be 
selected as the ''modeled modes" in the design of the modal dashpots using 
Actuators 7 and 8. 

Since Mode 3 is a torsion mode and is more appropriate to be controlled by 
moments than forces. The RMS values clearly suggest that Actuators 4 to 6 
(Group 3) will be more effective than Actuators 1 to 3 for controlling Mode 3. 
Although among the four groups, Actuators 9 to 12 did have the highest RMS 
values of control influences on Mode 3, we did not expect the proof-mass actu- 
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ators (9 to 12) to be capable of suppressing large torsional vibrations of any 
mode in such a very short time. 

The RMS control influence on Mode 4 is larger from Actuators 4 to 6 than 
from Actuators 1 to 3. Consequently , both Modes 3 and 4 should be selected 
as the ’’modeled modes'* in the design of modal dashpots using Actuators 4 to 
6 . 


All the four primary condidates have been selected as the ’’modeled modes" 
for the appropriate matching groups of actuators. A review of Table 5-1 will 
show that Mode 3, the secondary candidate, has a higher RMS value of control 
influences from the same Group 2 of actuators than both Modes 3 and 4. 
According to a previous study by Lin and Jasper [24], such a situation would 
result in large spillover of Mode 5, severe dynamic interactions of modeled 
modes with unmodeled modes, and significant degradation of damping performance 
if Mode 5 were not also modeled with Modes 3 and 4 for control by Actua- 
tors 4 to 6. 

To summarize, this analysis shows applying moments and forces at the Ref- 
lector end of the Mast beam will be more effective in controlling the excited 
vibrations" in SCOLE (and particularly Modes 1 to 5) than at the Shuttle end or 
at the intermediate points of the flexible mast. Instead of lumping up all 
candidate modes (1 to 3) to be controlled by Actuators 4 to 8 together, the 
designer for modal dashpots should match these modes with their most effec- 
tive or most appropriate actuators. Specifically, Modes 1 and 2 should be 
controlled by Actuators 7 and 8 and Modes 3 to 5 by Actuators 4 to 6 . 

There is no need to include more modes to each group since there are 
enough actuators to be distributed among all the 5 most important modes of 
the SCOLE flexible-body dynamics. Including more modes may not always help: 
it might simply Increase the magnitude of the feedback gains without any real 
benefit, particularly when the additional modes are of significantly smaller 
control influences by the actuators; the increased feedback gains might in- 
stead amplify various adverse effects of control spillover and system noises. 


5.3 Design of Modal Dashpot MD1 

The modal dashpot MD1 was designed for SCOLE for quick suppression of the 
excessive vibrations excited by the rapid BPB LOS pointing slew maneuver. It 
is composed of two parts. Part 1 is for applying forces at the Reflector mass 
center in the two transverse directions using a feedback of linear velocities 
at the Reflector end of the beam. Part 2 is for applying moments also at the 
Reflector about the three body axes but using a feedback of angular velocities 
instead. 

The location of these actuators are the same as specified by Taylor in 
Ref. 1 for the control forces and moments at the Reflector. The sensors were 
located where the "outputs" of Dr. Joshi's modal data set D3D585 had been cal- 
culated. Some of the control inputs (u^) and observation outputs (yj) were re- 
labeled for technical convenience. Sensors 1 to 8 are not really co-located 
with the correspopnding actuators, but note that their RMS values of modal 
observation influences (Table 5—2) exhibit virtually the same patterns as 
those of modal control influences (Table 5-1). 
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5.3.1 Part 1: Linear Velocity Feedback Force Control 

Two force actuators (or equivalently, a single force actuator capable of 
delivering separate forces in two independent axes) are assumed to be placed 
at the center of the Reflector. The force inputs* U 7 and ug (in the x and y 
directions, respectively) are each limited to 800 lb as specified. Two linear 
velocity sensors (or equivalently a single velocity sensors capable of measur- 
ing the rate of linear displacements in two independent axes) are assumed to 
be located at the Reflector end of the mast beam. The sensor outputs** yi 5 
and yi 6 (in the x and y directions, respectively) represent the time rate of 
deflection of the Mast beam at the Reflector end relative to the Shuttle end. 
Note that these sensors are only approximately co— located with actuators: 
they are apart by 18.75 ft and 32.5 ft in x and y directions, respectively, 
whereas the beam is 130 ft long. 

The design problem is thus to determine a 2x2 gain matrix Glvr for the 
following linear velocity feedback control law 


’ u ?' 

= -g 

’*15" 

_ u 8_ 

LVR 

_ y 16_ 


(5-1) 


The foregoing analysis of the control influences has suggested that only 
Modes 1 and 2 be selected as the "modeled modes” for this part of design. 
Accordingly, the control and observation influence matrices *M T Bp and Cy#M 
on the two modeled modes to be used in the modal dashpot design equation (4-8) 
have the following numerical values: 


To [ .19875923E+01 
*M B F " I .14599262E+00 

c f .18012760E+00 

V*M “ 1_ .55188192E-04 


.62669927E-0l“l 

-.20186396E+01J 

.21140305E-0ll 

-.18927317E+OOJ 


(5-2a) 


(5-2b) 


Before solving the corresponding design equation (4-8) for a specific gain 
matrix, we must specify the desired value for the additional damping matrix 
A|q. For technical simplicity, we choose it to be diagonal, so that its diago- 
nal elements SjJl and ®M 2 can Be use< i rather directly for guiding the modal- 
dashpot design. Since both modes 1 and 2 substantially dominate the vibratory 
response of the SCOLE configuration to the BPB point|yg maneuver, we wish to 
augment each with active damping as close to 70.7% of critical damping as 


* These correspond to U 4 * F rx and U 5 = F ry , respectively, in Dr. Joshi's nota- 
tion. 

** These are indirectly equal to the derivatives of the deflections yy m £ x anc * 
yg = £y in Dr. Joshi's notation. 

+ These terms represent the additional damping coefficients in the correspond- 
ing decoupled equations of motion; see Eqs. (4-14)-(4-15). In multivariable 
. oot-locus analysis, these values also represent the "rate of departure from 
the open-loop poles when the feedback loops are closed. 

tt 70.7% is an optimal value in the sense that the second-order system cor- 
responding to the single mode will neither be too sluggish nor have a large 
overshoot . 
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possible. Let us attempt a theoretical 2% settling time of 3 seconds for 
Mode 2 in estimating the desirable additional damping coefficient $ 2 - The cor- 
responding time constant is 3/4 sec; thus by definition 

^ 2^2 = 4/3 

where q denotes the closed-loop damping ratio desirable of Mode 2. Substi- 
tuting in the natural frequency w 2 = 1-97024 rad/sec yields 

l 2 = 0.6767, 

which is acceptably close to the optimal value. Whence, the closed-loop damp- 
ing coefficient desirable of Mode 2 is given by 

6 2 s 2 4 w 2 = 8/3 ’ 

Since an inherent damping of 0.3% has been specified in the data set 
M D3D585 M for each mode, the desirable additional damping coefficient^ desirable 
of Mode 2 is 

6 * = *2 " 2 C 2 w 2 = 8/3 " 2 x °' 003 x 1 * 97024 " 2.6548 

Next, we choose the additional damping desirable of Mode 1 to be 60%, 
i.e., ti* = 0.6, since Mode 1 has a smaller magnitude of vibration than Mode 2. 


In summary, the desired damping coefficients as in Eqs. (4-14)-(4-15) for 
the two modeled modes are then readily given as 

<5* = 6* = 2 c* u = 2 x 0.6 x 1.7470 = 2.0964 (5-3a) 

6*2 = = 2 -6548. (5-3b) 

Now the feedback gain matrix is readily obtained from solving (4-8) as 


g lvr 


.58420630E+01 .43392044E+00“ 

.42038249E+00 .69796355E+01_ 


(5-4) 


3.3.2 Part 2: Angular Velocity Feedback Moment Control 

Three torquers (or equivalently, a single torquer capable of delivering 
separate touques about three independent axes) are assumed to be located on 
the Reflector. The torque inputs u^, U 5 , and ug (about the x, y, and z axes, 
respectively) are each limited to 10,000 lb-ft as specified. Three angular 
velocity sensors (or equivalently a single sensor capable of measuring separ- 
ately the rate of rotations about three different axes) are assumed to be 
located at the Reflector end. The sensor outputs** y^Q, yil> and y^ ( a ^ out 
the x, y, and z axes, respectively) represent the time rate of rotations of the 


// The corresponding additional damping ratio C2 f° r M°d e 2 is 0.6737. A few 
slightly modified values were also tried when the design of this part was rep- 
eated (see Section 6.2). 

* These correspond to Dr. Joshi's u£ = T rx , U 7 = T ry , and ug * T rz , respec- 
tively. 

** These correspond to Dr. Joshi's y^ = $r> yi3 = ® r , and yi 4 * i r , respec- 
tively. 
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Reflector end of the mast due to bending and torsion. Again, these sensors 
are only approximately co— located with the actuators. Note that the Ref- 
lector itself is a rigid body. 


The design problem is to find a 3x3 gain matrix G AVK for the 
angular velocity feedback control law 


" u 4 " 


y 10 

u 5 

“ “ g AVR 

y ll 

_ u 6_ 


_ y 12_ 


following 


(5-5) 


Similarly, as suggested by the previous analysis on control influences, we 
choose Modes 3, 4, and 5 to be the "modeled modes" for this part of the 
design. Accordingly, the specific control and observation influence matrices 
are given by 


M F 


C * 

V M 


•10761780E-01 

•37541741E-01 

.31739483E-01 

.10592386E-02 

•17367745E-02 

•52450669E-02 


.18101653E-01 

.22172032E-01 

.54643290E-01 

•35194610E-02 

.21135667E-02 

.13325330E-04 


•18012847E-01 

.45802862E-04 

•81303515E-03 

•30149737E-02 

.52175940E-02 

•23641118E-03 


(5-6a) 


(5-6b) 


Since the vibratory responses of Modes 3, 4, and 5 are much smaller in 
magnitude than those of Modes 1 and 2, it is reasonable to augment them with 
only a relatively small amount of active damping. We chose rather arbitrarily 
3% of critical damping for each. The diagonal elements of the desired addi- 
tional damping matrix are then given as follows: 


* 

Ml 

= 6 3 

» 2 

X 

0.03 

x u» 3 

= 0.3065 

(5-7 a) 

CM 

* S 

-< 

=* 2 

X 

0.03 

X 0)/ 

4 

= 0.4470 

(5-7b) 

* 

M3 

- s 

= 2 

X 

0.03 

X (i). 

4 

= 0.7742 

(5-7 c) 


Substituting (5-5)-(5-7) in Eq. (4-8) and solving the resulting equation, we 
get the following gain matrix 


G AVR 


.24172707E+04 

•15734103E+03 

.13433660E+04 


.16653096E+03 

.21781213E+04 

-.22055215E+04 


•45158162E+03 

-.72768193E+03 

.42951681E+04 


(5-8) 



6. PERFORMANCE OF VIBRATION CONTROL DESIGNS ON SCOLE 

To evaluate the vibration control performance of the modal dashpot design 
MD1, we incorporate the two feedback laws (5—1) and (5-5) into the same SCOLE 
flexible-body dynamic model as was used in simulating its vibratory responses 
to the BPB pointing maneuver. As stated in the beginning of Section 3, the 
undeformed SCOLE configuration was assumed to have been aligned with the spec- 
ific attitude of zero LOS error. Thus, the velocity-sensor outputs would con- 
tain only the flexible— body rates*, just as desired for feedback control of the 
excited vibrations. Recall that the model is of the "full order" in the sense 
that it includes all the ten modes as provided in the data set D3D585. In 
order that the control moments and forces do not exceed their specified 
limits, the computer program also simulates the saturation of the actuators at 
their respective limits. For example, if at any time the feedback control 
input, say, ug would command the actuator to exert more than 800 lb force to 
the Reflector, the actual force applied would be only 800 lb maximum. 

The feedback control consisting of the two parts of modal dashpot MD1 is 
turned on right after the completion of the BPB pointing maneuver. Thus the 
terminal state of the SCOLE vibrations (i.e., the LOS error, the deflection and 
its rate of change, the angular displacement and its rate, modal displacements 
and velocities,...) at the end of the maneuver become the initial conditions 
of the feedback controlled system. The vibration control is applied for five 
seconds, which is about the same duration as of the pointing maneuver. We in- 
tentionally use such a rather "long" period in order to check if instability in 
the closed-loop system might start to develop after the excessive vibrations 
has been rapidly forcefully suppressed. Various versions of the modal dashpot 
design MD1 (each with a slightly different value for the additional damping 
coefficient fi 2 *) were evaluated. Reported below are two representative cases. 


6.1 Simulation Results of Modal Dashpot Design MD1 

The specific values of the gain matrices G LVR given by (5-4) and G AV r given 
by (5-8) were incorporated with the control laws (5-1) and (5-5) respectively 
in the full-order dynamic simulation. The simulation results are summarized 
by time-history plots in Fig. 6-1. 

The history of the applied moments and forces (Fig. 6-la) shows that the 
applied moment about each axis never exceeded the limit of 10,000 lb-ft, nor 
did the applied force in each direction exceed the limit of 800 lb. Large 

moments and forces were needed only during the early portion of the control 

period, but did not exceed the limits because of "saturation”. All the applied 

forces and moments quickly reduced to minimum automatically because the 
sensed rates of vibrations rapidly became insignificant. 

Fig. 6-lb shows that the LOS error was rapidly subdued to 11.79° from 

71.43° where the pointing maneuver ended. Note that the initial LOS error con- 
tinued to rise to 85.29° (or 115.13“ = 180° - 64.87° if not taking the principal 


* If the configuration had not been so aligned, then rigid-body rate would also 
be present and some filtering or signal processing might be required. Alterna- 
tively, one could use relative sensors instead of inertially referenced sen 

sors. 
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Fig. 6-1 Simulation results of vibration control design MD1; 

a. Histories of applied moments and forces at the Reflector. 
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Fig. 6-1 Simulation results of vibration control design MD1; 

b. Histories of line-of-sight error and Mast tip deflection. 
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value of the arcsine function) because of its large rate of change at the very 
instant of switching from pointing maneuver to vibration control. Neverthe- 
less, the LOS error was suppressed down within 18.46* after only 2.5 sec of 
active vibration control and remained below 17.54* thereafter. Moreover, it 

was even reduced to 11.79* in 3.1 sec. 

The bending of the mast beam was very rapidly suppressed to virtually 
null. Observe in particular that even the y-directionally deflection at the 
Reflector end (i.e., yi4) continued to increase to 119.7 ft (again because of 
the large "initial" rate of change), it was suppressed down into the band of 
±5 ft in only 2 sec, and into the band of ±0.5 ft in 4.2 sec. It is interest- 
ing to notice that it took less than 2.9 sec* to settle within 2.4 ft, 2% of 
the peak value. Recall that a 2% settling time of 3 sec for Mode 2 was used 
in the design. 

Fig. 6-lc shows the rapid reduction of the initially large deviations in 
the Shuttle and Reflector attitude angles to zero in a very short time. The 
last peak diviation of the Reflector roll, pitch and yaw attitude angles is 
only 0.460*, 0.546* and 1.360*, respectively. 

Fig. 6-ld shows that the large-magnitude vibrations of first five modes 
all were rapidly suppressed to virtually zero in a very short time. Note in 
particular that this vibration control was very effective for quick reduction 
of the excessively large magnitude of Mode 2. Observe, on the other hand, 
that Mode 5, the secondary candidate, was reduced only in a moderate rate by a 
moderate amount, but recall that it was not really significant at first place 
with respect to excitation by the BPB maneuver nor its contribution to the LOS 
error. Mode 5 did not need much active control anyway, and hence only a very 
small additional damping was designed for it. 

The vibrations in other modes (i.e., Modes 6 to 10) remained virtually in 
the same insignificant levels as before, and hence their plots are omitted. 
Still their magnitudes were more or less decreased with time because of some 
concomitant additional damping as a side benefit of spillover.** 


6.2 Simulation Results of A Modified Version of MD1A 

We also tried a few other versions of the modal dashpot design MD1 by var- 
ying the additional damping coefficient 6 2 * desired of Mode 2. Because of the 
saturation of the actuators at the imposed limits, it is reasonable to consider 
some smaller feedback gains. The following is a typical case. 


* The peak occurred at t * 0.5 sec whereas the deflection was -2.24 ft at t * 
3.4 sec. 

** In the standard LQG design, one will generally try hard to reduce spillover 
because it has been well known to degrade performance and even to introduce 
closed-loop instability. With a modal-dashpot design, spillover can be bene- 
ficial instead, in leaking some active damping forces to unmodeled modes. 
Such is particularly the case when the design is not carefully focussed. The 
side benefit in our design was intentionally minimized because we tried to 
maximize our effort on the most important modes and matched them with most 
influential control actuators to minimize the leak. 
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This version, let us call it MD1A, is almost the same as before, except 
that the additional damping ratio c 2 * desired of Mode 2 was arbitrarily set 
equal to that of Mode 1; namely, 

C* - C* = 0.60. 

Therefore, 

6* = 2 5* (o 2 = 2 x °.6 x 1.97024 = 2.36429 (6-1) 

Mode 2 would then have a theoretical 2% -settling time of 3.38 sec. 


Using the new value for 6 2 * in (5-3b), and repeating Part 1 of the modal 
dashpot design, the following new value of the feedback gain matrix was read- 
ily obtained. 

f .58420557E+01 .45784262E+00“| (6 _ 2) 
G LVR “ L •42061494E+00 .62209375E+01J 

The same simulation was then repeated with these new values. Results are 
summarized by the plots in Fig. 6—2. The applied moments and forces shown in 
Fig. 6-2a are virtually the same as before (Fig. 6-la) with only some invisible 
differences. Some meaningful differences do exist in the histories of LOS 
error and beam deflection. 


Observe that the LOS error (in Fig. 6-2b) quickly reduced to about 9.57* 
from the same initial value (71.43°). Similarly, due to large initial rate of 
change at the end of the pointing maneuver, the LOS error also continued to 
rise to 85.61° (or 180° - 64.81° * 115.19° if not taking the principal value of 
the arcsine). The large LOS error was suppressed down to the level of 
16.66* in 1.8 sec, and remained under it thereafter. Moreover, it was reduced 
to 9.57° also 3.1 sec after the vibration control began. 

The bending of the Mast was also suppressed down very rapidly. Though it 
continued to increase to 120.28 ft, the y-directionally deflection at the Ref- 
lector end (i.e., yi4) was suppressed down into the band of 17.35 ft in less 
than 1.8 sec., and into the band of ±0.75 ft in 3.7 sec. It took less than^3 
sec for the large deflection to settle down to the band of 2 % of the peak , 
i.e., ±2.4 ft. Recall that 3.38 sec is the theoretical 2%-settling time used 
for Mode 2 in this modified design. 

The histories of attitude changes (Fig. 6-2c) and modal responses (Fig. 
6-2d) are again virtually the same as before (compared to Figs. 6-lc and 6-ld, 
respectively) with only some invisible or insignificant differences. The last 
peak deviation of the Reflector roll, pitch, and yaw attitude angles is only 
0.714*. 0.582*, and 1.399*, respectively. 


* The peak occurred at t * 0.5 sec., whereas the deflection was only -1.95 ft 
at t * 3.5 sec. 


229 






ERROR 



O 


o 

on 



TIME SEC 


F0010 


Fig. 6-2 Simulation results of vibration control design MD1A; 

b. Histories of line-of-sight error and Mast tip deflection 





Fig. 6-2 Simulation results of vibration control design MD1A; 

c. Histories of attitude deviations at Shuttle (S) and Reflector (R) ends 
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6.3 Comments 


6.3.1 The modal dashpot designs of vibration control met the vibration 
control design challenges fairly well and are effective and fast in sup- 
pressing excessive vibrations. Excited by BPB type rapid pointing slew man- 
euver, the flexible mast beam deflected between +114 ft and -113 ft, but such 
an excessive vibration was then quickly suppressed down to less than 0.75 ft 
in less than 3.7 sec after a modal-dashpot vibration control was turned on. 
The roll vibration of the Reflector between -86.96° and +88.35° during the man- 
euver was also quickly suppressed down to less than 0.72°. The large LOS 

error of 89.8° (or 133.3° if not taking the principal value of the arcsine) was 
also reduced quickly to less than 17.54°. 

6*3.2 The original version of the modal dashpot design MDl performed 

slightly better than the modified version MD1A in suppressing the deflection of 
the Mast beam and all the attitude deviations, but not so well in reducing the 
LOS error. The modified version used a slightly smaller additional damping 
ratio for Mode 2 in the design, i.e., 

^ 2 * 13 0.60 instead of “ 0.6737. 

6.3.3 When a velocity feedback control, whether it is of the modal-dash- 
pot type or not, is not properly designed, even feedback gains of an intermedi- 
ate magnitude can cause severe interactions between modeled and unmodeled (or, 
equivelent ly , between "controlled” and "uncontrolled") modes, and hence badly 
degrade the desired performance of active damping augmentation [24]. The res- 
ults of these two versions have shown, on the other hand, that if modal dash- 
pots are properly designed, both the modal interactions and the performance 
degradation are not problems. 

Thus, some of the additional damping can be as high as the optimal value 0.707 
if necessary, hence can have high feedback gains , to be really effective in 
quick suppression of vibrations. In other words, not all velocity output 
feedback vibration controllers are of low authority, low performance ! 

6.3.4 Now, not having to worry about the spillover and modal interaction 
problems, the feedback gains of properly designed modal dashpots ideally can 
be as high as the designer wishes. High gains can be as desirable for flexi- 
ble-body vibration control as they have traditionally been for effective con- 
trol of rigid bodies. 

High gains are desirable for generating comparable negative feedback to offset 
the vibrations. Theoretically, the higher the better. For example, the ver- 
sion MDl has a higher gain (because of higher t 2 *) than the version MD1A, the 
deflection and attitude deviations can be continuously suppressed down to 
smaller values (e.g., ±0.5 ft vs ±0.75 ft, in deflection; 0.46° vs 0.714° in 
Reflector roll angle,...). 

The size of the feedback gains for a properly designed modal-dashpot 
vibration control is virtually limited only by the force and torque capability 
of the actuators. Since the vibrations were initially very large, the high 

gains resulted in requiring larger forces and torques than their limits. The 
simulated saturation thus restrict the applied force/ torques to the limits. 
Therefore, there are no needs to be concerned with high gains as much as 
before, even the actuators may saturate at their force/torque limits. 
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6.3.5 Recall that the badly excited modes, i.e., Modes 2 and 1, were made 
most strongly controllable and obserable by carefully matching them with 
actuators and sensors with the strongest influences. Recall also they were 
strongly controlled by selecting them as modeled modes In the Part 1 design of 
the modal dushpots and by adding to them the highest additional damping 
ratios. The results of Sections 6*1 and 6.2 show that the resulting modal- 
dashpot designs are very effective for fast suppression of large vibra- 
tions. The results also show that spillover is minimum and that the unmo- 
deled modes receive some small concomitant additional damping because of spil- 
lover. 

6.3.6 Unlike all other vibrations (deflections, attitude deviations and 
modal responses), the LOS error was not reduced to a smaller value by the 
version MD1 than by the version MD1A. Also, the LOS error was not continu- 
ously reduced to near zero as were all other vibrations, although the reduc- 
tion from its excessively high peak was quite substantial. Nonlinear proper- 
ties of large Euler attitude angles, and the truncated forces and moments 
from the saturated actuators (due to high feedback gains) likely are the 
causes. We have no clear explanations at the present time. Nevertheless, 
observe Figs. 3-lc, 6-lc and 6-2c that the Reflector continued to have suffi- 
ciently large pitch and yaw ratations during the initial phase of the vibration 
control, in addition to the main (and larger) roll rotations. 

6.3.7 Figs. 6-1 and 6-2 show that after the excessive vibrations have all 
been suppressed down to sufficiently low levels, the time rates of change 
naturally start to become much less significant, and the modal-dashpot 
vibration control also starts to become less effecitve. Unless the feedback 
gains are increased thereafter, the vibrations may not continue to be reduced 
to the desired precision in a reasonably short time. One way to achieve the 
desired precision is to start to increase the modal-dashpot gains progressively 
after the vibrations become sufficiently small, e.g., after 2 seconds of the 
initial vibration control. 

Another way is to switch to some form of "modern control" for complet- 
ing the vibration suppression and precision pointing. When all the displace- 
ments and rates of change have become reasonably small, the whole dynamic 
system becomes legitimatel y linear , and the LOS error expression legitimately 
linearizable. The condition is very suitable for application of the modern 
optimal state-feedback control technique. 

Modern control using standard Linear-Quadratic-Gaussian (LQG) optimal state 
regulators and optimal state estimators has traditionally performed very well 
in precision pointing and attitude control of rigid— body systems, even using 
small signals. For application to a flexible-body system, the modern control 
must be very carefully designed, however; otherwise the notorious spillover 
problems may destabilize the system instead! 

6.3.8 Several major approaches to extend or adapt the LQG design tech- 

niques were proposed during the years of ACOSS (Active Control of Space Struc- 
tures) and VCOSS (Vibration Control of Space Structures) programs [38]-[45], 
[15], [46]-[50]. Either the weighting matrices in the control performance 

index is modified in some ways [51]- [52], or some positivity requirement is 
imposed on the design [53], or some pre-design compensation of the actua- 
tor/sensor influences is made [35]- [36]. All were successful to some limited 
extents in addressing the major challenge of spillover problems, but are not 
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readily applicable to realistic large flexible space structures. 

Formal applications of the robustness theory [54]-[55] were started rec- 
ently [56]- [59]. The method of loop transfer recovery (LTR) was also 
applied to recover sizable gain and phase margins of LQ regulators. The modi- 
fication recently proposed by Blellock and Mingori [58] appears to have made 
the LTR method more directly applicable to LQG controllers designed for large 
space structures, so far as the uncertainties in the modal frequencies of the 
plant are concerned. Recent results obtained by Sundararajan, Joshi, and Arm- 
strong [59] are rather encouraging. Based on their interpretation of spillover 
problems as additive uncertainty [60], [55], they were able to make an innova- 
tive application of the LTR method to overcome spillover problem with their 
LQG attitude controllers designed for the Hoop/Column antenna. This approach 
has a great potential for practical application to realistic large flexible 
space structures, since it appears to be able to overcome the spillover prob- 
lem of an unlimited number of unmodeled modes . 

Incorporating modal dashpots into a LQG or LQG/LTR design and following a 
similar sequence of careful pre-design analyses certainly will greatly enhance 
the stability and performance of the resulting LQG/MD or LQG/LTR/MD vibra- 
tion controller. The two proof-mass actuators placed on the mast beam may 
be used together with all the force and moment actuators on the Reflector and 
the Shuttle for such a low-power but high-precision control. 
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7. CONCLUSIONS 


7.1 The two-stage approach is a feasible and promising one for rapid slewing 
and precis iorT pointing/ retarge ting - ^ f larg e fl exible spac e syst ems and, in par- 
Qc uTar", ~t he orbital SCOLE c onfiguration . It is capable of rapidly slewing 
the line-of-sight and settling the excited vibrations in a minimum time. 
The resulting control design, in general, will consist of the following three 
parts in cascade: 

Stage 1: a bang-bang type rapid slew maneuver based on the rigid-body dynamics 
for pointing/ retargeting in a minimum time; if excessive vibrations may 
be excited, using smaller forces and moments should be considered. 

Stage 2, Part 1: a high-power modal-dashpot design of velocity output feedback 
control based on the flexible-body dynamics for fast and effective 
reduction of large excited vibrations to a small magnitude; 

Stage 2, Part 2: a LQG/LTR design of optimal state feedback control augmented 
with a broad-band low-power modal-dashpot design of velocity output 
feedback control, also based on flexible-body dynamics, for (i) achiev- 
ing the specified pointing accuracy in a short time and (ii) maintaining 
the precision and closed— loop system stability. The LQG/LTR design 
may be incorporated or integrated with an appropriate modal-dashpot 
design. 


7.2 Not a ll bang-bang (BB) type of time-minimized slew maneuvers will excite 
large structural vibrations. When large forces are used up to their extremes 
(for example, 800 lb on the Reflector) to complete the specified slew angle 
(20°) of the rigidized configuration in the shortest time, the excited vibra- 
tions can be excessively large in magnitude (e.g., a 114-ft peak deflection of 
the 130-ft Mast beam), even only moderated maneuvers of the bang-pause-bang 
(BPB) type is used instead. On the other hand, when properly selected small 
forces, e.g., 25 lb, of the kind of vernier RCS thrusters onboard the Space 
Shuttle, are used, even BB-type maneuvers will excite very little vibra- 
tions (e.g., 0.3 ft peak deflection of the Mast beam). 


If the excited vibrations are excessive, a "high-power modal-dashpot design 
of velocity output feedback control can be used in the first part of the Stage 
2 to suppress the vibration down to a reasonable small magnitude quickly and 
effectively. If the excited vibrations are relatively small, or have already 
been suppressed to a small magnitude, some modified form of linear-qua- 
dratic (LQ) optimal state feedback control augmented with a "low-power" 
design of modal dashpots can be used in the Stage 2 to achieve the desired 
pointing precision. 


7^2.1 The vibration modes of the SCOLE configuration were excessively 
excited when an 800-lb force was applied on the Reflector in the y direc- 
tion during a BB type slew maneuver. When the best Stage-1 design, i.e., 
the BPB roll-axis slew having the best LOS pointing accuracy with a minim- 
ized slew time (4.89 sec) and the least sensitivity to nonzero products of 
interia) was applied to the SCOLE flexible-body dynamics, the Reflector end 
of the mast vibrated between +114 ft and —113 ft, the Reflector rolled 


237 



between -86.96° and +88.35°, and the line of sight jittered between 89.8° 
(or 133.3° if not taking on the principal value of the sine function) and 
14.7°. 


Our carefully designed modal— dashpot type of velocity output feedback 
control was able to suppress the excessive vibratiions quickly and 
effectively: the Reflector end deflection down to ±5 ft in 2 sec , and to 
±0.5 ft in 4.2 sec; Reflector roll to ±3.48* in 7.1 sec, and to ±0.54° in 4.4 
sec; and the LOS error down to 11.79° in 3.1 sec. In other words, after 
only about 2 to 3 seconds of applying the ”high-power" modal dashpots, the 
vibrations were reduced to a region where some form of linear-quadratic 
optimal ’’state" feedback control (properly augmented with "low-power" 
modal dashpots) would be effective in futher reducing the vibrations and 
LOS errors to the desired precision. 


7.2.2 The large magnitude of the force, i.e., 800 lb, applied on the Ref- 
lector was responsible for the excessive excitation of vibrations in 
the 3C0LE configuration. Whether bang-bang type time-optimal slew man- 
euvers would excite excessive vibrations or not depends on the allowable 
maximum magnitude of the applied forces. When the limit of the force was 
decreased to only one tenth (i.e., 80 lb) but the pointing slew maneuver 
was still performed in a similar time-optimal bang-bang manner for the 
same 20° angle, the excited vibrations were significantly decreased. The 
maximum LOS error was 24.7°, comparable to the specified initial value (20°) 
due to the initial misaligment of the SC0LE configuration. The maximum tip 
deflection (20.6 ft) of the mast beam was also quite reasonable compared 
to the length of the Hast (130 ft). When no additional forces were 
applied, however, the vibrations excited by the applied moments alone in 
creased, instead* 


We found that if the applied force on the Reflector was about 25 lb , 

i.e., in the range of the vernier RCS thrusters used on the Space Shut- 

tle, the corresponding time-minimized bang— bang pointing slew maneuver 
would excited very little vibrations in the SC0LE configuration. The 
Reflector end of the mast vibrated only between +0.25 ft and -0.30 ft , the 

Rellector rolled only between +0.16® and -0.30°, and the LOS error was aj: 

most 0.51°. If the BB slew maneuver was followed Immediately by some 
form of linear-quadratic optimal "state” feedback control (properly aug- 
mented with "low-power" modal dashpots), such small vibrations and LOS 
errors would be easily reduced to the desired precision. 


7.2.3 During Stage 1, the BB maneuver using a 25 lb force on the Reflector 
required 10.96 seconds to complete the 20° slew while the BPB amneuver 
using a 800 lb force on the Reflector required only 4.89 seconds. A 
"high-power” modal-dashpot design of velocity output feedback control 
required additional 2.5 to 3 seconds to bring the excessive vibrations 
excited by the 800-lb maneuver down to the same order of magnitude as the 
vibrations excited by the 25-lb maneuver. Therefore, the total time 
required for both Stage 1 (slew) and Stage 2 (stabilization and pre- 
cision pointing) is likely to be around 10 and 12 seconds , respectively, 
for the two cases. 
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The two stages of the 800-lb BPB maneuver will probably require the least 
total time, but the excessive vibrations during the maneuver are impracti- 
cal and undesirable. 


The 80-lb BB maneuver requires a similar high-power design of modal- 
dashpots for quick and effective suppression of the moderately large vibra- 
tions to the same order of magnitude as the case of 25-lb maneuver. The 
total time required for the two stages is likely to be also around 12 
seconds or a little less. 


7.3 Although modal-dashpot type of velocity output feedback control can be 
designed as a usual diffuse (or ”broad-band M ) low— power (or "low-authority ) 
control, the simulation results of our careful designs have shown that modal 
dashpots can als o be a concentrated high-power ("high-authority") control for 
fast and effective suppression of large vibrations. Careful pre-design ana- 
lyses made it possible to do so for SCOLE. 


7*3.1 Our pre-design analysis on the vibration modes of the SCOLE configu- 
ration shows that modes 2,3, 4, 1,5 are the five most important modes 
requiring for vibration control and LOS error reduction, with mode 2 need- 
ing active control the most. 


7.3.2 Our Pre-design analysis on the modal control influences of the actu- 
ators shows that: two force actuators on the Reflector in x and y 

directions, respectively, are most effective for controlling modes 1 and 2; 
three moment actuators also on the Reflector about the x, y, and z 

(i.e., roll, pitch, and yaw) axes, respectively, are most appropriate for 

controlling modes 3, 4, and 5. 


7.3.3 For quick effective suppression of the excessive vibrations in the 
SCOLE configuration excited by the time-minimized BPB slew maneuver, it is 
more appropriate to design the modal dashpots into separate parts than to 
lumping up all the 5 most important modes to be controlled by all the five 
actuators together. High gains not only do not create spillover and 
interaction problems as uaual but rather make the resulting modal dash- 
pots truly powerful and effective for quick suppression of excessive vibra- 
tions. 


7.4 In general, modal dashpots when properly and carefully designed, can add 
desirable amount of active damping to modeled (or "controlled") modes. Unmo- 
deled modes can also receive some concomitant active damping, as a benefit of 
spillover to complement their inherent damping. 
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8. RECOMMENDATIOHS 


We recommend that: 


1. the two-stage approach be accepted as a promising one, and included in Part 
Two of the Design Challenge, for validation using the hardware SCOLE labora- 
tory facility and for comparison with other approaches, and 

2, theoretical and simulation studies on the two— stage approach be continued 
using the mathematical models of both the orbital and the laboratory SCOLE 
configurations for further development of the technology. 


8.1 Careful scientific studies have been successfully conducted on the two- 
stage approach to rapid pointing and vibration control of the flexible orbital 
SCOLE configuration, and the results have been very encouraging. Now that the 
physical SCOLE labortory facility is operational, we recommend that the 
design techniques developed and the technical knowledge gained on the two- 
stage approach be translated to the tethered laboratory SCOLE configuration 
and be tested and validated by the experimental apparatus. Specifically: 

(1) Design a rapid time-minimize bang-pause-bang line-of-sight pointing slew 
maneuver (Stage 1), and a fast effective modal-dashpot type of vibration 
controller (Stage 2), using the mathematical model of the tethered confi- 
guration and the actuators and sensors actually available on the labora- 
tory article. Test the designs on the SCOLE facility in real time. 

(2) Then, conduct a comprehensive sequence of experimental evaluations similar 
to Steps (a) through (e) below. 


8.2 To further develop the technology associated with the promising pratical 
two-stage approach and to gain additional technical knowledge, we recommend 
that studies be conducted on the use of MD — augmented LQG/LTR design of 
vibration control for attaining the specified LOS pointing accuracy. We also 
recommend that the limit on the applied force at the Reflector of the 
orbital SCOLE configuration be lowered by one order of magnitude from 800 lb 
to between 100 and 200 lb, or alternately between 20 and 30 lb, in each direc- 
tion. 

Specifically, we recommend that: 

(1) a series of design, simulation, study and evaluation be carried out on two 
representative cases, 

(2) the total time required from the beginning of the LOS pointing slew man- 
euver to the end of stabilization with the desired 0.02 9 precision be deter- 
mined for each case, and 

(3) a trade-off study be conducted. 


240 


Case 1. Limit set at 150 lb* 


(a) Use a Stage 1 design similar to the one described in Section 3.2.2 for the 
time-minimized pointing slew maneuver. Simulate such a BB slew maneuver on 
the 3-dimensional nonlinear rigid body dynamics of the SCOLE configuration 
first; evaluate the LOS accuracy, and assess the effects of nonzero pro- 
ducts of inertia during the rapid maneuver; compare the results with BPB 
slew maueuver with the 800-lb limit. 

Then simulate this slew maneuver on the flexible-body dynamics of the 
configuration as if it were a time-varying disturbance, and analyze the 
vibrations thus excited. 

(b) Design a similar high-power modal-dashpot type of velocity output feedback 
control (following the same design proceedure as in Section 5). Such a 
vibration control design is to be used, as the first part of Stage 2 for 
suppressing the (moderately) excited vibrations quickly and effectively to 
some desirable low levels. 

(c) Design a "low-power" modal dashpot (MD) type of velocity output feedback 
control first. Augment the SCOLE configuration with the resulting modal 
dashpot design. Then design a LQG/LTR type of optimal state feedback con- 
trol. Such a LQG/LTR/MD control design is to be used as the second part 
of Stage 2 for continuing on suppressing the vibrations quickly to the 
desired LOS pointing accuracy of 0.02°. All force and moment actuators, 
including the two proof-mass actuators, are to be used in both the MD and 
the LQG/LTR/MD designs. 

(d) Simulate the entire Stage 2 design on the SCOLE flexible-body dynamics and 
evaluate the vibration control performance numerically. 

(e) Integrate the Stage-1 and Stage-2 designs (for a continuous operation of 
both pointing slew and vibration control), simulate their application on the 
coupled SCOLE dynamics (i.e., flexible-body dynamics kinematically coupled 
with rigid-body dynamics); evaluate the total LOS pointing and vibration 
control performance and determine the total time required for achieving 
the desired precision. 


Case 2. Limit set at 25 lb 


Use the same Stage 1 design as described in Section 3.2.3, instead of Sec- 
tion 3.2.2, for the time-minimized pointing slew maneuver. Conduct all the 
corresponding sequence of design, simulation, and evaluation as Case 1 
except step (b). 


* In the laboratory SCOLE configuration, the equivalent torque the thrusters on 
the Reflector can generate is about two times the torque producible by the CMG 
on the Shuttle. For the same ratio, the applied force on the Reflector of the 
orbital SCOLE configuration is approximately 160 lb. 
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ABSTRACT 

In this article, the dynamics of slew maneuver of NASA Spacecraft COntrol 
Laboratory Experiment (SCOLE) test facility are developed in terms of an arbi- 
trary maneuver about any given axis. The set of dynamical equations incorporate 
rigid-body slew maneuver and three-dimensional vibrations of the complete 
assembly comprising the rigid shuttle, the flexible beam, and the reflector with an 
offset mass. The analysis also includes kinematic nonlinearities of the entire assem- 
bly during the maneuver and the dynamics of the interaction between the rigid 
shuttle and the flexible appendage. The final set of dynamical equations obtained 
for slewing maneuvers are highly nonlinear and coupled in terms of the flexible 
modes and the rigid-body modes. 

The equations are further simplified and evaluated numerically to include the 
first ten flexible modes and the SCOLE data to yield a model for designing control 
systems to perform slew maneuvers. 
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1. INTRODUCTION 


The primary control objective of the Spacecraft Control Laboratory Experi- 
ment (SCOLE) is to direct the RF Line-Of-Sight (LOS) of the antenna-like 
configuration towards a fixed target under the conditions of minimum time and 
limited control authority [1]. This problem of directing the LOS of antenna- like 
configuration involves both the slewing rmineuver of the entire assembly and the 
vibration suppression of the flexible antenna-like beam. The study of ordinary 
rigid-body slew maneuvers has received considerable attention in the literature 
[2,3] due to the fact that any arbitrary large-angle slew maneuver involves 
kinematic nonlinearities. This is further complicated in the case of SCOLE by vir- 
tue of a flexible appendage deployed f rom the rigid space shuttle. The dynamics of 
arbitrary large-angle slew maneuvers of SCOLE model are derived in this report as 
a set of coupled equations with the rigid-body motions including the nonlinear 
kinematics and the vibratory equations of the flexible appendage. 

The dynamical equations of slewing maneuvers of this large flexible spacecraft 
are developed by writing the total kinetic and potential energy expressions for the 
entire system. The energy expressions are further utilized in formulating 
Lagrange’s equations which are expressed in terms of non-generalized co-ordinates 
using an inertial co-ordinate system and a body-fixed co-ordinate system at the 
point of attachment of the flexible beam to the shuttle. The generic model used f or 
this analysis consists of a distributed parameter beam with two end masses. The 
three dimensional linear vibration analysis of this free-free beam model with end 
masses [4] is incorporated together with rigid-slewing maneuver dynamics which 
are written in terms of four Euler parameters [5] and angular rotation about an 
arbitrary axis of rotation to yield the final set of highly nonlinear and coupled 
equations. In the derivation of the equations, it is assumed that the vibratory 
analysis is for small motions. 



2. LIST OF SYMBOLS 


a (z ) Position vector of mass element on the beam from the point 

of attachment 
B Damping matrix 

C Inertial frame to body-fixed frame transformation 

c_ Position vector f rom the point of attachment to the mass center 

of the beam 

D Mass density of the beam 

d_(z ,t ) Displacement vector of mass element in the body-fixed frame 

E Modulus of Elasticity 

(f ) Force applied at the orbiter mass center 
F 2 (r ) Force applied at the reflector mass center 

Qo (r ) Moment applied about the orbiter mass center 

G y, Modulus of rigidity for the beam 
I Beam cross section moment of inertia 

l x Beam cross section moment of inertia, roll bending 

I y Beam cross section moment of inertia, pitch bending 

I 0 Equivalent mass moment of inertia 

1 x Mass moment of inertia matrix of the shuttle 

/ 2 Mass moment of inertia matrix of the reflector 

J Mass moment of inertia matrix of the beam 

K Stiffness matrix 

L The Length of the beam 

M Angular velocity vector transformation 

M ^ Effective moment applied at the reflector c.g. 

m Total mass of the flexible beam 

m t Mass of the orbiter 

m 2 Mass of the reflector 
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n The maximum number of modes considered 

Q The generalized f orce vector 

q t Generalized coordinates 

R Position vector of the mass center of the orbiter in the inertial 

frame 

r_ Position vector from the orbiter mass center to the point of 

attachment 

r x x co-ordinate of the reflector mass center in the body-fixed 

frame 

r y y co-ordinate of the reflector mass center in the body-fixed 

frame 

T Total Kinetic Energy 

U Total Potential Energy 

u x (z ,t ) The beam deflection in x direction referred to the body- 

fixed frame 

Uy iz ,t ) The beam deflection in y direction referred to the body- 

fixed f rame 

u ^(z ,t ) The torsional deflection about z axis in the body- 

fixed f rame 

V Velocity vector of the mass center of the orbiter in the body- 

fixed frame 

Yo Velocity vector of the point of attachment in the body- 

fixed frame 

p Mass per unit length of the flexible beam 

A. Vector representing the axis rotation during the slew 

maneuver 

<t> xi i th Eigenf unction corresponding to u x 

<t> yi i th Eigenf unction corresponding to u y 

i th Eigenfunction corrsponding to 
9. The attitude of the orbiter in the inertial frame 

f Slew Angle 


The angular velocity of the orbiter in the inertial 
frame 
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ft The angular velocity of the reflector in the inertial 

frame 

4 Damping ratio 
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3. ANALYTICS 


Co-ordinate Systems 

The motion of SCOLE assembly when considered as a rigid body in space has 
six dynamic degrees of freedom: three of these define the location of the mass 
center, and three define the orientation (attitude) of the body. The motion of this 
rigid body is goverened by newtonion laws of motion expressed in terms of 
changes in linear momentum and angular momentum. These relationships are 
valid only when the axes along which the motion is resolved are an inertial frame 
of reference [9,10]. To define the orientation of the orbiter in space, a set of orthog- 
onal axes fixed in the body is utilized. Then the attitude of the orbiter is defined in 
terms of the angles ( 0 1 , 0 2 .^ 3 ) between the body- fixed axes and the inertial co- 
ordinate axes. The body-fixed f rame origin is located at the point of attachment of 
the flexible appendage with the rigid shuttle for this analysis (Fig. 1). 

The transformation from the inertial frame to the body-fixed frame is given 
by the matrix, C as developed in figure 2 where if i , j , k represent the dexteral 
set of orthogonal unit vectors fixed in the body- fixed frame and 0 x is the rotation 
about i , 0 2 is the rotation about j and 0 3 is the rotation about k . These rotations 
are carried out successively as shown in figure 1 and the matrix C is given as 


C = 


cos 0 3 sin 0 3 0 
— sin 0 3 cos 0 3 0 
0 0 1 


cos 0 2 0 — sin 0 2 

0 1 0 
sin0 2 0 cos0 2 


1 0 0 
0 cos 0 x sin 0 ! 
0 — sin 0 j cos 0 1 


Thus C 7 is obtained as 


( 1 ) 


C r = 


cos0 2 cos0 3 — cos0 2 sin0 3 

sin0 3 sin0 2 cos0 3 +sin0 3 cos0 j — sin0 x sin0 2 sin0 3 +cos0 3 cos0 x 
— cos0 jsin0 2 cos0 3 -l-sin0 3 sin0 3 cos0 jSin0 2 sin0 3 +cos0 3 sin0 j 


sin 0 2 

— sin0 icos© 2 
COS0 3 cos0 2 
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( 2 ) 



In order to completely define the attitude (orientation), it is needed to relate 
the rotation angles 0 lt 9 2 , and 0 3 to the angular velocity components (o lt o> 2 . o 3 ) 
of the orbiter. One way of obtaining the required relations is via body-three angles 
method [5] which was utilized in developing C matrix in equation (1) and these 
relations are 


0 1 = ( c^cos© 3 — o> 2 sin03)/ cos0 2 

0 2 = ( OjSin© 3+o 2 cos0 3 ) (3) 

0 3 = ( — o> 1 cos03+o 2 sin03)tan0 2 +o) 3 

Thus, the angular velocity of the orbiter can be obtained in the inertial frame 
by means of the following transformation 


to = M t 9 


(4) 


where the transformation M T is given as 


M t = 


cos0 2 cos0 3 sin0 3 0 
— cos0 2 sin0 3 cos0 3 0 
sin0 2 0 1 


(5) 


Although the body-three angles method is used here for obtaining the 
transformations C and M , there are three other methods which can be used to 
obtain the same transformations. A detailed discussion of all the methods is given 
in reference [5] and a summary of the transformations using the remaining three 
methods is given in the Appendix. 


Kinetic Energy 

If the position vector of the mass center of the orbiter in the inertial frame 
(Fig. 3), R, is given as 
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Rx 

R= Ry (6) 

Rz 

then the velocity of the mass center in the inertial frame is 

Rx 

v 7 (r)= r y (7) 

Rz 

This velocity can be transformed in the body-fixed frame as 

Rx 

V(t ) as C Ry (8) 

R Z 

The velocity of the point of attachment in the body-fixed frame is 


Vo = K+uxr ( 9 ) 

where t_ is the vector f rom orbiter mass center to the point of attachment. 

Defining the position vector (Fig. 4), a., of a mass element on the beam from 
the point of attachment (origin of the body-fixed frame) before deformation as 

0 

a_— 0 ( 10 ) 

z 

and the displacement vector of this mass element as 


u x (z ,t ) 
d_ (Z ,t ) = Uy (z ,t ) 
0 


( 11 ) 
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the position vector after deflection is given as a+d_. The kinetic energy in the beam 
[6] is 


T x = (1/ 2)mVjV 0 + (1/ 2) a T [J]^mVo T \tte+(ll 2) f dj ddm 


+ V j J ddm+o > r J addm +(1/ 2) J u x ‘ Uy ' u ^ j dl 


U r 


U. 


y 

u + 


( 12 ) 


where the vector c_ is f rom the point of attachment to the mass center of the beam 
and if it is assumed that the beam is a thin rod, then it is given as 


c = 


= ( 1 / m ) J adm = 


0 

0 

-LI 2 


(13) 


and using the skew symmetric form for the vector cross product for any two vec- 
tors c_ and to (in the same reference frame) as 

c.x 6> = [ c_ ] 6> 


0 -c 


c = 


— c„ 


x 

0 -C x 

: r 0 


also, the moment of inertia matrix is given as 


J =(1/ 3) pL 3 


1 0 
0 1 
0 0 


0 

0 

0 


(14) 


(15) 


where p is the mass per unit length of the beam. The last term in the equation (12) 
corresponding to torsional motion is given as 


= ( 1 / 2 )/ | u,- i y if 


(.11 2) / | u,‘ u,- if \u 

1/ 2 ipds)s 2 


u T 


zz 


y 

u + 


0 0 
0 1/ 2 ipds^s 2 0 

0 0 0 


u r 


u 


y 

u ^ 


The kinetic enei^y equation (12) can be simplified as 


T i = Cl/ 2) pLVjV 0 + ill 6) pZ, 3 I o, 2 +6>, 2 1-qLV/c 6>+oZ, £ 2 

' ' i = 1 


where 


+V^a+co r j+(l/ 4) p 


Zps. tf+LPeiti 7 

i =1 i =1 


(16) 


(17) 


“x = Z <t> xi (s )?i (f ) 

i = 1 
n 

U y = L <t>yi k )?i (? ) 
i = 1 

V = Z<M*)fi(r) 

i=l 

n 

Uy 2! )?j (f ) 

i — 1 
n 

u^= £ (fripi is )? t (r ) 
1=1 
z 

Pit = f<t> xi (s)ds 
0 
z 

P2i — f <t>yi (s )ds 
L 

Pm = f s<t> xi is)ds 
o 
z 

/>4< = f s<t>yi is )ds 
0 
z 

p 5< = / (^^r] 2 ^ 

= iki 


( 18 ) 


P 6i 


ds 



and 


<kO ) = P 


'LPuti 

t=l 

n 

ZP2i<ii 

i = 1 

0 


(19) 


n 


1=1 


kit) = p 


n 

Zp mQi 


1=1 


0 


( 20 ) 


The expressions for p lt , p 2i , p 3i , p 4i , p 5i , and p 6i are developed as follows. Note 
that 


<t> xi is) = A xi sin/3; s +B xi cos/3; s +C xi sinh/3; s +D xi cosh/3; s 


where /3; = 


0); 2 p 


1 / 4 


El 


Since for SCOLE configuration El x = EI y and = & iy , El and /3; are used 
for both 0 xi is ) and <P yi (s ). However, this may not be true for other 
configurations. 


Pu = 


PiL- 


Pu = f 4> xi is)ds 
0 
L 

Pli = f <t>yi is )ds 
0 

-A xi cos/3; L +B xi sintt; +C X< cosh/3; L j 


+Z) X ; sinh/3; L +A X ; — C x ; 


(21 A) 


Defining a; = 3; L 
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—A xi cosa* +B xi sinai +C xi coshai +D xi sinhaj +A xi —C xi 


(21B) 


Pu = 


1 


ai L 
similarly, 


Pu — 


P 2i 


— — j —Ayi cosft L +Byi sinft L 

ft L, I 


oil L 


+Cyi cosh^j L +Z)yi sinhft Z +Ayi — Cy | 
cosa; +B yi sinai +C y < coshai +0^ sinhai +Ayi —C yi 


and these can be given as 


Pm — f s< t>xi Ms 
0 
L 

Pm — f s<t> yi (s')ds 


Pm 


= A 


xi 


300==t« 


sinft L 

L cos/B< Z 


cosft Z 

Z sinft A i 

1 _ 

ft 

Xi 

oo , 

ft 2 

r _ 

+ A A 2 


C* 


Z cosh ft Z 

3 oo 

sinhai Z 

+Aci 

Z sinhft Z 

I 

cosh/3 f L 

1 

ft 

ft 2 

ft 

ft 2 

ft 2 


Pm - A 


xi 


Z, 2 sina^ L 2 cosa^ 


+B xi 


ap a t 

L 2 coshai L 2 sinha ( 


L 2 cosai L 2 sinai £2 


«< 


a> 


a 


L 2 sinhai L 2 coshai £ 2 


Similarly, 


Pm = 


A^ 


xi 

[ "* 

» 2 

j - 

<*< 

<* 2 

sinft Z 

Z cosft L 

+ 5yi 

cosft L 

Z sinft L 

1 

ft 2 

ft 

ft 2 

+ ft ' 

"ft 2 


a 


+ 


r . 

^yi 


Z cosh ft L 

sinft Z 

4- D 

L sinh ft Z 

cosh/3i Z i i 

ft 

ft 2 

-rVyi 

ft 

ft 2 + ft 2 


Pm = A 


yi 


Z 2 sinai 

Z 2 cosai 

4_ D 

L 2 coso^£ 

Z 2 sinai £2 

“ 2 



a 2 

r 

<*i a 2 


Z 2 coshai 

Z 2 sinha ( 

_l_ n 

Z 2 sinha 4 

Z 2 coshai £ 2 

<*< 

<* 2 



2 ^2 
a i <*< J 


2f9 


(22A) 

C22B) 


C23A) 

(23B) 


(24A) 


(24B) 



and these can be shown to be 



PSi = A xi 


1 

6 2 
—A x i B xi 

+A xl C xi 


2 c°s23i L +( 2 


)sin23* Z 


-^-sin23* Z/ + 
2 


1 Z 2 3* 

4 3* 2 


cos2/3 * L ■ 


4 3* 


3*Z ; 


(cos3* L sinh/3[ Z )+(sin3* Z cosh3* Z ) 


-2Z 


sin3* Z sinh3* Z 


ft 


^’-Axi Z> x < 


3*Z : 


-2Z 


sin3* Z cosh3* Z 

tfL 3 i 


+ 5 X ? 


(cos3< Z sinh3* Z )— (sin3* Z cosh/3 * Z ) 
(cos3i Z cosh 3* Z )+(sin3* Z sinh3< Z ) 

(cos3i Z cosh3* Z )— (sin3* Z sinh3* Z ) }+ 

)sin23* Z 


& 


^cos2 3*Z+(- L 1 


— ^xi Cxi 
+2Z 


3*Z 2 


2 ‘ 2 43* 

(sin3* Z sinh3i Z )— (cos 3i Z cosh3* Z ) 


cos 3i Z sinh3i Z 


_1_ 

ft 


(cos3j Z cosh 3i Z ) 


ft 


■ftci Aci 
+2Z 


+ci 


+CwZ) : 


+(sin3i Z sinh3i Z ) 

3* Z 2 (sin3i Z cosh3* Z )— (cos 3* Z sinh3* Z ) 

cos3* Z cosh 3i Z W (cos 3i Z sinh3i Z ) 

— (sln3i Z cosh3i Z ) 

3*Z 2 , 1 


3* 2 Z ; 


xi-^xi 


3* Z cosh23i Z +( + — )sinh 2 3< Z 

-?ii^cos23< Z -^sinh23< Z +-^-cosh23/ Z — 


43* 


43* 



i B L 2 i B 2 L : 

+Dj J fi ( L COS2 3, L +(- ! -y-+^-)sinh 2 /3 i L + 


P6i =A> ^f-+j yCOs23 i L+C^--^-)sin2|8 f L 


£^ft_ 1_ 

2 4ft 


AyiByi 2 S ^ n 23( L +( 2 )cos2 3;Z, ^ 

+Ayi C yi 3; L 2 1 (cosft L sinhft L )+(sinft L coshft L ) 1—2 L I sinft L sinh ft L 


(cos/3; L sinhft L )— (sinft Z, coshft L ) 


+A yi Z> yi 3; L 2 (cos3, Z. cosh3; Z, )+(sin3; L sinh3; L ) |— 2 L sin3; Z, coshft L j 

— -j— (cosft Z, cosh 3; Z, )— (sin3; L sinh3; Z. ) + -j- 
3; 3; 


+V ^T~-J fcos 2 fti+(Z^i— _i_)sin 2 fti. 


— Cy; fiiL 2 (sin 3< Z. sinh3; L )— (cos3; L cosh3; Z. ) +2L cos3; L sinh3; L 


— (cos 3; Z, cosh 3; L )+(sin3; L sinh3; Z, ) +-J- 
3; 3; 

1 

— -ft* Z)^ fiiL 2 (sin3 ; L cosh3; L )— (cos 3; L sinh3 ; L ) +2Z. cos3; Z, cosh 3; Z, 

— -j— (cos 3; Z, sinh3; Z. )— (sin3; Z, cosh3; L ) 

3; | I 

1 I 3 L 2 i fi 2 L 3 I 

Cyi — fi t L cosh2 fiiL+{— l — — +-^j— )sinh 2 3;Z. l — — 


Cyi Dyi 


fiiL 2 


'JLt 1 1 1 

— — cos2 3; L — -^-sinh2 fi t L +-^— cos2 3; L — 


1 ft r 2 i o2r3 ) 

■=■ 3; L cos23; Z, +( +~r4 )sinh 2 3; Z, + ~r — . 


2 4 ft 


The equations (25) and (26) can alternatively be derived by replacing 3; = — — • 

i * 


The kinetic energy of the reflector is 


2G0 



T 2 = (1/ 2) m 2 VjV X) -m 2 yJ^L)^m 2 }^d_{L )-(l/ 2) m^dJiL )2(L )o 

+m 2 co 7 ’aa)i.a)+Cl/ 2) TtijdJ {L )d_{L )+(l/ (27) 

where m 2 is the mass of the reflector and 1 2 is the mass moment of inertia matrix 
of the reflector. The deflection vector d_{L ) at the mass center of the reflector is 
given as 


d{L) 


u x (L) — r y u y >( L ) 

Uy(L) + r x u^ 
u x XL )r x + n y tL)r y 


(28) 


and the position vector from the point of attachment to the reflector mass center is 
given by 


aCL) = 



(29) 


Thus, 


iu) 


u x [L) — r y ii L ) 
Uy(L) + r x ii ^(L ) 
u x '( L )r x + ii y XL )r y 


(30) 


The angular velocity of the reflector in the inertial co-ordinate system ft can be 
shown to be 


ft = 01 + 


Ux 

V 

U0 


(31) 


The equation (27) can be simplified as 

r 2 = (l/ 2) m 2 VjV c —m 2 V^{L)^+m 2 yJd_iL )+(!/ 2) m 2 L z [o>i 2 +<o 2 2 ] 


+m a_(.L ')d_{L )+( 1 / 2) m 2 


Z Z^(^)<^xy(^>7i^ + 

<=U=l 
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+(1/ 2)P 7 7 2 P+(1/ 2)(Jl&L 


(32) 


n n 

II 


t=u=i 


(t>yi (7, )0y j (7, )q t qj 


where 


P T = 

I $xi ("^ )?i (? ) £ 0yi (7. )^i £ (Z, )<7^ (f ) 

< =1 t=l i=l 


The kinetic energy of the shuttle, T 0 , is given as 

T 0 =(1/ 2)m 1 V r V:+(l/ 2) o/ |/jjw 


where m j is the mass of the shuttle and 7 : is the mass moment of inertia 
of the shuttle. 

The total kinetic energy is given as 

T = To+Ti+rj 


This can be simplified as 

T = (1/ 2) m 0 V T V+oS [zzjv+d/ 2) to r [z o jtt+pZ.f^+V^a 
+6> r r ct- -l-fa> r B+m j (L )+m lo/rd (L )+ 

m 7 <j} T a ( L )d (7. )+(!/ 2) m 2 

+(1/ 2)P T l2PHll 4) p 


£ <*> 2 U )+</>£ (L) 


(=i 

n 


I/>5itf< +LP6i<ii 
(=1 i=l 


<?( 


where 


m 


H = 


7 0 = 7 !+(l/ 3) pL 3 


•o = m i+pL +m 2 
pL +m 2 Jr.+m jS (L )+pZ£ 


1 0 0 
0 1 0 
0 0 0 


+7 2 +J o7,r? —oLrc —m -,rr — m ? ra (L 


2S2 


(33) 


(34) 

matrix 


(35) 


(36) 



The term J 2 in. this equation can be shown to be: 


/-> = m- 


(r 2 +L 2 ) — r x , 


r x L 


~r x r y C rf+L 2 ) r y L 


r x L 


r y L (r x 2 +r y 2 ) 
The total kinetic energy expression can be further simplified as 


where 


T = (1/ 2) m 0 V T V+ 


hJk+Ci/ 2 )w 7 ’ 

+< o r 1 ^ 2 1 / 2)_£ 7 | a 3 


/ 0 Jco+v: 7 * [AjJi C37) 


|a 1 Jg. = a+mjdCL ) 

A 2 Jd = r q+ E+m ^fd (L )+m 7 a (L )d.(L ) 


A, = 


X 0 
pL +m 2 +p 5 i +^6< 
0 \ 


[ (0'a)) r / 2 [0'a))] 


In this equation 


<t> u '{L ) 0 0 

0 <t> ly '(L ) 0 

0 0 <l>i$(L ) 


0'U) 


d>u '(L) 0 0 

0 4> iy 'U ) 0 

0 0 <*> iV ,(L) 


Here i-2,3 n. The number n indicates the total number of flexible modes con- 

sidered. 


Equations of motion 
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Lagrange’s equations of motion for the case of independent generalized co- 
ordinates q k are 


±JSL _ J3I. = <v _ ML 

& ftq k 0?* 0?* 


Gt =1.2, n) 


where, T = is the kinetic energy 

U = £/ ( 2 .) is the potential energy, and 

Q k are the generalized forces arising from nonconservative sources. 


(38) 


The generalized co-ordinates are: 

R x ,Ry >Rz ~ position of orbiter mass center relative to inertial frame origin. 
0 1,0 2 $ 3 “ roll, pitch and yaw angles of orbiter. 

9 1,92 9 n ~ modal deformation co-ordinates for the beam. 


The previous kinetic energy expression developed in equation (37) is given in 
terms of nonholonomic velocities V and and generalized velocities g_- Using the 
notation T(V ,(o,q ) for this kinetic energy expression and T for kinetic energy 
expression in terms of generalized velocities, the equations of motion are developed. 
Thus, equation (37) is rewritten as 


T = (1/ 2) m 0 V T V + 0 / 

+ fi> r | A 2 


h]v + 
i + ill 


(1/ 2)o/ |/ 0 ]<o + Vj A 

2)i. r 


i 


(37) 


(a) Translational Equations 

From the chain rule applied to equation (37) using equation (8), one gets 




_a l. 

&Rx 


a^i 

J£- 

= C T 

_aL 

ftRy 


0^2 

J£_ 


sL 

bRz 


bv 3 


Z64 


( 39 ) 



Also, the generalized forces are CF_ it ) where 


£CO = £UHF 2 (r) (40) 

Fg (r ) represents the foroe applied at the orbiter mass center and F z it ) represents 
the force applied at the reflector mass center. From Lagrange’s equations 


d 

af 

+ cc T 

bT 

dt 

bZ 


bZ 


= Fit ) 


(41) 


and from equation (37) 


ST 

bZ 


— m 0 V—Hai+A 


(42) 


Substituting equation (42) in (41), 


m 0 V-Hto+A i ii = -CC T im 0 V-H(±+A l g_) + Fit ) 


(43) 


This can be rewritten as 


m 0 V_ — H in -¥ A {d = tL\ + F it ^ (44) 

where the nonlinear term Ni is given as 

N x = -CC r im 0 V-H(o+A i i') (45) 

= — oj( m 0 Y_—H A \Q ) 

Here, oj = CC T . 

(b) Rotational Equations : 

From equation (4) 


to = M r 6 

Again using the chain rule 


3T' 

= M 

< 

af 

ai 


6“ 


26 S 


( 46 ) 


Also 


bT 


0v: r 


0a r 

00 i 


00i 


00i 

bT 


0V r 

af 

+ 

0« r 

00 2 


002 

0Z 


002 

ar 


0V 7 


&L 

00 3 


003 


003 


ar 

a« 


It can be shown that 


A vL = V T C &L 
a o, - 00, 


.i =1,2,3,. 


and 


=u t M~' (=123 

ae, - ae. 


and 


al 


V r c - 
- 00 1 

i l 

001 


af 

+ 


002 

az 


00 2 

V r C-^— - 


C) r M -i^M 

- 003 


003 


From equation (37), 


bT 

0<i> 


= 77V + I 0 o) + A iQ 


and as before 


af 

0o> 


(47) 


(48A) 


(48B) 


(49) 


(50) 


bT 

bv 


— m oY- ~ 77 o) + A ig, 
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(42) 



Using the Lagrange’s equations 


d_ 

dt 


ST 

61 


-f =M£ 


(51) 


where Q, is the net moment about the mass center of the orbiter with respect to 
the body-fixed frame. It is given as 


a = Go + (r+o.)x£ 2 (52) 

Qt, is the external moment applied about the mass center. Eqation (51) can be 
simpified by substituting equations (42), (49), and (50) together with the relation- 
ship developed in (46) as 


HV +/ 0 w + A& = Q_ + N 2 
where the nonlinear term N 2 is given as 



* TcJ k 






1 

11 

^1 

zTc! k 

ST 

SZ 

+M" 1 



—M 


F c -aM 







Cc) Vibration Equations of the Beam 


ST 


(53) 


(54) 


Since T in equation (37) is given in terms of l which is a vector of general- 
ized velocities, 

6 ^ _ ST 
Si Si 


and 


SL 

Si 


= a{v + aI& + a& 

2C7 


( 55 ) 


The potential energy in the beam is given by 


U =(1/ 2 )gj Kg_ 


where the stiffness matrix K is given as 


\° 

K = k u 


k u = El /3j 4 / (s )rfs + f )ds + G f <t)}i (s )ds 


, ,i 

Z) COi 2 2 

G j, represents the modulus of rigidity of the beam and /3,w = — — — where D 

I 

is the mass per unit volume (mass density) of the beam. Thus, 


f- -ja. . 

32. j 

Using the Lagrangian Equations (38) and assuming that F 2 = 0. , 


A jV + A 2 <A + A $£_ = —Kg. 


matlons 


If it is considered to perform a slew maneuver about an arbitrary axis X. and 
the slew angle to be then the slew maneuver can be expressed in terms of four 
Euler parameters. These four Euler parameters are defined as 


_L = *2 = A, sin- 
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( 61 ) 



and their derivatives with respect to time are given as 


if = i(€4 ffi + X*!*) 


d €4 
~dt~ 



(62) 

(63) 


The four Euler parameters can be related to the angular velocity components 
of the rigid assembly as 




*1 €4 —€3 €2 


0 

*2 


e 2 € 3 € 4 — € 1 



€3 


e 3 — € 2 € 1 e 4 


<0 2 

*4 


€4 — fi —€2 ~ € 3 


0)3 


(64) 


If a slew maneuver is considered to be purely rotational, then the transla- 
tional velocity and acceleration can be shown to be negligible during the slew 
maneuver and only the rotational and vibration equations are reqired for the 
analysis and they are simplified by setting V = 0. in both (53) and (59) and are 
written as follows 


I 0 o) + A 2 'd_ — G. (f ) 3" tL 2^0) (65 ) 

Al<k+ A& + A^ = a(r) . (66) 


where, 

£(r ) is the net moment applied about the mass center of the orbiter and is 
given by the following equations (figs. 1& 2) 

G(r) = G 0 (r) + (r + a)xF 2 . (67) 

Also, Q(r ) represents the generalized force vector which is given by the following 

2G9 



equation 


Z c <w* ) + Qjyfi )) + Q*, + Q y , + Q* 

Z C Qjx 2 (t ) + Q^yjC 7 )) + Qx 2 + Qy 2 + Q^J 


aco= 


Z c Q yXj (r ) + QjySt )) + Q xi +Qyi+ Q 0 , 


( 68 ) 


where, the generalized force components are given as 


L 

Qjx, = f F jx (z >t )S(z -Zj )<t> xi (z )dz (69) 

0 
L 

Qjy, = f F j y (z ,t )S(z —Zj )<t>yi (z )dz (70) 


and 

Qj+ l (T) = 0 . (71) 

Here, Fj X (z t ) is the x component of the concentrated force applied at location j 
on the flexible antenna and F Jy is the y component of that force. 

Also, 


Q*i(f) = F 2x (f )<t> xi (L) 

Qyi(t) = F 2 y U)<t> yi U) (72) 

Q ^(f ) = M ^ (.t (.L ) • 

Here, F 2 is the force applied at the reflector C. G. 
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Thus, 

M^(f ) = Fix + F 2y r JC + M 2 ^ • (73) 

The location of reflector C. G. is given by coordinates (r x ,r y ) and M 2 ^ 
represents the external moment applied at the reflector C. G. 

Thus equations (62) - (73) completely represent the dynamics of the slew 
maneuver. These equations are nonlinear and coupled including both the rigid- 
body dynamics and the dynamics of the flexible appendage with kinematic non- 
linearities. It is important to note that the nonlinear term N 2 (co) is dependent on 
the rotational velocity and as a result determined by the slew maneuver rate. Thus 
the basic slew maneuver stretegy has to be developed before this term can be 
linearized. 


Vibration Enuattons of the Beam with Damping 


If damping is included in the derivation of vibration equations of the beam, 
then the damping effect can be expressed in terms of frictional forces. These are 
nonconservative, retarding forces and are assumed to be proportional to the gen- 
eralized velocities. In deriving the vibration equations by means of Lagrange s 
equations, the following function is introduced 


z (= 1;=1 


It also has a positive definite quadratic form similar to the kinetic and poten- 
tial energy expressions. 


With this definition, Lagrange’s equations assume the form 

= q. - -ML (* =1,2 n) (75) 

dt ftq k Qq k d ik St* 
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Again, as before 


$r 

bi 


= A\V_ + AjO) + A 


(55) 


and 


&L 

ba. 


= k<l 


and it can be seen from (74) that 


b£d 

bi 


= B l 


where the damping matrix B is symmetrical and is given as 


B = 


*n b 12 . . b ln 

*21 *22 • • *2n 

• • i • • 

• • * « • 

*n 1 *n 2 • • *nn 


The vibration equations are given as 


(58) 


(76) 


(77) 


A{V +A$u + A& + Bl = -K& + Q(t) (78) 

The slewing equations (65) and (66) would be modified as 

I 2<k ~l- A 2<L = G (f ) + N ->( ta,<7 ) (79) 

A 2 u> + + Bg_ = — Kg_ + Q.(r ) (80) 

Nonlinear Term in the Rotational Equations 


The nonlinear term N 2 in the rotational equations (65) and (79) during the 
slewing maneuver is simplified as 
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— A/f-1 


n 2 = m 


* M ' l X 


- M 


/ O co + A 22. 


(81) 


where 


_ 0 0 0 
9 y i 1 


(82) 


/ir'M = 1 

00 2 COS0 2 


(— o^sinG 2 cos 2 0 3 +ci> 2 sin0 2 sin0 3 cos0 3 ) (ojjSinO 2 sin0 3 cos0 


— ci) 2 sin0 2 sin 2 0 3 ) (cojeosG 2 cos0 3 — o) 2 cos0 2 sin0 3 ) 


(83) 


u? M 1 (co 2 cos0 2 ) (— oi 1 cos0 2 ) 0 

Q(7 3 COSC7 2 


(84) 


Since the transformation matrix, M , is a function of 0 2 and 0 3 , the time 
derivative of M can be expressed by the chain rule as 




(85) 


From equation (5) 




(— sin0 2 cos0 3 )0 2 (sin0 2 sin0 3 )0 2 (cos 0 2 )02 
0 0 0 
0 0 0 


( 86 ) 




(— cos0 2 sin0 3 )0 3 (— cosO 2 cos0 3 )0 3 0 
(cos0 3 )0 3 (— sin0 3 )0 3 0 

0 0 0 


Substituting these equations (86) and (87) in (85) 


( 87 ) 
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(— sinG 2 cos0 3 )0 2 +(— cos 0 2 sin0 3 )0 3 (sin0 2 sin0 3 )0 2 +(— cos0 2 cos0 3 )0 3 (cos0 2 )0 2 
= (cos0 3 )0 3 (— sin0 3 )0 3 0 

0 0 0 

( 88 ) 

From equation (4), this can also be expressed as 


(— sin0 2 cos0 3 )(w 1 cos0 2 sin0 3 (sin0 2 sin0 3 )(co 1 cos0 2 sin0 3 

+o 2 cos0 2 cos0 3 )+(— cosG 2 sin0 3 ) +o> 2 cos0 2 cos0 3 )+(— cos0 2 cos0 3 ) 


1 

COS0 2 


(— WjSinG 2 cos 0 3 +co 2 sin0 2 sin0 3 
+o 3 cos0 2 ) 

( cos 0 3 )(— ojjsin© 2 cos0 3 
+o> 2 sin0 2 sin0 3 +o> 3 cos0 2 ) 


(— OjSinG 2 cos0 3 +<o 2 sin0 2 sin0 3 
+co 3 cos0 2 ) 

(— sinG 3 )(— OjSinG 2 cos0 3 
+<o 2 sin0 2 sin0 3 +o 3 cos0 2 ) 


cosG 2 (o> 1 cos0 2 sin0 3 +<a 2 cos0 2 cos0 3 ) 


Also, M 1 is given as 


cos0 3 cos0 2 sin0 3 — sin0 2 cos0 3 
M~ l = — ^ — sin0 3 cos0 2 cos0 3 sin0 2 sin0 3 

0 0 COS0 2 


Thus, the nonlinear term N 2 can be rewritten as 
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N 2 = a ’ 3 (q,(L) I 0 cd+A 2Z 


Where the term A '3 is 


A ',(*£) = M" 1 


0 



- M 




N 2 = A '3(<o,0_)/ o cd + A '3(0), 6jA 22. 

= A 4 CcD,e3 + A 5 CcD,e;i. (9i) 

where A 4 depends on the rigid-body slewing and is nonlinear in terms of cd and Q_. 
The second term relates the coupling between the rigid-body slewing and the flexi- 
ble modes. This equation can be further simplified in terms of Euler parameters by 
relationships developed in appendix II as 


N_2 = A 6 ( CD. 6 ) + A 7 ( GiX ) i (92) 

where jl is the Euler vector comprising all four Euler parameters, 

From equations (65) and (66) and by defining A = A 2 I 0 1 A 2 + A3 , the 
following equations are obtained 


cd = h 1 


A 2 A l Bg_ + A 2 A l Kg_ + 


A,A" 1 A 2 7 '/ + I 


Git) 


A 2 - 1 A 2 r / 0 “ 1 +/ 3 


HiifAs) 


(93) 


£ = A -1 Bg_ - A -1 Kg_ - A~ 1 A 2 1: I 0 _1 a(t ) + A~ 1 A 2 T I 0 1 


+ A -1 Q. it ) . 

27T 


( 94 ) 


It is assumed that control foroes applied for vibration suppression has negligi- 
ble effect on rotational maneuver of the spacecraft in developing equations (93) 
and (94). Also, / 3 represents 3x3 identity matrix in these equations. 


The analytics developed in the previous section are utilized together with the 
basic SCOLE data [l] and the three dimensional linear vibration analysis [4] to 
generate the following numerical data. 

= 6366.46sh/gs. ; m 2 = \2A2slugs . ; p = 0.0955 slugs I ft . ; L = 130/f. 

G^ = 7.2E+&lb I ft 2 ; (El ) x = (EI) y =(£/) = AE+Vb-ft 1 -, 

0.036 
jr as -0.036 

-0.379 

0 

c_ = 0 

-65.0 

905443.0 0.0 145393.0 

7j = 0.0 6789100.0 0.0 

145393.0 0.0 7086601.0 

18000.0 -7570.0 0.0 

1 2 = -7570.0 9336.0 0.0 

0 0.0 27407.0 

The three dimensional vibration analysis is given in terms of the first ten 
modal frequencies and mode shapes in table 1. Here, 

, . . aiS aiS otjS ac,s 

<t>xi (s ) = A xi sin— — \-B xi cos— — I- C xi slnh— — \-D xl cosh—— 

j—* Ju JU 

Z76 



(*£ ) ^ ypi sinof^j 4 " B cosof^j 

JL 4 / i 


a 


i ~ 


(js t z pL 4 


j_ 

4 


£7 


a V-i = 


DL 2 <i >, 2 
G 


2 


Using these data the following matrices are obtained. 


/o 


1216640 -1.530307 175667.1 

-31.66433 7082976 -52474.84 

175690 -52503.9 7131493 
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TABLE 1 


FIRST TEN FLEXIBLE MODES OF SCOLE MODEL 


THREE DIMENSIONAL MODE SHAPE CHARACTERISTICS 

MODE No. 

1 

2 ; 


0.27804240E+00 

0.31357296E+00 I 

a 

0.12012084E+01 

0. 127565 18E+01 

a 

0.16282665E+00 

0.38855291E-02 

A 

-0. 1 96 70286E+00 

-0.14998387E-01 

c x 

-0.16983450E+00 

-0.43321018E-02 

A 

0.19616259E+00 

0.14985820E-01 

Ay 

-0.10274618E-01 

0.14219781E+00 

By 

0.575791 33E-02 

-0.22695797E+00 

c, 

0.11810057E-01 

-0.19283105E+00 

Dy 

-0.57220462E-02 

0.22644561E+00 


0. 193609 55E-01 

0.21835058P-01 

A + 

-0.50748354E-01 

0.31115282E-01 

B* 

0.13978018E-04 

-0.75992337E-05 

MODE No. 

3 

4 

FREQ. (HZ.) 

0.81300189E+00 

0.11856099E+01 

a 

0.205403 8 7E+01 

0.24804687E+01 

a 

0.40868 188E-01 

0.80641794E-01 

A 

-0.61958845E-01 

-0.67233377E-01 

C x 

-0.41309992E-01 

-0.80913938E-01 

A 

0.61880796E-01 

0.671063 16E-01 

Ay 

-0.22438404E-01 

0.13728679E+00 

By 

0.36509234E-01 

-0.11 746932E+00 

Cy 

0.24390447E-01 

-0.14085209E+00 

Dy 

-0.36464758E-01 

0. 1 1 72505 7E+00 


0.5661 1842E-01 

0.82557693E-01 

! A * 

0.92698901E-01 

-0. 16 158934E-03 

B* 

-0.87320799E-05 

0.10437718E-07 

MODE No. 

5 

6 


0.20536300E+01 

0.497 16090E+01 

a 

0.32645546E+01 

0.497 16090E+01 

A x 

0.99278129E-01 

0.45739784E-01 

A 

-0.92344553E-01 

-0.46365581E-01 

A 

-0.99442145E-01 

-0.45763106E-01 

A 

0.92225801E-01 

0.463296 76E-01 

Ay 

-0.57396019E-01 

0.78612940E-01 

By 

0.53976008E-01 

-0.79952853E-01 

Cy 

0.58114853E-01 

-0.78914485E-01 

D y 

-0.53906980E-0 1 

0.79891039E-01 


0.14300062E+00 

0.33165303E+00 

A + 

-0.165 886 14E-02 

-0.93 394833E-05 



0.61861804E-07 

0.15017211E-09 
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| THREE DIMENSIONAL MODE SHAPE CHARACTERISTICS 

MODE No. 

7 

8 


0.55157833E+01 

0.12281 249E+02 


0.53501560E+01 

0.79833305E+01 

A x 

0.81311804E-01 

0.44835061E-01 

B x 

-0.82056569E-01 

-0.44834914E-01 

c x 

-0.81344923E-01 

-0.44840508 E-01 

D x 

0.81997259E-01 

0.4481 3 000E-O1 

a 7 

-0.47145439E-01 

0.77404 756E-01 

B y 

0.47703590E-01 

-0.77465629E-01 

C y 

0.47289807E-01 

-0.77475327E-01 

D y 

-0.47669155E-01 

0.77427782E-01 


0.38408 110E+00 

0.85518 143E+00 

a\ 

-0.23855560E-02 

0. 1583037 IE-05 

B* 

0.33122 041E-07 

-0.98715017E-11 

MODE No. 

9 

10 


0.12890442E+02 

0.23679520E+02 

■ 

0.81789349E+01 

0.11085347E+02 

a* 

0.78 743 585E-01 

0.44348498E-01 

B x 

-0.78755259E-01 

-0.44367373E-01 

c x 

-0.78752483E-01 

-0.443505 11 E-01 

D x 

0.78717693E-01 

0.44351 763E-01 

A y 

-0.45569244 E-01 

0.76707490E-01 

By 

0.45609474E-01 

-0.76762782E-01 

c. 

0.45607884E-01 

-0.76733612E-01 

D y 

-0.45587726E-01 

0.76735779E-01 


0.89760 145E+00 

0.16488784E+01 

A 

0.94995 483E-03 

-0.5 1 105957E-06 

_J*± 

-0.56437766E-08 

0.16528495E-11 
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0.45879E+2 
0.36305E-1 
-0.89042E-1 
-0.14067E0 
— 0.1457E0 
0.1914E-1 
0.845 97E-1 
— 0.6893E —2 
-0.4269E-1 
0.4204E-2 


0.36305E— 1 
0.6211E+2 
0.11263E0 
-0.1471E0 
— 0.5518E— 1 
0.19839E— 1 
0.3935E-2 
— 0.7165E— 2 
0.5969E-2 
0.41227E— 2 


— 0.89042E— 1 
0.11263E0 
0.32 73 7£ +2 
— 0.6392E— 1 
-0.14526E0 
0.7925E— 2 
— 0.8369E— 1 
— 0.2829E— 2 
0.8976 7E—1 
0.1866E-2 


-0.14067E0 
— 0.147 IE 0 
— 0.6392E — 1 
0.254 7E +3 
0.1908EO 
— 0.4278E— 1 
-0.761 15E-1 
0.1543E-1 
0.285 9E—1 
-0.906 7E -2 


0.1914E-1 
0.19839E— 1 
0.7925E— 2 
-0.4278E— 1 
-0.257QE-1 
O.23209E+5 
O.10383E— 1 
-0.2089E— 2 
-0.3955E— 2 
0.1227E— 2 


0.84597E— 1 
0.3935E— 2 
— 0.8369E— 1 
-0.761 15E-1 
— 0.12912E0 
0.10383E— 1 
0.55561E+5 
— 0.37286E— 2 
— 0.3859E— 1 
0.23 97E -2 


-0.68 93E -2 
— 0.7165E— 2 
— 0.2829E— 2 
0.1543E— 1 
0.92 22E -2 
— 0.2089E— 2 
— 0.37286E— 2 
0.1342962E+8 
0.1421E— 2 
— 0.4427E— 3 


— 0.4269E— 1 
0.5969E— 2 
0.89767E— 1 
0.2859E— 1 
0.461 IE -1 
— 0.3955E— 2 
-0.385 9E-1 
0.1 42 IE -2 
0.20956 72E+ 8 
— 0.9108E— 3 


— 0.2133821E0 
0.380892 IE +3 
— 0.1808478E+3 
0.1423380E+3 
— 0.2416743E+2 
-0.68022 73E0 
0.2784792E+2 
0.78428 18E + 1 
— 0.2694455E +2 
— 0.9225328E — 1 


— 0.368705 7E+3 -0.7253901E- 
— 0.3030935E +2 -0.8427658E- 
— 0.1318596E+3 -0.125799E0 

-0. 11 3585 IE + 1 -0.236735 IE - 
0.574383E +2 -0.9150328E- 

0.3104929E 2 -0.384306 2£- 

0.6651585E+2 0.596075E-1 

— 0.1930097E+2 -0.4363533E- 
— 0.5544252E +2 -0.4200623E- 
0.1594045E+2 -0.1626004E- 


0.1457E0 
-0.5518E — 1 
— 0.14526E0 
0.1908E0 
0.8103E +3 
-0.2570E— 1 
-0.12912E0 
0.9222E— 2 
0.461 IE -1 
-0.594 7£ -2 


0.4204E— 2 
0.4 12 7E -2 
0.1866E— 2 
-0.906 7E -2 
-0.594 7E -2 
0.1227E— 2 
0.2397E— 2 
-0.442 7E -3 
— 0.9108E— 3 
0.8662547E + 10 

1 

1 

1 

1 

1 


1 

1 
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The stiffness matrix K is calculated using equation (57) and the mode shape 
coefficients given in Table 1. This matrix is a diagonal matrix and is represented in 
terms of the diagonal elements as 


K = 


* u = 0.282021 7£0 

* 2,2 = 0.3574692E0 
* 3 f3 = 0.241 2 807£ 1 
* 4,4 = 0.5285 11 6£ 1 
* 3>5 = 0.1588654£2 
* 6 , 6 = 0.8573860E2 
* 7 , 7 = 0.1146118£3 

* 8 ,8 = 0.5686 10 1£ 3 
* 9 , 9 = 0.6254598£3 
^ 10,10 = 0.2114612£4 


The damping matrix B used for this analysis is a diagonal matrix and for 
damping ratio 4 = 0.003, it is calculated to be 


B = 


b 1A = 0.9685964£-3 
b 2 ,2 = 0.1088608£— 2 
b 3,3 = 0.2834016£— 2 
5 4 ,4 = 0.4256808£— 2 
b 5,5 = 0.7387177£— 2 
b 66 = 0.1719014£— 1 
* 77 = 0.1984237£— 1 
^, 88 = 0.442 1234£—1 
b 9,9 = 0.4633434£— 1 
b io,io = 0.8527647£— 1 
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APPENDIX I 

The following is a summary of transformations between inertial frame and 
body-fixed frame. Here, s t and c t (i = 1,2,3) denote sinO* and cos0( (i =1,2,3) 
respectively. 

(a) Space-three Angles 

c 2 C 3 c 2 s 3 ~ s 2 

C = S1S2C3— s 3C 1 S X S 2$ 3"bC 3C 1 S\C 2 

c I s 2 C 3+ s 2 s 1 c I s 2 s 3~ c 3 s 1 c l c 2 

1 0 S 2 

M T = 0 C ! SjC 2 

0 -s x c 1 c 2 

(b) Space-two Angles 

c 2 s 2 s 3 ~ s 2 C 3 

C = 5 1 s 2 — SjC 2 S3+C3Ci $30^3+5305 

c l s 2 — c l c 2 s 3 — c 3 s 1 C 1 C 2 C 3 — 5 3 5 1 

1 0 c 2 

M T = 0 C ! SiS 2 

0 -s x c x s 2 

(c) Body-two Angles 

C 2 s l s 2 ~ c l 5 2 

C = s 2 S 3 — ^ iC 2 S 3+C 3C 1 CjC 2 S 3+C 3S 1 
s 2 C 3 — 5 l c 2 C 3 — 5 3 C 1 c l c 2 C 3~ s 3 s 1 

c 2 0 1 

M t = s 2 s 3 c 3 0 
• S 2 C 3 ~ s 3 0 
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APPENDIX II 

The transformation that relates the orientation angles 0_ to Euler parameters ± 
is a nonlinear transformation. This transformation is developed for body-three 
angles representation in this appendix and similar transformations can be derived 
for other three representations, namely space-three angles, space-two angles, and 
body-two angles. 

(a) For sin0 2 ^ 1 : 

If-— < 9 2 then 

2 2 2 


02 = 

If ( cos0 x cos0 2 ) ^ 0, then 

9 1 = sin -1 

If ( cos0 1 cos0 2 ) < 0, then 

9 1 = tt — sin' 

If ( cos0 2 cos0 3 ) ^0, then 


sin 1 2( 63 c! + e 2 € 4 ) . 


—2 (€ 2 € 3 — €1 €4 ) 

COS 

sin 1 

2 ( €361 + €2^4 ) 




-2 (c 2 €3 - €4 ) 

COS 

sin 1 

2( €3 € X + € 2 € 4 ) 

^1 



(A.l) 


(A. 2 ) 


CA. 3 ) 


0 3 = sin 1 


cos 


~2 ( € 2 ~ € 3 € 4 ) 

sin -1 j 2C €363 + € 3 € 2 ) 


If C cos0 2 cos0 3 ) < 0, then 


(A. 4 ) 


2BS 


-2 (€^2 ~ €3*4 ) 

COS Sin -1 2( € 3 €! + € 3 £ 2 ) 

(b) For sin0 2 = ± 1. 0 2 is a constant. For sin0 2 =l, 0 2 =-^. However, if 

sin0 2 = — 1, then 0 2 = — For this case, if ( sin0 1 sin0 2 sin0 3 + cos0 3 cos0 ^ ) ^0, 
then 

9 1 = sin -1 2 ( e 2 e 3 + €j€ 4 ) j . (A.6) 

If ( sin0 L sin0 2 sin0 3 + cos0 3 cos0 1 ) < 0, then 

0 i = IT — Sin -1 | 2 ( € 2 € 3 + €!£ 4 ) J . (A.7) 

For this entire case, 0 3 = 0 . 


(A.5) 
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FIGURE 1 
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Figure 3- Position Vectors in Inertial Frame 
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Figure 4- Vectors in Body-fixed Frame 
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Optimal Planar Slewing of the Flexible Orbiting SCOLE 


P. M. Bainum and Feiyue Li 
Department of Mechanical Engineering 
Howard University, Washington, D.C. 20059 


Abstract 


The nonlinear equations for the planar motion of the 
flexible SCOLE are derived by using Lagrange's equations. The 
displacements of the flexible parts are assumed small as 
compared with the SCOLE dimensions. The linearized version of 
the motion equations is obtained. The Maximum Principle is 
applied to the planar maneuver of the SCOLE to obtain the 
associated optimal control. The resulting nonlinear two-point 
boundary-value problem is solved by using the quasilinearization 
method, in which the solution of the linearized version is used 
as a starting solution. Some numerical results are presented to 
show the application of this method. 
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Introduction 


Pontryagin's Maximum Principle has been used for the optimal 
attitude maneuvers of a spacecraft in the following cases: 

1) Rigid (3-D, multi-control source) 

(Junkins, J.L. and Turner, J.D., 1978) 
(Bainum, P.M. and Li, F. , 1986-1988) 
(Lin, Y.Y. and Kraige, L.G., 1987) 


2) Flexible 

a) Linearized Eq. ( single-axis rotation, rigid-hub-2 

symmetric beams; Breakwell, J.A., 1979) 

b) Nonlinear Eq. (single-axis rotation, rigid-hub-4 

symmetric beams, 1 control torquer; 
Junkins and Turner, 1980) 

c) Nonlinear Eq. (single-axis rotation, rigid-hub-4 

symmetric beams, 5 control torquers; 
Turner and Chun, 1984) 

d) The Present Problem: 

Nonlinear Eq. (single-axis rotation. Shuttle-beam- 
reflector, control torquers and 
control forces) 


293 



The methods used in solving the two-point boundary-value 


problem are : 

1) For Rigid Spacecraft: 

a) Differential Correction (of unknown initial costates) 
and Relaxation Process (to increase the participation 
of the nonlinearity in the solution) 

(Junkins and Turner) 

b) Hybrid Approach (direct gradient method and the method 
of particular solutions) (Lin and Kraige) 

c) Quasilinearization Method (Bainum and Li) 


2) For Flexible Spacecraft: 

a) Linear Eq. 

Transition Matrix Method (closed form solution of the 
unknown initial costates) 

b) Nonlinear Eq. 

Differential Correction and Relaxation Process 
(Junkins, Turner, and Chun) 

Quasilinearization Method (Bainum and Li) 
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MAXIMUM PRINCIPIF 


STATE EQUATIONS 

x - f(x) ♦ B(x)u, x(0)-x Q , x(t f )-x f ( 1 ) 

PERFORMANCE INDICES 


Jj=(l/2) f (x T Qx ♦u T Ru)dt (2) 

0 

J 2 = J (1)dt= tj |U|U Ujjj, i=l ... n (3) 

0 

NECESSARY CONDITIONS 

H | =( 1 /2)(x T Qx ♦ U T Rll) ♦ A T (f(x) * Bu) (4) 

A = - (dHj/dx), A (0) unknown (5) 

OH 1 /^u)=0, Ru=-B T A (6) 

H 2 = I ♦ A T (f(x) ♦ Bu) (7) 

* 

A = - (dH 2 /dx), A (0) unknown (8) 

up - u ib sign(B T A ) , i=1 ... n (9) 

TPBVP 

z = g(z) , z=[ x . A ] T = [ Z| , z 2 1 T (10) 

Zj(0). z,(t f ) known. 
z 2 (0) f z 2 (t f ) unknown. 
z 2 (0) to be determined. 
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Equations op Motion 


M(y,) £ = Fly,, &>) + Bly,) v , 

l) — Control {/actor- 

Ml* O’* 

"\ 

1 + ¥^2*1 ~ 2^fyn 0 CoS(f> r t-2^Si'n(p r S^nune-trtc 

?”z-t Mfl ccscp r - yn s srnf> r M 3 -M 3 s>h(f> r J 

El tit; Vl) ^ 

-2Bi T (M4+ rr w wf> r -M l qs,n? r ) + j T [n 2 ccSf r + aJ, »Jq f r ] j 
6 a ( M/j + -ni A ct,sf r - M,r]srnf t ) - 2BH^ j ■Mf> r +( j T iff J ppjtesf^ ~K'l 

&a,) = 


I X l <^if l tu l Sinf l X 2 u>Sf 2 + U i l.„f 2 Z } uaf }i -u 3 Sm<f> 


L * 


t <*$% 


^ u,Sf 2 

'tf'} <M $ < P 3 1 

J 

lobesg H a 

H 

2 ; 

H </ 

M s ^ y/ 

, *2 , K 

— GChsta*vtr 1na±rtces 


7T!z y 

*1 o j 

™a. _ 

COlSt/crFC 

jcctors 

1 

/ 

771 % ' 

— C-CK St&*%/£s 
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QUASI Lj NEAR I ZAT ION AL60RITHM 


(A) LINEAR DIFFERENTIAL EQUATION: 

Nonhomogeneous: z = Az ♦ B. z=lzj. z 2 J T , (11) 

Zj(O), Z](t f ) known, z 2 (0) to be determined 
Homogeneous: z = Az (12) 

(a) n solns. of (12) ♦ 1 particular soln. of (11) 

(b) n+1 particular solns. of (1 1 ) 


(B) NONLINEAR CASE: 


Linearized equation of (10): 

^ =Og/dz) z^ k * * ^ ♦ h( z (k) ) (13) 

where 

z^ is the k^ approximate solution 

of the nonlinear equation (10), 
z (k+ 1 ) =z (k) +AZ (k) 

z-lz|, z 2 1 t , 

z j 1 ^(0), z 1 *(tf ), known 
z 2 (k+ 1 ) (o) to be determined 
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PLANAR SLEWING OF FLEXIBLE SCOLE 


LINEARIZED EQUATION OF MOTION: 


1 « T ' 


0 


0 

e T 

e 


i 

z \ 


L 

u s 




*#■ 




= 





U 1 

m II 


7 


o 

* 

X. 

- 

L \ 


o 

A 

*2 

& 

J 

u 2 


l U 3 


where 

6 is the angle of rotation, 

7 nX | is the amplitude vector of the flexible modes, 

n is the number of mode used, 

1 is the moment of inertia about the axis of rotation 
m. M are the Inertia parameter vector, matrix. 

K is the stiffness matrix, 

<£>(z) is the mode shape function vector, 

0j=^(^j) , zj is the coordinate along z axis, 

L is the length of the beam, 
u s is the control torque on the Shuttle, 

Uj are the control actuators on the beam and the 
reflector. 
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STATE EQUATIONS 


s - As ♦ Bu 


s = 


S 1 


e 


, si = 


s 2 


'I 


BOUNDARY CONDITIONS FOR s 


s(0)= 


e< 


o 

o 


s(t f )= 


0 

0 

0 

0 


2(n+1 


where n is the number of mode shapes used. 
PERFORMANCE INDEX 


J =(1/2) J (x t Qx +u T Ru)dt 
0 


TPBVP 

Z = CZ . z=I s . A ] T = ( Z|. 

A is the costate vector, 
Z|(0), Z|(t^) known; 
z 2 (0) to be determined. 
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. 16E+02 
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“MODE LINEARIZED CASE: U=( 1000), Q=DIAG{ 


. 10E+01 



-MODE LINEARIZED CASE: U*( 1 OOO ). 0-DIAG( 


M00E2 - ( FT ) 
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-MODE LINEARIZED CASE: U*( 1000), 0°D1AG( 



. 10E+04 
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-MODE LINEARIZED CASE: U=( 1000). Q=DIAG( 



MODE 1- (FT) 



3-MODE NONLINEAR CASE: U«< 1000). Q*OIAG< 


MODE 2 - ( FT ) 
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3-MODE NONLINEAR CASE: U=( 1000). Q=DIAG( 



. 10E+04 
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-MODE NONLINEAR CASE: U=( 1 OOO ), Q=DIAG( 



MODE 1 - ( FT ) 



-mode linearized CASE: U=( 1004 ), q-diag(0, 0, 0, 0), r-diag( 1 . 



MODE 2 * ( FT ) 



-mode linearized CASE: 



. 13E+04 



3/0 


mode linearized CASE: u-( 1 004 ). q=diag(0. 



. 12E+01 
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-mode linearized CASE: U - ( 1004 ), q=dlag(0, 0, 0. 0), r«d1ag( 1 . 


MODE 1- (FT) 



312 


-mode linearized case: u( 1 0 0 4). q=diag(0. 10. .1. .1, .1, .1, . 1 , .1), r=diag( 1 . 



M0DE2- ( FT ) 
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mod. linearized case: u ( 1 0 0 4 ). q.dlag(0. 10 





0.27E+01 



-mode linearized case: u( 1 0 0 4), q=diao(0, 10, .1, .1, .1, .1. .1. .1). r s diag( 1 . 


Concluding Remarks 


1) For the examples given here, the solutions of the 
linearized equation and the nonlinear equation are close 

2) Use of the Maximum Principle can make the states satisfy 
the boundary conditions very well. 

3) Due to the fact that the costates must be used in the 
method, the number of equations of the system is doubled 
and more time is needed for the computation. 

4) Further work on more complicated models (3-D) is needed. 

5) Need to consider different cost functions and perform 
parametric studies . 
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MOTIVATION 


• Modeling, analysis, and control of multi-body systems 
with flexible elements is of interest in many techologies 
including robotics, space structures, and high speed 
mechanisms. 

• While the basic problem is well studied and 
understood, there are some important aspects requiring further 
research. These include full nonlinear coupling, prismatic joints, 
controller design, and impact. 

• In all the research areas cited above, there is a need for 
experimental as well as analytical and numerical studies. 
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University of Michigan, Ann Arbor, Michigan 

316 


BACKGROUND 


Modeling: 

• Song and Haug 1980 

• Shabana and Wehage 1983 

• Ryan 1985 

Analysis and Simulation: 

• ADAMS, DADS, TREETOPS, etc 

• Sunada and Dubowsky 1983 

• Wang and Wei 1987 

• Wehage and Haug 1982; Khulief and Shabana 1984 

Control: 

• Book, Maizzo-Neto and Whitney 1975 

• Meckl and Seering 1983 

Experiments: 

• Zalucky and Hardt 1982 

• Cannon and Schmitz 1983 
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RECENT RESEARCH AT MICHIGAN 


• Control of a spherical coordinate robot arm with one 
flexible link using only the joint actuators, but with both joint 
and end-of-arm sensors. Modeling, analysis, and controller 
design. Modeling of the rigid/flexible manipulator including the 
non-backdrivable leadscrew transmission mechanisms. 
Comparison of simulation and experimental results. Evaluation 
of a rigid body motion controller and a rigid and flexible motion 
controller. 

• Modeling and simulation of robots with rigid and/or 
flexible links and revolute and/or prismatic joints. Employs a 
Lagrangian formulation with kinematic constraints and a finite 
element discretization. Euler-Bemoulli beam theory is employed, 
but the axial stiffening effect is included. Experimental 
evaluation of the modeling approach and solution method. 

• Modeling of impact in systems with flexible links. 
Various impact modeling methods are compared with each other 
and with the results of experiments for a radially rotating elastic 
beam. The impact models employed all provide reasonably 
accurate results when appropriately employed; even the 
momentum balance (coefficient of restitution) method which 
strictly speaking is not applicable to flexible systems. 


M Mechanical Engineering and Applied Mechanics 

University of Michigan, Ann Arbor, Michigan 

320 


FLEXIBLE MANIPULATOR CONTROL 


Lead Screw 



Schematic of the experimental setup. 
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FLEXIBLE MANIPULATOR CONTROL 


• The laboratory robot is used to compare 
the performance of a rigid body motion 
controller with that of a rigid and flexible 
motion controller. 

• The rigid body motion controller uses 
only the joint motion measurements and joint 
actuators. The rigid and flexible motion 
controller also uses the end of arm motion 
measurements, hut no additional actuators. 

• The leadscrew transmission 
characteristics as well as observation and 
control spillover are considered. 

• The numerical and experimental results 
show good agreement, and indicate that 
significant reductions in arm vibration are 
possible through use of the rigid and flexible 
motion controller. 
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FLEXIBLE MANIPULATOR CONTROL 
Equations of motion: 

+ E(* , i) = f'(T) 

— ~ ^ ) f// > ; ‘frzt > 7 

I T = C-r, ^ 7 

Linearized equations: 

■* ^ y + E a 

f=Lk T *i T ]> a- Sr 
S* T - Zir,Se t S*,S iH ,S itt 2 

Integral plus state feedback controller: 

ir t 

; y«= jC«fi-Rx)dt ; y iz -- 

a - - y 
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FLEXIBLE MANIPULATOR CONTROL 



Block diagram of the Integral plus state feedback controller. 
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FLEXIBLE MANIPULATOR CONTROL 


Standard Set of Physical System 
Parameters 

VALUE 

Mass of the first beam (m,j 

0.454 Kg 

Mass of the second beam (m.) 

0.816 Kg 

Mass of the Payload (m, ) 

Cross sectional area of the second 

0.07 Kg 

beam (A 5 ) 

0.000151 m : 

Length of the first beam (Lj) 

0.233 m 

Length of the second beam (L 7 ) 

2 ra 

Gravitational acceleration (g ) 

9.81 ni/tec 2 

Aluminum density (p) 

2707 Kg/ m > 

Flexural rigidity (El) 

770.87 Pa 

Reference position for r 

1.85 m 

Reference position for 9 

0 rad 

Reference position for 4 

0 rad 

Desired reference position for r 

2m 

Desired reference position for 9 

0.5 rad 

Desired reference position for 4 

0.5 rad 

Servo natarai frequency for r (w„ ) 

4 rad /sec 

Servo natarai frequency for # (w, ,) 

4 rad /Sec 

Servo natarai frequency for 4 (u, ,) 

8 rad/sec 

Flexible motion gain, K £ 

•0.000178 

Flexible motion gala, K 3 

-0.084 

Flexibk motioa gaia, K{ M 

1.568 
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FLEXIBLE MANIPULATOR CONTROL 



TVn^t iMBOOtf 

9 response obtained from the rigid body controller In 
the experimental work. 


f 



16.00 24.00 


QaaJ 

Tottl vertical deflection In responae to the rigid 
body controller In the experimental work. 
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FLEXIBLE MANIPULATOR CONTROL 



T>na_ f IttQttl 

• response obtained from the rigid and flexible 
motion controller In the experimental work. 



TOMA iMBQOft 

Control signal for the second Joint obtained from the 
rigid and flexible motion controller In the 
experimental work. 
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FLEXIBLE MANIPULATOR CONTROL 



Total vortical deflection In response to the rigid and 
flexible motion controller In the experimental work. 




maximum deflection 
(peak to peak) 

rigid body controller 

11.0 

7.5mm 

rigid and flexible 
motion controller 

3.0 

2.7mm 
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FLEXIBLE M VNIPULATOR CONTROL 
Simulation 

* Contro! spillover effect can be observed, but 
does not cause significant deterioration. 

Control and observation spillover can 

destabilize the residual mode. However, a 

small amount of damping (0.0145) eliminates 
the problem. 

* Settling time is reduced from 3.5 to 1 .07 

seconds, and maximum vibration amplitude is 
reduced by 5058. 

Experiment 

* With low pass filtering and light structural 
damping, no detrimental spillover effects were 
observed. 

* Settling time is reduced from 1 1 to 3 
seconds, and maximum vibration amplitude is 
reduced by 75 ?•. 
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FLEXIBLE SYSTEMS WITH PRISMATIC JOINTS 


• Robots with both rigid and flexible 
links attached with revolute and/or 
prismatic joints can be modeled and 

analyzed. , . . 

• The equations of motion are derived 

using Lagrange's equations. The prescribed 
motion, and prescribed torque/force cases 
can both be handled. 

. Flexible elements are represented as 

Euler-Bernoulli beams, and the axial 
shortening effect is also included. 

. Finite element analysis is used tor 

the discretization of the resulting hybrid 

equations of motion. 

• Constraints are handled using 

Lagrange multipliers. 

. The resulting algebraic-differential 

eauations are solved numerically using 
constraint stabilization methods. 
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FLEXIBLE SYSTEMS WITH PRISMATIC JOINTS 



Axial shortening of a beam under plane transverse deflection. 



Schematic of revolute joint i. 
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FLEXIBLE SYSTEMS WITH PRISMATIC JOINTS 




Schematic of prismatic joint L 
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l ip displacement (m ) 


FLEXIBLE SYSTEMS WITH PRISMATIC JOINTS 




M Mechanical Engineering and Applied Mechanics 

University of Michigan, Ann Arbor, Michigan 

33 + 


Tip displacement (m.) Tip displacement (in.) 


FLEXIBLE SYSTEMS WITH PRISMATIC JOINTS 




Tip displacement in "fast pull" case with 
Ci = 0.025 m, C 2 = -0.7 m and T = 0.7 sec. 
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Tip displacement (m ) Tip dis|ulccmeiil (m.) 


FLEXIBLE SYSTEMS WITH PRISMATIC JOINTS 



Vertical elastic tip displacement of the last link in the two-dimensional 
maneuver. 
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Elastic lip displacement (m ) 


FLEXIBLE SYSTEMS WITH PRISMATIC JOINTS 



A. Galip Ulsoy 

Mechanical Engineering and Applied Mechanics 
University of Michigan, Ann Arbor, Michigan 


Q| 




FLEXIBLE SYSTEMS WITH PRISMATIC JOINTS 


• A general modeling procedure for robot 
arms consisting of rigid and flexible links 
connected by revolute and/or prismatic joints 
has been developed and experimentally 
validated. 

• The significance of full coupling (effect 
of flexible motion on rigid body motion) has 
been demonstrated. 

• The axial shortening effect is shown to be 
significant for high speed operation of 
lightweight manipulators. 
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FLEXIBLE SYSTEMS WITH IMPACT 


The effect of impact on mechanical systems 

Impact 


• gives rise to impulsive forces 


These impulsive forces in turn 


• Induce high stress levels at different joints 
damage to the mechanical components of the system 

• Cause higher modes of vibration to be excited 
deviation from a desired performance 
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FLEXIBLE SYSTEMS WITH IMPACT 


Spring-Dashpot Models 

Dubowsky and Freudenstein (1971) Impact pair model 
Crossley and Hunt (1975) Viscous damping model 

Lee and Wang (1982) Damping functions 

Dubowsky and Gardner (1975) Impact beam model 

• Based on a force-displacement law, and a form of damping 

• Does not neglect the contact duration 

• It is possible to predict the contact forces directly 

• Requires the determination of stiffness and damping 

• Computationally expensive 
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FLEXIBLE SYSTEMS WITH IMPACT 


Methods based on Momentum Balance 

(Coefficient of Restitution Model) 


Haug and Wehage (1981) For rigid body impact 
Khulief and Sfiabana (1984) Extension to flexible bodies 


• Simplest model for impact 

• Computationally efficient 

• It is possible to predict the contact forces directly 

• Does not require determination of stiffness and damping pa- 
rameters 

• Neglects the period of contact 

• A direct prediction of stresses and frequency content of impulse 
during contact is not possible 

• Mathematically discontinuous 
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FLEXIBLE SYSTEMS WITH IMPACT 


Procedure for Nlomentum Balance Method. 

1- Natural frequencies and mode shapes are determoned for 
each flexible body in the system; analytically, experimentally, 
or by using a finite element approximation. 


2- Discretized equations of motion are generated and integrated 
forward in time using this modal information. 

3- Impact conditions are checked, if an impact is detected to 
occur, integration is stopped, and the equations of momentum 
balance are generated using the coefficient of restitution and 
solved for jump discontinuities in velocities. 


4- The velocity vector is updated 

5- Integration is started once again with the new initial condi- 
tions found from above. 

(It is assumed that the system configuration does not change 
with impact.) 
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FLEXIBLE SYSTEMS WITH IMPACT 


Coefficient of Restitution, e 

Newtonian Concept: 
e = f (material pair) 


Modern Concept: 
e = f (severity of the impact) 


for rigid bodies 

severity of impact = impact velocity 

e = f (material pair. impact vel.) 

experimental data or analytical/emprical formulas 

for flexible svstems 

— x. 

the severity of impact = f (impact vel. .flexibility, configuration) 

No experimental data or analytical formula available 
for e 
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FLEXIBLE SYSTEMS WITH IMPACT 
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FLEXIBLE SYSTEMS WITH IMPACT 


Equations of Motion 

— Euler- Bernoulli Beam Theory 

— Hamilton’s Principle 

— Neglect higher order elastic terms 

— Galerkin's Method 

w(x,t) = 

j = i 

— Neglect higher order elastic terms 


N 


(Ir + J)0 + E SjQj = M(t) 

j = 1 


N 


S r 9 + E \jn r j cjj + 9 2 qj(c r j — m r j ) 4- k r jqj] = 0 
i=i 


where 

Ir : Inertia of the Rigid Shaft 
J i Inertia of the Flexible Be am 
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FLEXIBLE SYSTEMS WITH IMPACT 


Simulations 

- A variable order variable time step integrator 

- A velocity dependent coefficient of restitution 

e = f(v) (experimental data) for low velocities 
e = (analytical data) for high velocities 

Contact algorithm 

based on the distance between impact point and the impact 
surface 

— define a clearance zone 

— monitor the location and the velocity of the impact point 

— first penetration into the clearance zone — ► back up one time 
step and reduce the time step. 

— second penetration with the smaller time step — > impact if 
the velocity is toward the surface 

— keep the smaller time step size as long as the impact point 
is within the clearance zone 
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FLEXIBLE SYSTEMS WITH IMPACT 
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FLEXIBLE SYSTEMS WITH IMPACT 

The momentum balance (coefficient of restitution) method 
has been demonstrated, through experiments on a specific sys- 
tem, to be capable of adequately predicting the dynamic be- 
havior, with impact, of systems which consist of both rigid and 

elastic links. 

It has been demonstrated experimentally that using a con- 
stant coefficient of restitution value (chosen for a particular 
initial impact velocity) throughout the simulation does not sig- 
nificantly affect the accuracy of the model, even in the presencof 
multiple impacts. 

Sensitivity studies were used to show that the momentum 
balance (coefficient of restitution) model will also work reason- 
ably well for a wide range of system parameters. 


Employing high speed video techniques the existence of mul- 
tiple impacts which appear to the naked eye as a single contact 
were demonstrated. Contact algorithms have been developed 
which captures these multiple impacts in the simulation. 
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SUMMARY AND CONCLUSIONS 


• This presentation has reviewed the results of some recent 
studies, experimental and theoretical, in the modeling, analysis, 
and control of flexible multi-body systems. The control of a 
leadscrew driven robot with a flexible link was considered. The 
modeling of robots with rigid and/or flexible links connected by 
revolute and/or prismatic joints was developed. Studies were 
conducted to evaluate various competing impact models for use 
in simulations of flexible systems. 

• The major conclusions to be drawn from these studies are: 

For robots with sufficient actuator bandwidth, the use of end- 
of-arm sensors can enable effective control of the joint motions 
as well as active damping of the end-of-arm vibration. 

A general modeling and simulation approach for robots 
consisting of rigid and/or flexible elements connected by 
revolute and/or prismatic joints has been developed and 
experimentally validated. 

Impact models for systems with flexible elements have various 
problem dependent characteristics, however, even the simplest 
of these can give good results in many problems of engineering 
interest. 
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MINIMUM ATTAINABLE RMS ATTITUDE ERROR 
USING CO-LOCATED RATE SENSORS 


A. V. Balakrishnant 


Abstract 


In this paper we announce a closed form analytical expression for the minimum attain- 
able attitude error (as well as the error rate) in a flexible beam by feedback control using 
co-located rate sensors. For simplicity, we consider a beam clamped at one end with an 
offset mass (antenna) at the other end where the controls and sensors are located. Both 
control moment generators and force actuators are provided. The results apply to any beam- 
like lattice-type truss, and provide the kind of performance criteria needed under CSI 
Controls-Structures-Integrated optimization. 


t Research Supported in part under NAS1-18585 Task Assignment 49. 


Paper presented at the 3rd Annual Conference on Aerospace Computational Control, Oxnard, 
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1. Introduction 


One of the challenges in the Design Challenge For Flexible Flight Structure Control 
System Design formulated in the inaugural paper on SCOLE [1] was to hold the antenna 
pointing error within ±0.02 degrees after slewing by appropriate feedback control. In this 
paper we derive a closed form expression for the minimal achievable mean square pointing 
error using co-located rate sensors. A slightly simplified form of the SCOLE article (which 
eliminates rigid-body modes) is used: a cantilevered beam with an offset mass where the 
controls — both c.m.g.'s and force actuators — and the rate sensors are located. Our results 
are in terms of continuum model parameters — the uniform Bernoulli version is used. The 
beam dynamics are given in Section 2. The main results are in Section 3. We note that a 
technique for deriving equivalent Bernoulli beam parameters for various types of trusses is 
described by Noor and Anderson in [4]. Recently Noor and Russell [5] presented equivalent 
anisotropic Timoshenko beam models for beam-like lattice trusses with an arbitrary degree of 
modal coupling, which appear to yield excellent agreement with modal frequencies derived 
from finite element models. Our theory is able to handle these Timoshenko models, and 
moreover we can also use it for rigid-body modes, although they are not included here. Thus 
our results can be used for any beam-like lattice truss structure. 



2. The Model 


We consider a uniform Bernoulli beam clamped at one end with an offset mass (antenna) 
at the other end which also houses the sensors and actuators. See Figure 1. We allow for both 
force actuators and moment actuators. The sensors are rate gyros. Because of the clamping at 
one end, no rigid-body modes are involved and hence no attitude sensors are needed. 

We allow bending in two mutually perpendicular planes containing the beam axis, as 
well as torsion in the plane perpendicular to the beam axis, all uncoupled. The continuum 
model (uniform Beroulli beam) dynamics can then be described by the following partial differ- 
ential equations (similar to those in [2, 3]). Let the beam extend along the z-axis, 0 < s < L, 
and let u 0 (s, t), u e (s, t), denote the bending displacements and u^(s, t) the torsion angle 
about the beam axis. Let in the usual notation (cf. [1]), El 0 , E/ 0 denote the flexural stiffness 
and G/ v the torsional rigidity. Let p denote the mass per unit area and A the cross-sectional 
area. Then we have: 


pA 


pA 


t) 

dt 2 

3 2 m 0 (j, t ) 

dt 2 


+ El, 


d\( s , t ) 
♦ ds* 


+ El, 


d*u Q (s, t) 
& Bs* 


= 0, 


0, 


9 2 « U) (5, t) 

p/y — — - GV;(S, 0 = 0, 
with the clamped boundary conditions at 5 = 0: 


0<^<L; 0<f 


0<s<L; 0<f 


0 < s < L \ 0<f 


u 0 (O, r) = u 0 (O, t) = Uy(0,t) = 0 


«*(0, t) = Uq (0, t) = 0. 


The antenna center of gravity is located at 


(r x , r y ,L) . 

The distance from the beam tip to antenna center of gravity is denoted by 

\r\ = Vr? + r 2 . 




Figure 1 
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: Shuttle/Antenna Configuration 
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The force balance equation at s = L yields 



1 0 r x 


M* (L, t) 


m 


El^u;\L, r) 

m 

0 1 r v 


Uq (L, t) 

4 

fi(t) 


EIqUqXL, t) 


y 


My (E> 0 





where m is the antenna mass and /i(-), hi') are the applied control forces. The torque 
balance equations yield 


0 = 

EI^i'(L,t) 
EIq u&(L, t) 

+ l a (d + M(t) + r® 

flit) 

flit) 


G/yM^L, t) 


0 


u^(L, t) 4- r x tiy (L, t) 


where the superdots indicate time derivatives and the primes the derivatives with respect to 
the spatial variable s; ® denotes the vector cross-product and (0 the angular rate vector 


(0 


“o' ( L - r ) 

Mq (L, f) , 

My(^>0 


i a denotes the moment of inertia of the antenna about the beam tip (s - L) and finally, M{t) 


denotes the applied control moment. 

It is convenient to denote by b{t) the boundary vector 


b(t) 


u 0 (L, t) 
u;(L,t ) 
UQ (L, r) 

My (E, 0 


The boundary rate vector would thus be b(t ) . Hence our sensor model is: 

v(f) = b(t) + N 0 (t) 

where we assume that N 0 (t) is white Gaussian noise with spectral density matrix dj, where 
/ is the identity (5x5) matrix. Similarly we assume that the control actuators are also 
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characterized by additive white Gaussian noise. Denoting the applied control vector by u(r): 


u(t) 


«i (L,t) 
t) 

ui(L, t) 
u 4 (L, t) 
u s (L,t) 


we have 


/ 1(0 
/z(0 

a m 


u(t) + NJt) 


where N s (t) is white Gaussian with spectral density d s . We shall also use M b to denote the 
actuator mass/inertia matrix 

I 

m 0 0 

0 m 0 

0 0 

0 0 


where 


M b = 


0 mr. 


mr„ 


mr x mr y 




~ r x r y 


-r x r y 


0 

0 


where ! a is the antenna moment of inertia about its center of gravity. For any control input 

u( ) (which must perforce be a “feedback” control, based on the sensor data v( )) the mean 

square pointing error is then expressed by: 

• f T T t 

1 im = \ J Uq(L, t) 2 dt + J u e (L, t) 2 dt + |r| 2 f u (L, t) 2 dt 
T -*°° to 0 0 


and the mean square pointing rate is given by 



{ T T T 

J «<,(£, t) 2 dt + J u 0 (L, t) 2 dt + \r\ 2 J u^(L, t) 2 dt ■ . 

0 0 0 

From the results in [6] it follows that the minimal attainable mean square pointing error is 
given by 

aM~ b 1 a* 

where 

a = row vector ( 1 , 1,0,0, |r|) 
a* = transpose of a 


3. Main Results 

We need some notation first. The mean square attitude response, whatever the feed- 


back control used is defined by 


1 im 5; { f u 0 (f, l) 2 dt + J u 6 (t, l) 2 dt + \r\ 2 J « v (r, l) 2 dt ' . 
1 l o 0 0 ^ 


(3.1) 


This is recognized as the mean square displacement of the center of gravity of the antenna 
which is then also proportional to the mean square “pointing” error — see [1] for the 


relationships. 

Next let u denote any (vector) of control inputs — a constant “step” input: 


u 


«i 

m 2 

«3 

«4 

“5 


(3.1) 


Solve the equations 


EW'(s) = 0 ' 
£/ e «e”(s) = 0 
GIyu;(s) = 0 , 


0 < s < L , 


(3.2) 



subject to the end conditions 


EI^u;\L) = Ul 
EWiL) = u 2 
El $ Uq(L) + «3 = 0 
EIqUq(L) + u 4 - 0 
GlyUy(L) + U 5 = 0 ■ 


(3.3) 


Note that the solution can be recognized as the steady-state response of the system to the 
step-input u , assuming that there is some damping. We only need to calculate the response to 
three specialized inputs: 

Calculate the response to Wp(L) to the special case. Case 1, where: 


«i = 1 


u, = 0 , 2 < i < 5 . 

Calculate the response u e (L) to the special case. Case 2, where: 

«i = 0 

«2 = 1 

«3 = u 4 = « 5 = 0 . 

Calculate the response ( L ) to the special case, Case 3, where 

«i = u 2 = K3 = « 4 *= 0 


«6 = 1 . 

Then the minimal achievable mean-square response whatever the choice of the feedback and 
whatever the mean-square control effort, is given by 

(«0 ( L > 2 + “e (E) 2 + ^«y(^) 2 ) • (3.4) 

This is our main result. Unfortunately the derivation is beyond the scope of this report and 


364 



will be published elsewhere. To proceed further with (3.4) we calculate the solution of (3.2), 
(3.3) explicitly. Thus for any u, the solution is of the form 


where 


Thus for Case 1 we have 


and for Case 2 we have 


and for Case 3: 


(■*) 


a 3 




a 2 




C 1 


- a 3 s 3 + a 2 s 2 , 
= 63 j 3 + t^s 1 , 


C\S , 



U \ 





U 2 


6£/ e 


1 


U\t' 

2 

i EI * + 

A. 

1 


U\L‘ 

2 

[EIq + 

EIq, 


_ i_ 

" G/ v • 



« e (^) 2 


L 3 

3£/ e 


S ,^) 2 = 


L 

G/ v ' 


Hence the mean-square attitude error 


0 < s < L 
0 < s < L 
0 < s < L 



JLL . \r\ 2 L] 

3 EI q G/ v J • 


(3-5) 


Note the appearance of the noise parameters in (3.5) in product form. 

The technique for calculating the minimal mean square atttiude error in more complex 
models than that illustrated is the same: calculate the mean square step response (assuming 
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some damping) to unit step inputs. 

In conclusion we suggest this result (3.5) can be the basis for combined structures- 
controls optimization — CSI, since the required structural parameters can be calculated for a 
lattice truss from the material gage and physical dimensions as in [4, 5]. We omit the details 


of these calculations. 
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A WAVE EQUATION FORMULATION FOR THE SCOLE 

Roberto Araya 
AutoMi nd 

Providencia 591 D.51 

Santiago Chile. 


This paper review the modeling and control strategy proposed by 
the author [1][2] -for the slewing and stabilization problem o-f the 
Space Shuttle Orbiter with a large Antenna attached to it 
through a long flexible mast. 

The main distinguished feature of the problem treated here as com- 
pared with the standard problems studied in this area is the require- 
ment to find optimal solutions for a nonrigid spacecraft configu- 
ration which is subject to nonlinear kinematic forces. 

A distributed parameter model of the Space Shutt 1 e/Antenna Confi 
guration is derived from first principles of rigid body dynamics and 
elementary beam theory. The model is then put in the form of a compact 
semi linear abstract wave equation. Within this framework, the slewing 
and stabilization problem is then formulated as a nonlinear infinite- 
di mensi onal control pr obi em. 

A linear feedback control law is proposed to simultaneously solve 
the slewing and stabilization problem. This control law is a generali- 
zation of standard position-plus-rate feedback controls used for simi- 
lar purposes in rigid spacecrafts and is also an extension of linear 
feedback controls used in stabilization theory of linear infinite di- 
mensi onal systems. 
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STATEMENT OF THE PROBLEM 


Difficulties in controlling lightly damped Large Space Structures 
has motivated NASA to offer a design challenge called Spacecraft 
Control Laboratory Experiment <SCDLE)[5] . The aim of this "design 
challenge" is to evaluate control laws for flexible spacecraft. 

The SCOLE configuration consists of a flexible beam with two rigid 
bodies at either end. One body represents the Space Shuttle orbiter, 
the other is the Antenna reflector. 

I here are two types of control : control moments and the control 
force. Control moments denoted by M, , are applied to the Shuttle 
orbiter and the reflector, respectively. The components of these mo- 
ments for each axis are limited to lU,o00 ft— 'lb. Control force F is 
applied at the center of the reflector in two perpendicular directions 
along the plane of the reflector. The control force in a particular 
direction is limited to 80 h lbs. 

Measurements consist of the inertial attitude direction cosine 
fPcit.f ices for the Shuttle body and the Antenna, angular velocities o , 
for the Shuttle and the reflector body respectively and accelera- 
tions measured by three axis accelerometers located on the Shuttle and 
the center of the reflector. 


A certain r <ny is emitted towards the center 
the feed located on the Shuttle. The direction 
ray is called the 1 i nc-cf-si ght . Initially the 

1 i ne —of — si ght and the target is 20 degrees and 
is at rest. 


of the reflector from 
of reflection of this 
error between the 
the SCDLE configuration 


The problem is to minimise the time required to slew or change the 
1 i ne— of — si ght towards a fixed target and to settle or stabilise the 
induced structural vibrations of the Shuttle/Antenna configuration to 
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the degree required for precise pointing. 


THE EQUATIONS OF MOTION AND THE CONTROL SATRATEGY 


of 


The derivation of the equations of motion is based on 
rigid body dynamics and elementary beam theory. 


basic principles 


Under the assumption that the flexible beam is subject to only 
small deformations, a linear Eul er /Bernoul 1 i beam model is utilised. 
Nonlinearities are due only to kinematic effects caused by the large 
angular velocities that the whole configuration can reach in a time 


OD '■ 


t i mal maneuver . 


S i n c e i n h e r © n t s *t r u c t u r a 1 
i s n o t. w 0 1 1 l< n o w n , n o d a in p i n g 
dissipation is caused only by 


damping is very weak and the phenomenon 
is assumed in the model. Therefore energy 
the control action. 


G:i. ven that 
on 1 y r ot a t i on a 1 


the main concern is spacecraft 
dynamics is considered- 


ori entati on , 


A coupled system ot hyperbolic partial differential equations and 
ordinary differential equations is obtained. The variables considered 
are the elastic deformations u* and u* , the torsion u* and the 
longitudinal deformation u 2 along the beam axis. 


To obtain these equations, the 
configuration about the center of 
momentum p^ ( p^) and the angular 
the beam ( Antenna ) are computed. 


angular momentum p ( of the whole 
gravity of the Shuttle and the 
momentum p (p 1 of each element 

3 5 


of 


After deriving the equations of motion, 
equation can be obtained: 


the following Lyapunov 
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d_ [D(t)+T(t)+V<t> ] 
dt 


« 

= e* (IW X jDj 8Rj ) +•,, ♦IM^+v., «F 

i= i 


where 


T is the kinetic energy, 

V is the potential energy 

D is a measure of the distance to the target orientation 

D - h E x i°i • D i 

i ~ i 

D i =IR ^ — S j i = i , . . , j 

(Si , i” is an orthonormal triplet in the inertial -frame 

(IRi , i= is an orthonormal triplet fixed to the Shuttle 

» angular velocity of the shuttle 
e, is the angular velocity of the antenna 

v, is the linear velocity of the mass center of the antenna 

This equation provides a clear insight into our problem . It 
suggests a strai ghtf oward and i mpl ementable feedback control law : 
choose Hull, and F so that the right hand side in becomes strictly 
negative . 

Thus, the problem is broken down into two parts 

i) Orientation of the whole configuration as a rigid spacecraft. 

This is achieved by choosing ID, to make XiBjQRj) negative. 

X *= I 

ii) Stabi 1 i zati on of the induced structural deformations 
Choosing M, and F to make o, *IHI. ( +v >1 *F negative 


The first control strategy is well known for rigid spacecreaft 
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stabilisation and slewing. The second one, is the classical 
stabi 1 i zati on control for rigid and flexible structures. 


THE ABSTRACT WAVE EQUATION SET UP 


The equations of motion can be written as an abstract semi linear 
wave equation: 


x (t) = Ax ( t ) + 


(t ) ) + K, (x (t) ) + Bu(t) 


x(0) = x, 

in a Hilbert space R . Here, A is an infinitesimal generator of a group 
with compact resolvent, B is a linear operator with finite dimensional 
range, K, and K z are nonlinear but conti nous operators : 

K, with finite dimensional range and representing the orientation of 
the configuration as a rigid spacecraft; 

K z represents the nonlinear kinematics and satisfies 

[ K z < x ) , x ] = 0 Vxe*. 


More specifically, 


H = D ( A 1 ' r ) x D ( A • ) 

0 0 


W 


ith the inner product 


[ 1 

} f y * ) 

[ X S 

V: 1 


C A 


/ Z v. A I / * v 1 

, " • ’ H o y • j 


[Hx, , y x ] 


i/,! h fp r e 
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D ( A ® ) = R'x ( L z ( -L , 0 ) )’x*' 




X,D, 

B z 


X«»i 

U. 8 


eiJLV 

a? 1 * 

U * 


Ell-V 

as" 

u- 


-eaII 1 -' 2 
as 1 

u. * 


-GllV 

as 1 

r 1 

55 

Ei a’u 8 (-L) 

as 5 

r : 


ex a 3 u 8 ( — L ) 

as 5 

r , 


— EA 8t -' 2 ( ~ L ) 

as 

® i 


ei i:u*(-d 

as 1 

« z 


_ EI a *u 8 ( -l ) 
as 1 

a 3 


_ei au ♦ <- l) 


as 


and the domain cf A 0 is 

D(A„) = ( (I) , , U, , u. 8 , u * , li z , u * , r , , r z , r , , et , , c z , c , ) eD ( A jj ) such 
u 8 ,U*E H ’ (-L,0) , u z ,u*eH* <-L,0> 

u 8 (o) =u ♦ (0) =u s (0 )=li* (0) =iy 8 (0) = au ♦ <o)-o 

as a? 
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that 


u 8 ( -L ) =r , ,U*(-L )=r, ,u z (-L)=r, , 

lu 8 (-L) =a lu*(-L) =c ,u*(-L)=a, ) 


M is the mass moment operator, defined from D(AJ) into itself by. 


0 

€» 

O 


PA ( * > 0 V 
PI 


0 # O 

1 f 0 ( • ) PAd 5 
-L 


m , 


( « ) 0 ( H - + f ) 


I .,+m, P0 ( ( * ) 0|», ) 


PA 

O 

O 

o 


0 


o 


o 


Jpi + (*)zxU m,<P,+ F>0<*> I ,+m,p,0( ( *> ®P> 

•~L ... n 


(.) 

PI * 

e 

o 


o 
o 

m H 

m , f 0 ( • > 


O 

o 


m n ( • > 0 1 


In this abstract set up, the norm induced by the inner product is: 
<1/2)H x (t) »e - D ( t ) +T(t) + V<t) 


and therefore the slewing and stabilisation problem can be viewed as 
driving the initial state x 0 towards the origin. 

This compact and geometric formulation of the problem is very 
fruitful. It allows a detailed mathematical analysis of the evolution 
equation and the control strategy. This gives existence results 
of mild and strong solutions under suitable conditions. It gives also 

the correct framework to do modal analysis: 


A ,x = pMx 

This compact formulation gives also the insight to understand the 
proposed control laws as a natural general i sation from a finite 
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dimensional to an infinite dimensional set up. 

The proposed control stategy can be expressed as: 


) h!x : & standard linear feedback control used for slewing of 
rigid spacecraft. 

i i > : a standard feedback stabi 1 i zat i on control used in in- 

finite dimensional control theory and widely known 
for stabilization of flexible beams. 

Both components have simple physical interpretation and can 
be computed directly from the sensor data. 

The Lyapunov equation and the proposed control law, is expressed 
d i x ( f ) i * — 

= - »♦# ~ - v, *v, 

Thus, * x ( t ) i is decreasing. The original evolution equation is 
converted to the autonomous system: 

x = Ax + V ( x > 
x(0) » x, 

or, the integral equation version: 

P t 

x( t) - e^'t-x 0 + e A (t-s) y (x (s) ) ds 

J 0 

The major result obtained with this approach is the ability to 
show that if the system is initially at rest and the 1 i ne-of -si ght erro 
: : s whithin 20 degrees, then the linear feedback position-plus-rate 


control law : 



u = -B*x 


+ Rx 


or equivalently 

j 



will drive >x(t) * to zero as t goes to infinity- Thus, in finite time 
the 1 i ne-of -si ght will be pointed to the target up to an error j and 
simultaneously the oscillations of the flexible configuration will be 
settled down, in the strong sense, to the degree required to precise 

poi n t i n g - 

To achieve this result is necessary to prove the cont r ol 1 ab i 1 1 ty 
of the.' flexible modes* That is, when linearising and considering 
only flexible modes, the obtained linear model has to be controllable. 

This abstract wave equation formulation is a basic framework for 
the complete problem analysis. It also offers a correct and compact set 
up to study the different control strategies. In particular, non linear 
feedback controls on the angul ar and linear velocities [3 ] [4 1 [6 ] . 
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